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1.0  INTRODUCTION 


For  the  past  dglit  years,  the  Wright  Laboratory  Armament  Directorate  at  EgUn  Air 
Force  Base  has  coi^ucted  a  Lethality/Vuhierability  program  to  evaluate  the  effectiveness  of 
kinetic  energy  weapons  (KEWs)  against  ballistic  imssiles.  This  program,  part  of  the  Lethality 
and  Target  Hardening  (LTH-S)  Program  of  the  Strate^c  Defense  Initiative,  has  focused  on 
the  re^nse  of  ballistic  missile  boosters,  post-boost  vehicles,  and  their  associated  waiiieads  to 
KEW  impacts.  The  evaluation  and  selection  of  the  systems  to  advance  from  the  conceptual 
phase  of  design  to  demonstration/validation,  aigineering/manufacturing  development, 
production,  and  deployment  requires  the  assessment  of  candidate  weapons  effectiveness 
against  their  intended  threat  spectrum. 

The  response  of  a  target  to  a  KEW  impact  can  be  said  to  consist  of  two  basic  and 
distinct  types  of  response:  local  response*  and  'global  response*.  For  KEW  impacts,  material 
damage  assodated  with  local  response  occurs  very  quickly  O-e.  within  the  first  100-200  psec) 
and  is  linuted  to  a  volume  immediately  adjacent  to  the  impact  rite.  At  sufiBciently  high  impact 
vdodties,  shatter,  melting,  and/or  vaporization  of  the  materials  can  occur.  For  an  aluminum- 
on-aluminum  impact,  the  projectile  and  target  materials  will  begin  to  shatter,  melt,  and 
vaporize  at  impact  velodties  of  approx.  3.2,  S.6,  and  10.4  km/sec,  respectively  [1,2]. 

Gobal  response  can  refer  to  any  one  of  a  number  of  global  phenomena  that  occur  over 
a  longer  period  of  time  (on  the  order  of  milliseconds),  under  less  intense  loads,  and  over  a 
much  largo’  area  of  the  target.  In  KEW  impacts,  one  or  more  dd>ris  clouds  are  created  during 
the  initial  impact  on  the  outer  wall  of  a  target.  These  dd>ris  clouds  spread  out  as  they  move 
through  target  voids  and  eventually  impact  an  inner  wall  or  interior  component  of  the  target 
structure.  Depending  on  the  impact  velodty  and  the  relative  material  properties  of  the  projec¬ 
tile  and  target,  these  debris  clouds  can  contain  solid,  melted,  and  v^rized  projectile  and 
target  materials.  The  levels  of  melt  and  vaporization  within  the  d^ris  clouds  in  turn 
determine  the  nature  of  the  loads  transmitted  to  various  target  components.  Typical  global 
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realises  include  the  denting,  buckling,  tearing  or  catastrophic  dismemberment  of  internal 
misale  conq)onents. 

To  accurately  determine  total  target  damage,  a  lethality  assessment  methodology  (see, 
e.g.  [3-12])  must  include  the  effects  of  discrete  impacts  by  solid  d^ris  cloud  fragments  as 
wdl  as  impulave  loadings  due  to  molten  and  vaporous  debris  cloud  material.  Clearly,  the 
amount  of  debris  cloud  material  in  each  of  the  three  states  of  matter  must  be  known  to 
accurately  assess  total  target  damage  and  break-up  due  to  a  KEW  impact.  This  rq)ort 
presents  a  first-principles  method  to  calculate  1)  the  amount  of  material  in  a  dd>ris  cloud 
created  by  a  perforating  hypervdodty  impact  that  is  solid,  molten,  and  vaporous,  2)  the  ddiris 
cloud  leading  edge,  trailing  edge,  center-of-mass,  and  expanrion  velocities,  and  3)  the  angular 
spread  of  the  dd>ris  cloud  material.  The  method  presented  can  be  used  for  single-  and  multi¬ 
material  solid  rod  projectiles  impacting  an  array  of  target  plates.  The  methodology  presented 
in  this  report  includes,  improves,  and  expands  upon  the  debris  cloud  characterization  scheme 
presented  in  WL-TR-93-7028  [13].  As  such,  the  information  this  report  is  intended  to 
supo^sede  the  information  in  [13]  and,  with  the  exception  of  the  portions  that  have  been 
modified,  the  information  in  [13]  is  reproduced  in  its  entirety  in  this  report  for  completeness 
and  to  m^tain  continuity.  At  this  point,  no  adjustments  have  been  made  to  account  for 
diff^ences  in  response  due  to  projectile  yaw  or  impact  obliquity.  The  predictions  of  this 
methodology  are  compared  against  those  of  empirically-based  lethality  assessment  schemes  as 
well  as  against  numerical  and  empirical  results  obtained  in  prevdous  studies  of  hypervelocity 
impact  dd>ris  cloud  fonnation. 
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2.0  lj:THiiUTY  ASSESSMENT  MODEL  REQUIREMENTS 


The  key  to  conducting  an  accurate  lethality  assessment  is  the  use  of  a  robust 
assessment  methodology.  The  methodology  should  incorporate  all  the  agnificant  response 
and  damage  mechanisms  which  result  from  all  hypervelocity  weapon-target  interactions.  To 
accurately  determine  the  total  danuge  level  sustained  by  an  impacted  target,  a  lethality 
assessmoit  methodology  must  include  the  effects  of  discrete  and  simultaneous  dd>ris  cloud 
fragment  impacts,  as  well  as  impulave  target  dd>ris  doud  loadings.  Discrete  or  amultaneous 
impacts  by  individual  fragments  can  pose  a  lethal  threat  to  the  inner  wall  or  to  an  interior 
component  of  a  taiget,  dq)ending  on  the  fiagmoits'  q)eed,  density,  and  tnyectoiy,  and  on  the 
density  and  strength  of  the  target  inner  wall  or  interior  component  material.  Individually,  the 
molten  and/or  vl^)orous  fragments  in  a  debris  cloud  may  not  do  significant  damage;  however, 
as  a  vdiole,  they  can  produce  a  significant  impulnve  loading  over  a  relatively  large  area  imdde 
the  target.  This  in  turn  can  result  in  fiirth^  damage  to  the  target  at  later  times.  Clearly  then, 
to  accurately  assess  the  total  damage  to  a  target  impacted  by  a  KEW,  the  amounts  and  types 
of  dd)ris  in  a  dd)ris  cloud  produced  by  a  hypovelodty  impact  must  be  known. 

A  number  of  empirical  and  semi-analytical  procedures  have  been  developed  over  the 
past  decade  to  determine  the  lethal  effectiveness  of  KEW  systems.  While  these  procedures 
are  capable  of  assisting  engineers  and  system  architects  in  optimizing  weapon  deagns  and  in 
performing  cost  trade-  off  studies,  they  are  agnificantly  limited  in  their  characterization  of  the 
material  in  the  dd>ris  clouds  created  by  ttypervelodty  impacts.  Unforturutely,  very  little 
impact  test  data  for  relativdy  massive  projectiles  (on  the  order  of  10  gms  or  more)  is  available 
at  greeds  above  8  km/sec.  This  makes  it  difficult  to  property  characterize  the  nature  of  the 
mat«ial  in  the  dd>ris  clouds  over  the  entire  impact  velodty  regime  of  interest.  Electrostatic 
devices  which  can  launch  small  particles  to  speeds  as  high  as  100  km/sec  exist,  but  these 
tystems  can  only  launch  micron-^  particles  [14,15].  Other  electric  gun  systems  have 
launched  Ki^)ton  flyer  plates  to  speeds  of  1 1  km/sec,  but  cannot  reach  that  velocity  with 
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diunky  projectiles  [16].  Thus,  existing  lethality  assessmoit  models  must  be  used  with  a  &ir 
amount  of  caution,  especially  in  scenarios  involving  impact  vdocities  greater  than  those 
attainable  in  experiments. 

Currait  semi-analytical  lethality  assessment  models  .usually  fidl  into  one  of  two  broad 
groups:  discrete  partide  modds  [3-12,17,18]  and  expanding  shell  models  [19-23].  Discrete 
partide  modds  typically  account  for  only  solid  fragments  [3-8,17,18],  or  track  only  a  small 
numbo*  of  discrete  fragments  [9-12]  in  the  ddiris  doud  created  by  a  high  speed  impact. 

These  modds  are  best  suited  for  applications  in  which  the  debris  clouds  generated  by  the 
initial  impact  contain  only  a  rdatively  snail  number  of  fragments  and  in  which  melting  or 
vq)orization  of  the  projectile  and  taiget  materials  do  not  occur. 

The  expanding  shdl  modds  typically  assume  that  all  of  the  dd>ris  cloud  material  is 
homogeneously  distributed  ovor  a  uniformly  expanding  spherical  shell.  These  models  are 
^plicable  only  in  those  impact  situations  where  complete  projectile  and  target  materid 
vaporization  occurs.  Flash  X-ray  photographs  of  the  dd)ris  clouds  created  in  lead-on-lead 
impacts  at  speeds  high  enough  to  cause  melting  and  vaporization  do  show  that  the 
assumptions  of  the  spherical  sheU  model  are  valid  at  least  for  the  leading  portion  of  the  dd>ris 
clouds  [24].  However,  when  a  debris  cloud  is  comprised  primarily  of  solid  fragments,  then 
simOariy  obtained  photographs  show  the  ddrris  clouds  to  be  elliptical  with  an  eccentricity  of 
tqjproximatdy  1.6  [24]. 

It  is  evident,  therefore,  that  the  need  exists  to  bridge  the  gap  between  the  discrete 
particle  models,  which  consider  only  a  finite  number  of  solid  fiagments,  and  the  expanding 
shell  modds,  which  are  most  accurate  when  complete  vaporization  occurs.  Spedfically,  a 
lethality  assesnnent  modd  that  considers  the  creation  and  subsequent  effects  of  dd)ris  clouds 
containing  all  three  matter  states  is  needed.  FATEPEN2  [3-6],  PEN-4  [7,8],  and  KAPP-n 
[10-12]  are  discrete  particle  lethality  assessment  models  which  can  be  modified  to  include  the 
effects  of  non-solid  ddrris  cloud  constituents. 
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The  FATE  ftmily  of  codes  was  developed  for  the  Naval  Sur&ce  Warftre  Qeiitn' 
(NSWC)  for  analyzing  the  impacts  of  warhead  fragments  against  airaidl  structures  over  a 
range  of  inqpact  velocities  from  2.5  to  5.0  km/sec.  Initially  called  FATE  [3],  lata*  FATE-2 
[4],  and  now  FATEPEN2  [5,6],  the  code  has  bem  modified  to  include  projectile  tip  erosion 
even  at  impact  vdochies  below  shatter  vdodty.  The  equations  wnthin  the  FATEPEN2  code 
predict  the  number  of  plates  perforated  in  a  multi-plate  target  configuration  as  well  as  the 
holes  in  the  perforated  plates.  In  addition,  FATEPEN2  also  predicts  the  number,  size, 
tnyectories,  and  velocities  of  the  fragments  in  the  various  dd)ris  clouds  created  as  the  projec¬ 
tile  or  its  remains  move  through  a  multi-plate  target. 

The  PEN-4  lethality  assessment  model  was  developed  for  the  NSWC  in  an  attempt  to 
model  fiagment  impact  against  thin  plates  over  a  wider  range  of  impact  velocities  [7,8].  This 
model  is  similar  to  the  FATEPEN2  model  in  that  the  equations  used  in  the  modd  require  a 
number  of  simplifying  assumptions  and  experimentally  derived  fiictors.  By  restricting  the 
lower  limit  of  the  impact  velocity  to  3.6  km/sec,  PEN-4  is  able  to  neglect  shear  frilures  in  the 
projectile  materid;  by  restricting  its  upper  limit  to  7.6  km/sec,  PEN-4  n^ects  materid 
mdting  and  vaporization.  In  more  recent  veidons,  PEN-4  has  been  updated  to  incorporate 
advanced  fragmentation  models.  These  frt^entation  models  are  considerable  improvements 
over  the  modds  used  in  earlier  versions  of  the  code. 

KAPP-n  was  developed  for  the  Defense  Nuclear  Agency  to  predict  damage  to 
complex  three-dimensiond  aerospace  targets  impacted  by  multiple  hypervelodty  projectiles, 
induding  chunky  fragments,  rods,  and  hollow  cylinders  [10].  It  is  the  fiidon  of  the  previously 
devdoped  KAPP  and  KNAPP  computer  codes  [1 1,12]  and  has  been  calibrated  with  an 
experimentd  database  covoing  an  impact  vdodty  range  of  approx.  1-9  km/sec.  The 
empiricd  rdationships  within  KAPP-II  allow  the  user  to  characterize  the  state  of  the  projectile 
as  it  passes  through  the  target  as  well  as  the  response  of  the  target  system  to  the  impact 
loadings  of  the  inhid  projectile  and  the  dd)ris  created  by  the  initid  impact.  Unless  otherwise 
noted,  verdon  1 . 1  of  KAPP-n  was  utilized  in  the  study  whose  results  are  presented  herdn. 
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3.0  SHOCK  LOADING  AND  RELEASE  ANALYSIS 


3.1  Introductoiv  Comments 

Consider  the  normal  impact  of  a  rig^  drcular  c^indrical  projectile  on  a  flat  target 
plate.  In  this  study,  the  inojectile  is  assumed  to  be  made  of  one  or  more  perfectly  bcmded 
dissimilar  layers  or  disks.  Upon  impact,  shock  waves  are  set  up  in  the  projectile  and  target 
materials.  The  pressures  associated  vdth  these  shocks  typically  exceed  the  strengths  of  the 
projectile  and  target  nuUerials  by  several  ordoa  of  magnitude.  For  example,  in  an  8  km/sec 
ahiminum-on-aluminum  impact,  the  ratio  of  the  impact  pressure  (1 16.S  GPa=l  .15  MBar)  to 
the  streitgth  of  the  material  (3 10  MPa  for  aluminum  6061-T6)  is  approximatdy  375,  or 
roi^ify  2.5  orders  of  magnitude.  As  the  shock  waves  propagate,  the  projectile  and  targ^ 
materials  are  heated  adiabatically  and  non-isoitropically.  The  release  of  the  sh(v-?' pressures 
occurs  isoitropically  through  the  action  of  rare&ction  waves  that  are  generated  as  the  shock 
waves  interact  with  the  free  surfaces  of  du;  projectile  and  target.  This  process  leaves  the 
projectile  and  target  materials  in  high  energy  states  and  can  cause  dther  or  both  to  fragment, 
mdt  or  vqrorize,  depending  on  the  material  properties,  geometric  parameto^,  and  the  velocity 
of  impact.  At  very  early  times  during  the  impact  event,  only  the  area  in  the  immediate  vicinity 
of  the  impact  rite  is  affected  by  the  impact.  For  the  projectile  and  target  geometries  con- 
ridered  in  this  study,  the  shock  waves  can  be  conridered  to  be  initially  planar.  This 
simplification  allows  one-dimenriorud  relationslups  to  be  used  for  analyzing  the  creation  and 
release  of  shock  pressures. 

The  shock  pressures,  energies,  etc.,  in  the  projectile  and  target  materials  were 
calculated  using  the  three  1-D  shock-jump  conditions,  a  linear  relationship  between  the  shock 
wave  vdodty  and  partide  vdodty  m  each  material,  and  continuity  of  pressure  and  vdodty  at 
the  prcjectOe/target  interfru:e.  If  we  conrider  the  1-D  impact  of  a  projectile  with  velodty  Vq 
on  a  stationary  target,  conservation  of  mass,  momentum,  and  energy  across  the  shock  fronts 
in  the  projectile  and  in  the  target  yields 
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Projectile 

Iscgfil 

“q/^op  “  (“sp-WppVVp 

«st^ot*(“sr«ptyVt 

(la,b) 

Pj^  *  “^“pp^op 

PHt““st«pAot 

(2a,b) 

■■  ^*I^(Vop"Vi^))/2 

EHt-PHt(VorVHty2 

(3a,b) 

where  V-l/p  is  spedfic  volume,  Ug  and  Up  are  shodc  and  partide  vdodty,  respectivdy,  and 
P({  and  Ef]  are  the  pressure  and  energy  state  assodated  with  the  initial  inq>act.  In  equations 
(1-3),  the  subscripts  'p',  and  t*  refer  to  projectile  and  target  quantities,  respectivdy. 
Furthermore,  in  the  devdopment  of  equations  (1-3),  the  initial  conditions  ahead  of  the 
projectile  and  target  shodc  waves  were  takoi  to  be  zero  (with  the  exception  of  dendty  vdiich 
is  Po*l/Vo)  and  the  shock  vdodty  in  the  projectile  is  taken  rdative  to  a  'stationary*  projectile. 

The  linear  shock  vdodty-particle  vdodty  relationships  for  the  projectile  and  target 
materials  are  in  the  form 

Us«=Co  +  kup  (4) 

vdiere  Co^(KVo)  is  the  material  bulk  speed  of  sound,  K»E/3(l-2v)  is  the  adiabatic  bulk 
modulus,  E  and  v  are  Young's  modulus  and  Poisson's  ratio,  respectivdy,  and  k  is  an 
empirically-derived  constant.  At  the  projectile/target  interfiice,  pressure  equilibrium  implies 
that 

PHp  =  PHt  (5) 

while  materid  continuity  at  the  interfiice  implies  that 

Vo“«pp  +  «pt  (6) 

Solving  equations  (1-6)  simultaneously  yidds  exprestions  for  projectile  and  target  partide 
vdodties  vdiich  can  then  be  used  to  calculate  shodc  vdodties,  pressures,  internal  oietgies, 
and  material  dendties  after  the  passage  of  a  shock  wave. 

The  shock  loading  of  a  material  is  an  irreverdble  process  that  results  in  an  increase  of 
the  intonal  enogy  of  the  shocked  material.  However,  the  release  of  a  shocked  material 
occurs  isentropically  along  an  'isentrope'  or  'release  adiabat'.  The  difierence  between  the  area 
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under  the  isentrope  and  the  eneigy  of  the  shocked  state  is  the  amount  of  residual  energy  that 
remains  in  the  material  and  can  cause  the  nuUerial  to  melt  or  even  vaporize.  A  sketch  of  a 
generic  Hugoniot  and  a  generic  rdease  isentrope  with  initial,  shocked,  and  final  material  states 
lughlighted  is  shown  in  Figure  1.  In  order  to  calculate  the  rdease  of  the  projectile  and  target 
materiak  fix>m  their  req>ective  shocked  states  (each  characterized  by  Pf],  £h>  and  V|]),  an 
q)propriate  equation^f-state  is  needed  for  each  material.  To  keq[>  the  analysis  rdativdy 
single,  tlw  NGe-Gruneisen  [2S]  and  Tillotson  [26]  equations-of-state  were  examined  for 
suitability  fi>r  use  in  this  study. 

3.2  Mie-Gruneisen  Eouation-of-State 

The  Nfie-Grundsen  equation-of-state  ^OS)  is  an  accurate  thermodynamic  description 
of  most  metals  in  the  solid  r^ime  and  is  rdativdy  easy  to  use.  It  has  the  form 

P  =  PH  +  pr(E.EH)  (7) 

vidiere  the  time-dependent  Grundsen  coefiSdent  F  is  given  for  most  metals  as 

r-FoPo/p  (8) 

where  Fo-K3/poCp  is  the  ambioit  Grunds«i  coefiBdent,  K  is  the  adiabatic  bulk  modulus, 
P=3a  is  the  volumetric  coefiBdent  of  thermal  expandon,  and  Cp  is  specific  heat  at  constant 
pressure.  Invoking  the  Second  Law  of  Thermodynamics 

dE  =  TdS-PdV  (9) 

along  with  the  isentropic  constraint  dS^  for  the  rdease  process  allows  us  to  construct  the 
rdease  isentrope  in  P- V  space  fi>r  a  matoial  referenced  to  the  material  Hugoniot  in  P-V  space 
and  a  givoi  initial  shocked  state  defined  P{{,  Vh>  Eh-  Using  the  procedure  outlined  by 
McC^een,  et.al.  [25],  the  pressure  Pj  and  internal  energy  Ej  at  a  specific  podtion  V  along  the 
isentrope  can  be  shown  to  be  ^en  by 

Pi  -  [PK  +  (FA0i(Ei.i  -  Pi.i  AV/2  -  EHi)]/[l+(rA0iAV/2]  (10) 

vdiere  AV  is  the  incremental  change  in  volume  used  to  oreate  the  release  isoitrope,  and  Pjs 
and  E|£  are  the  pressure  and  energy  along  the  Hugoniot  corresponding  to  the  i-th  podtion  in 
the  rdease  process.  The  rdease  process  is  continued  udng  equation  (10)  until  the  release 
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iaentn^  so  detemiiiied  crosses  the  V-axis  (i  e.  until  P|  becomes  zero). 

Based  cm  hs  thermodynamic  origins,  the  Nfie-Griineisen  EOS  cannot  be  expected  to 
give  accurate  results  in  the  oqMuided  liquid  r^ime  or  in  the  vapor  regime.  This  is  because  as 
impact  energy  increases,  the  assumption  tluit  the  Gruneisen  coeffident  is  a  function  of  density 
akme  is  no  loiter  valid.  At  high  impact  ener^es,  the  Grundsoi  coeflfidoit  is  a  fiuiction  of 
internal  energy  as  well  as  dendty.  Experience  has  shown,  however,  that  it  does  yidd  fiurly  ac¬ 
curate  end-state  results  even  when  there  is  a  small  percentage  of  molten  matmal  presem  [1]. 

The  TiUotson  EOS  has  a  more  complicated  form.  In  its  original  form  [26],  it  is  has 
two  parts.  The  choice  of  which  part  to  use  dqjends  on  the  location  of  the  rdease  isentrope 
within  P-V-E  space.  The  first  part  applies  vriioi  the  material  is  in  compresrion  regardless  of 
the  internal  energy  O-O-  for  V<Vo  and  for  all  E>0)  and  in  the  small  r^on  of  expansion  in 
which  Vo<V<Vg  (Mxrvided  that  E'<^'=Es-fHv  where  Eg  is  the  total  heat  needed  to  produce 
indpient  vaporization  and  Hy  is  the  latent  heat  of  vaporization.  The  quantity  Vg=l/pg 
correrqmnds  to  the  volume  (or  density)  of  a  material  that  completes  its  rdease  process  with  an 
internal  energy  E-Eg.  In  these  two  regions,  the  TiUotson  EOS  has  the  form 

Pj  “[a  +  b/fi3E,p)]Ep  + +  (11) 

where  P“Vq/V-1  and 

<IRp)-(Pffio)(jVp)2+l  (12) 

Equation  (1 1)  appUes  in  particular  to  shock  loadmgs  in  which  the  material  remains  a 
solid  after  it  isentropicaUy  returns  to  ambient  pressure.  In  equation  (1 1),  A^q^o^  ^ 
arfb^Fo.  For  most  metals,  a  value  of  a=0.S  wUl  yidd  satisfiu;tory  results.  In  his  report, 
TiUotson  states  that  the  constants  E^.  and  B  should  adjusted  to  gve  the  best  fit  for  the  EOS 
surfiue  [26].  However,  recent  efforts  by  MidUn,  et.d.  [27]  diow  that  the  constant  B  can  be 
approximated  reasonably  wdl  as 

B-PoCo2(2k-l-ro«)  (13) 

but  that  Eq  stiU  has  to  be  treated  as  a  curve-fitting  parameter. 
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One  of  die  dangers  of  improperly  gnesang  a  value  for  Eg  is  that  the  isentrqpe  would 
actually  curve  up  frcmi  its  starting  point  CPh>^h3h)  of  curving  down  as  would  be 
oqiected.  If  this  were  to  occur,  the  rdease  process  would  have  to  be  termiiMted,  another 
value  of  Eq  would  have  to  be  spedfied  (usually  a  lower  one),  and  the  rdease  process  would 
have  to  start  over  again.  The  following  empirical  rdationship  was  obtained  as  part  of  this 
investigation  for  Eq  as  a  fonction  of  other  material  parameters  to  serve  as  a  guide  in  the 
sdection  of  an  ai^ropriate  starting  value  for  Eq. 

Eq/Ej*  -  0.819ro-0  768k6.594(T^y)-0.021(Hf/Hv)®  572  (14) 

sriiere  Hf  is  the  latent  heat  of  fusion.  This  equation  is  based  on  the  materials  considoed  in  this 
study  (see  Tables  1  and  2  for  meduuiical  and  thermal  properties,  reflectively)  and  has  a 
corrdation  coeffident  of  87.21%.  When  compared  with  the  given  values  of  Eq  used  to  derive 
h,  equation  (14)  had  an  average  error  of  2.6%  with  a  standard  deviation  of  30%. 

In  a  highly  expanded  state  0.e.  for  V>Vs  r^nnlless  of  internal  energy)  or  if  the 
internal  energy  is  lugh  enough  to  cause  complete  valorization  even  in  a  moderately  expanded 
state  Q.e.  for  Vo<V<Vg  and  if  E>Es').  the  second  part  of  the  EOS  is  invoked: 

P2  -  «Ep  +  (II*p«(E.p)  +  Ape!tp[-MV/Vo-l)l)expt<VWo-l)2)  (1 5) 

Driiere  the  constants  a  and  p  are  adjusted  to  control  the  rate  of  convergence  of  the  EOS  to 
that  of  an  ideal  gas.  The  exponential  fiictors  force  the  second  term  in  equation  (15)  to 
approach  zero  at  large  eiqiandon  volumes.  Hie  remaining  first  term  is  then  equivalent  to  the 
ideal  gas  term  (Y-l)Ep  with  y^l.S,  which  is  a  reasonable  value  for  real  gases  [26]. 

In  this  two-part  form,  the  IiUotson  EOS  is  asymptotically  correct  in  die  conqiression 
and  eiqiansion  regimes  and  reproduces  many  of  the  isentropic  rdease  features  observed  with 
mudi  more  complicated  equations-of-state  [27].  It  should  be  noted  that  the  rdease  process  as 
described  by  the  Tillotson  EOS  does  not  always  terminate  in  a  simple,  dear  cut  manner  as  it 
does  with  the  Nfie-Grundsoi  EOS.  For  impact  conditions  in  which  the  matoial  remains  in  a 
solid  state  upon  rdease,  the  isentrope  genoated  vnth  the  Tillotson  EOS  will  in  fiict  cross  the 
V-axis  in  a  manner  analogous  to  that  vriiich  is  observed  when  using  the  Mie-Grundsen  EOS. 
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However,  for  impact  conditions  that  lead  to  material  melt  and  valorization,  instead  of 
crossing  the  V>axis,  the  isentrope  created  witii  the  llllotson  EOS  l^>p^oadles  the  V-axis 
asymptotically  and  never  crosses  it.  Therefore,  an  additional  user-supplied  parameter  imist  be 
a  cut-off  point  for  the  rdease  process  in  the  event  of  extreine  gaseous  expansion. 


Table  1.  Material  Mechanical  Properties 


Material 

k 

(gm/cm^) 

WESM 

(GPa) 

V 

Ahimimim 

5.380 

1.34 

2.71 

120 

71.0 

0.35 

7.975 

1.12 

1.82 

120 

290.0 

0.08 

Cadmium 

2.307 

1.64 

864 

24 

46.2 

0.33 

Copper 

3.940 

2.49 

8.93 

37 

131.0 

0.34 

Gold 

3.060 

1.57 

19.24 

33 

85.5 

0.42 

Iron 

4.580 

1.49 

7.87 

95 

200 

0.30 

Lead 

2.030 

1.47 

11.34 

7 

13.8 

0.45 

Magnesium 

4.490 

1.24 

1.74 

45 

44.1 

0.29 

5.173 

1.22 

10.20 

200 

317.2 

0.31 

Nidcd 

4.667 

1.53 

8.86 

200 

227.5 

0.30 

Platinum 

3.680 

1.50 

21.37 

70 

191.0 

0.39 

Silver 

3.230 

2.50 

10.49 

25 

82.7 

0.37 

4340  Steel 

4.570 

1.55 

7.83 

290 

200.0 

0.30 

Tantalum 

3.374 

1.20 

16.65 

200 

179.3 

0.35 

Tm 

2.560 

1.52 

7.28 

4 

41.4 

0.33 

'ntanhim 

4.786 

1.05 

4.51 

330 

124.1 

0.30 

4.150 

1.24 

19.17 

400 

406.8 

0.30 

Zinc 

3.042 

1.50 

7.14 

82 

74.5 

0.33 
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TaUe  2.  Material  Thermal  Properties 


Material 

To 

a 

(xio-</»a 

Cp 

(cal/smOO 

Tm 

(®C) 

Tv 

(®Q 

Hf 

(cal/cm) 

Hv 

(cal/fm) 

Ahimimim 

2.13 

0.240 

0.235 

660 

2450 

95 

2450 

1.16 

0.140 

0.570 

1281 

2884 

260 

8195 

Cadnuum 

2.27 

0.343 

0.058 

321 

765 

13 

212 

Copper 

2.00 

0.170 

0.097 

1083 

2590 

49 

1150 

Gold 

3.10 

0.161 

0.034 

1063 

2960 

16 

413 

Iron 

1.57 

0.120 

0.120 

1539 

3035 

65 

1591 

Lead 

2.77 

0.293 

0.031 

327 

1740 

6 

210 

1.50 

0.300 

0.295 

650 

1110 

88 

1326 

Molybdenum 

1.52 

0.061 

0.079 

2610 

5555 

70 

1242 

Nickel 

1.80 

0.143 

0.130 

1454 

2865 

74 

1523 

natimim 

2.94 

0.110 

0.037 

1769 

4349 

26 

632 

Silver 

2.50 

0.211 

0.062 

961 

2210 

25 

554 

4340  Steel 

1.67 

0.112 

0.110 

1510 

3070 

65 

1590 

Tantalum 

1.69 

0.065 

0.033 

2996 

5425 

38 

1007 

Tm 

1.85 

0.269 

0.058 

235 

2450 

14 

580 

lltanhim 

1.10 

0.100 

0.150 

1676 

3260 

99 

2182 

Tufwaten 

1.48 

0.035 

3410 

5900 

53 

1054 

2nc 

2.15 

0.274 

0.100 

420 

907 

25 

420 

Closed-form  expressions  for  Pj  along  the  isentrope  described  by  equations  (1 1)  and 
(IS)  can  also  be  obtained  using  the  procedure  described  in  [26]  and  used  in  deriving  Pj  for  the 
Mie-Gruneisen  EOS.  Three  different  variations  of  the  incremental  form  of  equation  (9)  with 
dS^  were  ccmsidered  in  the  development  of  the  expressions  for  Pj.  These  variations  are 
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V 


« 


(#1) 

Ej  -  Ej.!  ■  -(Pj  +  I^i)AV/2 

,(16a) 

(#2) 

Ei-Ei.i-.Pi.,AV 

(16b) 

(#3) 

Ei-Ei.i-PjAV 

(16c) 

These  three  forms  were  considered  in  an  attempt  to  simplify  the  fiiud  expression  for  Pj 

In  the  procedure  described  in  [25]  equation  (1 1)  needs  to  be  manipulated  so  that  the  unknown 


inessure  ^  at  the  current  increment  is  written  in  terms  of  quantities  at  the  previous  increment, 
including  the  previous  pressure  Pj.}.  This  is  rdativdy  ea^  to  do  using  variation  (#1),  the 
nsost  sensible  of  the  three,  for  the  Nfie-Gruneisen  EOS  because  the  pressure  terms  in  the  Mie- 
Gruneisen  EOS  are  easily  separable.  In  the  HUotson  EOS,  the  complexity  arises  from  the  &ct 
dE^PdV  is  used  in  the  denominator  of  only  one  term  on  the  right-hand-ade  in  equations  (1 1) 
and  (15).  This  makes  the  separation  of  the  pressure  toms  somewdut  more  cumbersome. 

After  deriving  the  expresrions  fi>r  I^  using  each  of  the  three  proposed  variations,  the 
predictions  of  the  three  variations  for  the  impact  vdodty  required  to  produce  mdt  and 
vaporization  in  matoials  for  whidi  such  quantities  were  known  wav  compared  against  known 
vdochy  values.  It  was  found  that  variations  (#2)  and  (#3)  did  not  reproduce  the  known 
values  very  wdl.  Thus,  variation  (#1)  was  sdected  fr>r  further  use  in  the  devdopment  of  the 
equations  for  Pj.  The  final  expressions  using  equation  (16a)  are  presented  below. 

a’l)i  =  (C2-W-«CiC3)y2C,Vi  (17) 

Cl  =  Vi(AVXl+a(AV’Ari)]  (18) 

C2  -  CiRi/Vi(AV)  +  (AVA^i)Ri'  +  QiVi2(AV0  -  Pi.i(AV)Vi2[l+a(AVA^i)] 


(19) 

C3  -  (•Ei.l-K!iVi)Ri  +  bEi.iE„Vo2  -  Pi.i(AV)I(l-*»)Ei.iVi2  +  (I+b)EoV<,2  + 
QiVi3]  +  (Pi.i(AV)]2Vi2  (20) 

Qj-Ani  +  Bm^  (21) 

Ri-Ei.|Vi2  +  EoVo3  (22) 

Ri'-.Ei.iVi2  +  bEoVo3  (23) 

and  AV  -  AV/2.  Although  a  substantial  amount  of  algd)ra  is  required  to  derive  equations 
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(17-23),  the  manipulations  involved  in  deriving  a  closed-form  expression  for  (P2)i  ^an  be 
reduced  significantly  if  equation  (IS)  is  re-written  in  the  following  form: 

P2-[a  +  b'/l(E,P)]Ep  +  Q*  (24) 

ndiere  f(E,p)  is  still  given  by  equation  (12),  b'>4)U  and  Q'^US  where 


S-A|*exp[-fi(V/Vo-l)l  (25) 

U-exp[-a(VA^o-l)^]  W 

Thus,  the  expression  for  P2  can  be  written  in  exactly  the  same  form  as  the  expression  for  P]. 
As  a  result,  we  can  use  the  expressions  that  w&re  derived  for  (Pi);  can  be  used  to  give  us 
(P2)i  as  weD  provided  that  in  evoy  instance  b  is  replaced  with  bU;  and  Q;  is  replaced  with 
lljS;  where  U;  and  S;  are  fixind  using  equations  (25)  and  (26). 

3  A  yprfified  riHnt.«inn  Rqii«rion-of-State 


If  we  examine  equations  (1 1,15)  in  more  detail,  we  note  that  they  are  continuous 
across  V«Vo,  vdiich  implies  that  the  TQlotson  EOS  is  continuous  across  V=Vo  for  very  high 
impact  energies.  However,  at  V-Vs.  there  is  an  abrupt  jump  in  the  release  isentrope  for 
moderate  impact  energies,  that  is,  i^ien  £$<£<£5*  at  V=Vg.  This  jump  occurs  because  accor¬ 
ding  to  the  original  formulation  proposed  by  Tillotson,  whenever  £<£5  equation  (1 1)  is  used, 
even  in  the  Vo<V<Vs  re^on  of  the  curve.  However,  once  we  move  across  V-Vg,  equation 
(1 5)  is  invoked  regardless  of  the  impact  energy.  Since  these  two  equations  are  not  continuous 
atV»Vs,  neither  is  the  isentrope.  Table  3  shows  values  of  calculated  uting  the  Tillotson 
EOS  and  the  EOS  parameters  used  to  obtain  them.  Examitution  of  the  last  column  in  Table  3 
reveals  that  the  ratio  V^/Vo  is  relatively  insoisitive  to  the  choice  of  material:  the  average  value 
of  Vg/Vo  is  1.138  with  a  standard  deviation  of  only  4.3%  of  the  average  value. 

The  effect  of  this  discontinuity  in  the  Ullotson  EOS  is  that  it  over  predicts  the  amount 
of  expansion  that  occurs  in  the  release  of  a  material  fi’om  a  moderately  energetic  state,  that  is, 
one  that  is  not  sufficiently  enogetic  to  cause  an  appreciable  amount  of  vaporization  to  occur. 
For  example,  in  the  case  in  v^diich  E  is  only  slightly  greater  than  Eg  at  V^Vg,  the  original  form 
of  the  Ullotson  EOS  dictates  that  the  release  isentrope  for  V>Vg  would  follow  a  path  similar 


ft 
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to  the  (Nie  in  the  event  of  complete  vaporization,  that  is,  it  would  become  asymptotic  to  tlw 
V-axis  and  would  terminate  at  an  unrealistically  high  value  of  ^>ecific  volume. 


Table  3.  Values  of  V5  for  Materials  Conridered 


Material 

Initially 

MttHlDlter 

a 

P 

wsm 

V,/Vo 

Ahnaiauai 

1.1 

5.0 

5.0 

10.2 

0.369 

0.424 

1.149 

mm 

5.0 

5.0 

17.3 

0.549 

0.620 

1.129 

PiMimiiiin 

1.0 

5.0 

5.0 

3.2 

0.116 

0.128 

1.106 

Caaoer 

10.0 

10.0 

7.1 

0.112 

0.130 

1.161 

Gold 

0.3 

10.0 

10.0 

5.3 

0.052 

0.060 

1.154 

Inn 

10.0 

10.0 

7.8 

0.127 

mm 

1.141 

Lead 

0.3 

10.0 

10.0 

3.5 

0.088 

0.101 

1.148 

Maanestom 

1.0 

5.0 

7.4 

IBHII 

1.089 

0.5 

10.0 

10.0 

9.4 

0.098 

0.109 

1.112 

Nickel 

■■ 

10.0 

10.0 

8.5 

0.113 

0.133 

1.177 

Platinum 

0.2 

10.0 

10.0 

6.1 

mm 

0.053 

1.128 

Silver 

■■ 

10.0 

10.0 

4.6 

0.095 

0.122 

1.284 

Tantalum 

0.2 

10.0 

10.0 

6.0 

0.060 

0.067 

1.116 

Tin 

1.0 

10.0 

10.0 

4.9 

0.137 

wm 

1.187 

Utanium 

0.3 

10.0 

10.0 

9.0 

0.238 

1.072 

Tunasten 

0.3 

10.0 

10.0 

6.6 

0.052 

0.057 

Zinc 

1.0 

10.0 

10.0 

mm 

0.140 

0.155 

2.207 

Note;  Initial  Eq  guess  based  on  Eq  (J/kg)  =  2.56x1(Ha®-^,  A=PoCo2  [56] 


Generic  release  isentropes  obtained  by  implementing  the  Tillotson  EOS  in  its  original 
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formulation  are  illustrated  in  Figures  2a-c.  In  Figure  2a,  the  energy  as  the  isentrop^  crosses 
V*>Vo  is  less  than  Eg.  No  vi^i^jrization  is  expected  to  occur  and  calculation  of  the  isentrope 
continues  uang  equation  (1 1).  The  isentrope  in  this  case  terminates  at  a  specific  volume 
Vf<Vg.  In  Figure  2b,  the  energy  as  the  isentrope  crosses  Vq  is  greater  than  Eg  but  less  than 
Eg^  Since  E  is  already  larger  than  Eg,  the  isentrope  in  this  case  must  terminate  at  a  value  of 
spedfic  volume  greator  than  Vg  as  shown  in  Figure  2b.  In  Figure  2c,  the  energy  as  the 
isentrope  crosses  Vq  is  already  greater  than  Eg .  In  tins  case,  a  significant  amount  of 
vaporization  is  expected  to  occur.  Equation  (15)  is  invoked  automatically,  the  isentrope  is 
continuous  across  V-Vq,  and  there  is  no  jump  at  V=Vg. 

A  modification  in  the  form  of  a  'Mixed  Phase  Formulation'  of  the  Tillotson  EOS  was 
proposed  in  an  attempt  to  lessen  the  effects  of  the  discontinuity  at  V=Vg  [28,29].  The  Muted 
Phase  Formulation  states  that  if  Eg<E<Eg'  as  the  release  isentrope  crosses  V=Vo.  then  for 
V>Vo  the  i^essure  is  to  be  calculated  using  the  equation 

P3  =  [P2(E.Eg)  +  Pi(Es'-E)]/(Es'.Es)  (27) 

This  ensures  that  the  EOS  and  the  release  isentrope  are  continuous  if  E=Eg  or  if  E=Eg' 
at  V-Vg.  This  modification  was  motivated  by  the  feet  that  if  E>Eg  as  the  isentrope  crossed 
V=Vo,  then  enough  energy  would  be  present  to  cause  partial  vaporization.  Thus,  rather  than 
continue  to  use  equation  (1 1)  if  Eg<E<Eg'  when  V-Vg,  equation  (27)  is  to  be  implemented  to 
account  for  some  additional  expansion  of  the  material. 

An  alternative  means  of  eliminating  the  discontinuity  in  the  Tillotson  EOS  then  V>Vg 
and  Eg<E<Eg'  O-e.  in  moderately  high  energy  impacts)  is  uniformly  subtracting  the  magnitude 
of  the  jump  at  V=Vg  fi’om  the  pressure  values  calculated  when  V>Vg  using  equation  (15),  that 
is,  the  original  Tillotson  EOS  equation  applicable  when  V>Vs  [13].  Thus,  if  Eg<E<Eg'  as  the 
isentrope  crosses  V=Vs,  then  for  V>Vg  the  pressure  is  to  be  calculated  using  the  equation 
P4  *  P2  -  [P2(V*Vs).P3(V=Vs)]  (28) 

in  i^ch  P2  is  calculated  using  equation  (15)  and  P3  is  calculated  using  equation  (27).  As  can 
be  seen  fi’om  equation  (28),  this  correction  is  not  intended  to  replace  the  Mixed  Phase 
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Fonnulation  of  the  TUlotson  EOS,  but  rather  to  complement  its  use  in  the  region  V>Vs. 

The  quantity  within  the  square  brackets  of  equation  (28)  is  the  amount  of  the  jump  in 
the  rdease  isentrope;  it  is  largest  if  £=(£$)'*'  at  V=Vs  and  decreases  as  E->Eg .  In  the  event 
that  l^Eg'  at  V’^Vg,  the  proposed  modification  in  the  TUlotson  EOS  disappears,  the  EOS 
reverts  back  to  its  original  form  (i.e.  P4=P2)>  uid  continuity  at  is  nuuntained.  If  E<Eg 
as  the  isentrope  crosses  V^Vq,  then  the  isentrope  never  reaches  V=Vg  so  that  in  such  cases, 
the  correction  is  never  invoked.  Thus,  the  proposed  correction  is  only  invoked  when  needed, 
that  is,  if  Eg<£<Eg  as  the  isentrope  crosses  V^^Vg. 

The  effect  of  implementing  the  subtraction  jump  correction  in  the  TUlotson  EOS  on 
the  luiture  of  the  release  isentrope  is  shown  in  Figures  3a,b,  and  c  for  impact  scenarios  in 
udtich  E’=(Eg)'*',  E  is  between  Eg  and  Eg',  and  E»(Eg')“,  respectively,  as  the  isentrope  crosses 
V>»Vg.  In  Figures  3a-c,  refers  to  the  final  specific  volume  obtained  using  the  jump 
subtraction  correction  formulation  of  the  TUlotson  EOS.  As  can  be  seen  in  Figures  3a*c,  the 
proposed  modification  gives  an  appropriate  amount  of  ocpansion  when  E  is  near  Eg'  and  does 
not  over  predict  the  amount  of  expansion  when  £  is  only  slightly  greater  than  Eg. 

In  this  report,  the  version  of  the  TUlotson  EOS  in  which  the  jump  at  V=Vg  is 
eliminated  by  the  uniform  subtraction  iq)proach  is  referred  to  as  the  TUlotson/SJC  fonnulation 
whUe  the  mbced  phase  formulation  of  the  TUlotson  EOS  is  referred  to  as  the  TUlotson/MPF 
formulation.  Unless  otherwise  specified,  the  TUlotson/SJC  fonnulation  was  used  in  the 
remainder  of  this  effort.  Table  4  presents  a  summary  of  which  equation  to  use  in  which 
r^ime  of  P-V-E  ^ace  to  generate  a  release  vnth  the  TUlotson/SJC  and  TUlotson/MPF 
equations-of'State. 
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Table  4.  Equatiom  for  Generating  Rdease  Isentrope  as  a  Function  of , 
Location  in  P-V-E  Space 


nnotson/SJC 

TiOotsoii/MPF 

_ Y<Vo _ 

allE>0 

(11) 

(11) 

E<E, 

_ (11> _ 

_ 01) _ 

HI^H 

E,'<E 

(15) 

_ (15) _ 

v,<v 

_ m _ 

mBsm 

v,<v 

3.5  Propagation  of  Shock  Pressures  in  a  Multi-Material  Projectile 

As  the  shock  wave  generated  by  the  impact  on  the  target  propagates  through  the 
projectile,  it  <»icounters  the  various  interfoces  between  material  layers.  At  each  interfoce 
between  two  dissimilar  materials,  a  transmitted  shock  wave  and  a  reflected  wave  are 
generated.  The  properties  of  the  reflected  and  transmitted  waves  are  found  using  a  technique 
based  on  the  method  of  impedance  matching  (see,  e.g.,  [30-32]).  In  this  technique,  continuity 
of  pressure  and  particle  velocity  are  enforced  at  each  interface.  If  the  reflected  pressure  is 
greater  than  the  incident  pressure,  then  the  reflected  wave  is  a  shock  wave.  Conversely,  if  the 
reflected  pressure  is  less  than  the  incident  pressure,  the  reflected  wave  is  a  rarefiu;tion  wave. 
The  equations  governing  the  reflection  and  transmisaon  of  shock  waves  at  projectile  material 
interfiices  are  derived  as  follows. 

Figure  4a  shows  an  incoming  shock  wave  in  material  A,  a  reflected  wave  firom  the  A/B 
inteifoce,  and  a  trananhted  shock  wave  into  material  B.  Shock  wave  I  is  shown  mo>dng  into 
undisturbed  material  (denoted  with  a  'O'  subscript).  The  shocked  state  of  the  material  in  its 
wake  is  denoted  by  a '!'  subscript.  The  reflected  wave  n  moves  back  into  this  shocked 
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material  and  leaves  bdiind  h  material  wdiose  state  is  denoted  by  a  subscript  of  *2'.  Xbe 
tnmsmitted  shock  wave  m  moves  into  undistributed  material  (denoted  by  a  '4*  subsoipt);  the 
condition  of  the  shocked  material  behind  it  is  labded  with  a  *3*  subscript.  Figure  4b  shows  tiie 
same  configuration  only  all  motion  is  shown  under  steady  conditions. 

Across  shock  fit>nt  I  we  have 

Pl(Ui-ui)-po^l-uo)  (29) 

Pl-Po  *  Po^  l-«oXui-uo)  (30) 

as  wdl  as  the  constitutive  relationship  between  the  shock  wave  speed  Uj  and  the  particle 
vdodty  ui  induced  in  the  shodced  material 

Ui-Co^  +  kAui  (31) 

In  equations  (29,30),  po  has  been  replaced  1^  Pq\  the  am^r^  density  of  material  A. 
Assunung  stationaiy  conditions  at  zero  pressure  ahead  of  the  shock  wave  (i.e.  uo=P(f^)  and 
that  P}  is  known  (as  it  will  be  in  impact  problons),  equations  (29-3 1)  can  be  used  to  solve  for 


ui  and  PI  as  follows: 

U|  -.(CoAflkA){l.[l+4kApi/poA(c„A)2]i/2)  (32) 

Pl”Po*Ui/(Ul-ui)  (33) 

This  completely  defines  the  state  of  the  material  bdiind  shock  wave  I.  For  shock  wave  n  we 
have 

P1(U2-Hii)  *  P2(U2-Hi2)  (34) 

P2-P1  “  Pl(U2+«iX«l-U2)  (35) 

and  for  shodc  wave  m  we  have 

P3(U3-V3)  =  Po®(U3-«4)  (36) 

P3-I'4-Po®(03-»4Xu3-»4)  (37) 
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U3-Co®  +  kBu3  <38) 

where  P4  has  been  rq>laced  by  pg^,  the  ambient  denaty  of  material  B.  Assuming  stationary 
conditions  at  zero  pressure  ahead  of  shock  wave  m  and  enforcing  pressure  and  vdochy 
continuity  at  the  A/B  nuterial  interface  0-n-  ^3=^2  U3=ni2)  reduces  equations  (36-38)  to 

P3(U3-«2)*Po®U3  (39) 

P2=Po®U3U2  (40) 

U3«Co®  +  kBu2  (41) 

Equations  (34,35,39-41)  are  a  system  of  S  equations  in  6  unknowns  (P2>P2>U2>U2> 

P3,U3).  The  dimination  of  one  unknown  is  illustrated  graphically  in  Figure  S.  In  Figure  S, 
Curves  A  and  B  are  the  Ibgoniots  of  matoials  A  and  B,  and  Curve  A'  is  the  Ehigoniot  of 
fiMteriiik  A  reflected  about  point  C  which  denotes  the  initial  shocked  state  in  material  A  (i.e. 
prior  to  the  passage  of  the  reflected  wave).  The  shocked  state  of  material  B  must  lie  at  the 
intersection  of  its  Hugoniot  (Curve  B)  ami  the  reflected  Hugoniot  for  material  A  (Curve  A*). 
This  state  is  denoted  by  point  D.  The  partide  vdochy  corresponding  to  point  D  is  the 
inter&ce  vdochy  02=03  while  the  pressure  corresponding  to  point  D  is  the  inter&ce  pressure 
P2=P3. 

Knowing  that  curve  A'  is  the  reflection  of  Curve  A,  that  is,  it  passes  throughout  the 
points  (up=ui4**Pi),  (up=2ui  J*=0),  and  (up=0,P=2po^l(Co^2kAui)),  allows  us  to 
obtain  the  following  functional  form  for  Curve  A'  in  P-Up  space; 

PA'  ='  2poAuj(coA+2kAuj)  .  PoA(cj,A+4kAuj.kAup)up  (42) 

Thus,  vriien  we  set  p^'  equal  to  the  functional  form  of  Curve  B,  we  have  an  equation  for  the 
partide  velodty  that  corresponds  to  point  D.  Solving  for  this  particle  velodty  yields: 

UD  =  [Np2-4aY)l^]/2a  (43) 

iidiere 
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a  “  -  Po®k®  (^) 

3  ■  Po®Co®  4po^^l  (45) 

y-2po^l(co^2kAui)  (46) 

in  U],  of  course^  is  known.  Setting  U2*vj)  in  equations  (34,35)  and  (39-41)  allows  us 

to  solve  for  all  the  remaining  quantities: 

U3  =  Co®  +  k®U2  (47) 

P3*Po®U3U2  (48) 

P3  *  po®U3/(U3-U2)  (49) 

U2  “  (Pl-P2ypl(“l-»*2)  -  “1  (50) 

P2  “  Pl(U2‘*^iy(U2'^2)  (5 1) 

An  example  of  this  technique  is  <hscussed  in  the  next  paragraph. 


Figure  6,  ^ch  conasts  of  three  Hugoniot  curves  drawn  in  pressure-paitide  vdochy 
(P-Up)  space,  shows  vdut  happens  vdien  a  shock  wave  travding  in  copper  at  12  km/sec 
encounters  an  aluminum  inter&ce.  Using  the  three  one-dimensional  shock  jump  conditions  and 
the  linear  P-Up  relationship  for  copper,  it  b  found  that  a  12  km/sec  shock  wave  in  copper 
oeates  a  pressure  jump  of  568  GPa  and  induces  a  partide  velodty  of  5.34  km/sec  in  its  wake 
as  it  moves  into  copper  at  ambient  conditions.  To  find  the  pressures  and  partide  vdodties  of 
the  reflected  and  transmitted  waves,  the  Hugoniot  for  copper  in  P-Up  space  is  reflected  about 
the  point  defined  by  Up-5.34  km/sec,  P=568  GPa.  Its  point  of  intersection  with  the  Hugoniot 
fi>r  aluminum  yidds  the  desired  pressure  (290  GPa)  and  particle  velodty  (7. 15  km/sec)  for  the 
wave  reflected  back  into  the  copper  and  transnutted  into  the  aluminum. 

Once  the  pressure  and  the  partide  vdodty  in  a  subsequent  material  layer  are 
determined,  the  one-dimendonal  shock-jump  conditions  are  used  to  calculate  the  spedfic 
volume  and  the  eneigy  of  the  shocked  material.  This  procedure  is  repeated  for  each 
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niccessive  (Mtyectile  material  layer.  Thus,  while  the  impact  conditions  are  used  to  fl^ne  the 
shodmd  states  in  the  target  and  first  projectile  layer  materials,  the  shocked  states  in 
subsequent  projectile  material  layers  are  obtained  uring  the  impedance  matdiing  tedinique  just 
described  and  iDustrated. 

3.6  Rateias.flf  Shock  PrcMuta 

The  target  shock  pressures  are  released  by  the  action  of  the  rarefiiction  wave  that  is 
created  by  the  reflection  of  the  shock  wave  in  the  target  fi-om  the  target  rear  firee  surfiu:e. 

This  rarefaction  wave  propagates  throu^  the  target  material  and  into  the  shocked  projectile 
li^  nutoials.  In  doing  so,  it  also  releases  the  projectile  materials  fi’om  thdr  respective 
allocked  states.  For  the  purposes  of  the  modd  devdoped  herdn,  this  process  of  shocking  and 
rdeasing  continues  until  the  rarefiiction  wave  overtakes  the  shock  wave.  After  this  point  in 
time,  it  is  assumed  that  no  additiond  shoddng  and  rdease  of  projectile  material  occurs.  In  this 
manner,  the  modd  coiunders  only  material  that  is  "fiiUy  shocked". 

As  mentioned  previously,  in  some  instances  the  rdative  impedance  of  two  adjoining 
projectile  layer  matoials  may  result  in  a  shock  wave  being  reflected  back  into  a  projectile 
material  layer  that  has  been  shocked  and  rdeased.  However,  it  is  assumed  for  the  purposes  of 
tlus  study  that  this  reflected  shock  wave  does  not  "re-shock"  the  projectile  material  and  that 
the  material  into  which  h  is  reflected  renuuns  released.  This  assumption  is  reasonable  since  as 
the  reflected  shock  wave  moves  bade  into  the  released  layer  material,  it  continuously  creates 
rarefiwtion  waves  at  the  projectile  edge  fi^  surfaces  which  release  any  material  shocking  it 
produces.  Thus,  in  the  modd  developed  herein,  any  projectile  layer  material  that  has  been 
shocked  and  rdeased  will  remain  rdeased  regardless  of  the  nature  of  the  wave  reflected  fi'om 
its  interface  with  an  adjoiiung  layer. 
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SPECIFIC  VOLUME 


Figure  1.  Generic  Hugonioc  and  Release  Isentrope 
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Figure  3b.  Modified  Tlllotson  Equ«tlon-of -State  with  and  E  <E<E 


Figure  3c.  Modified  Tillotson  Equation-of-SCate  with  E-(E 


P2*P3 

U2*U3 

Rgure  4b.  Shock  Wave  Reflection  and  Transmission  at  a  Material  Interface  —  Steady  Conditions 
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4.0  DEBRIS  CLOUD  MATERIAL  CHARACTERIZATION  , 


4.1  Computing  the  Percentages  of  Solid.  Liquid,  and  Gaseous  Debris  Cloud  Matgial 

Once  the  residxial  interaal  enerpes  in  the  shocked  and  rdeased  portions  of  the  pro¬ 
jectile  and  target  materials  had  been  obtained,  the  percentages  of  the  various  states  of  matter 
in  the  resulting  dd>ris  cloud  were  estimated  uring  the  following  procedure.  This  procedure 
requires  the  knowledge  of  the  materials'  solid  and  liquid  specific  heats  (Cps>^l)>  mdfing 

and  boiling  points  (T|||,Tv),  and  their  heids  of  fiiaon  and  viq>orization  (Hf^ly)  addition  to 

file  reridual  internal  energy  (Ef). 

If  Er<CpsTin,  then  all  of  the  shocked  and  released  materials  was  considered  to  renuun 
in  a  solid  matter  state,  that  is, 

Ps=10 

Pi  =  0.0  (52a,b.c) 

Pv*0.0 

If  CpsTm<Er<Cps+Hf,  then  the  quantity  (Er-CpgTmVHf  represented  the  fiaction  of 
the  shocked  and  released  material  that  was  melted,  vdiile  the  remaining  shocked  and  released 
material  was  assumed  to  be  in  solid  form,  that  is, 

Pg  =  1.0  -  (Ei-CpgTmj/Hf 

Pi  =  (ErCpsTmyHf  (53a,b,c) 

Pv  =  0.0 

If  CpsTm+Hf<Er<  CpsTm+Hf+CplCTy-Tm),  then  all  of  the  shocked  and  released 
material  was  considered  to  be  in  a  liquid  state,  that  is, 

Ps  =  0.0 

Pi  =1.0  (54a,b,c) 

Pv*0.0 

If  CpsTm+HfK:pl(Tv-Tm)<Er<  CpsTm+HfK:pl(Tv-Tni)+Hv,  then  the  quantity 
{Ef-  [CpsTm+Hf+Cpi(Tv^Tm)]}/Hv  represented  the  fiaction  of  the  shocked  and  released 
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mterial  that  was  v^wrized,  wdiile  the  ranaining  shodced  and  rdeased  material  wa^  con¬ 
sidered  to  be  in  liquid  form,  that  is, 

Ps-0.0 

Pi  -  1.0  -  {Er[CpsTm+HfK:pl(Tv.Tm)]}/Hv  (55a,b.c) 

Pv-  {Er[CpsTm+HfK:pl(TvTm)l}/Hv 

If  CpsTni'^Hf+Cpi(Ty-T(||)+HY<^,  then  aU  of  the  shocked  and  released  material  was 
vaporized,  that  is, 

P8*0.0 

Pi  =  0.0  (56a,b,c) 

Pv»1.0 

4.2  Computing  the  Ma-wes  nf  the  Solid.  Liquid,  and  Gaseous  Dd)ris  Cloud  Material 

The  material  in  the  dd>ris  cloud  created  by  the  initial  impact  consists  of  the  target 
material  removed  by  the  impact  and  the  impacting  projectile  mass.  While  the  mass  of  the 
projectile  material  in  the  d^ris  cloud  was  known  a  priori,  the  mass  of  the  target  material  in 
the  dd>ris  cloud  had  to  be  detemuned  by  multiplying  the  target  hole-out  area  by  the  target 
thickness  and  the  target  material  denrity. 

4.2.1  Target  Plate  Hole  Diameter 

The  diameter  of  the  hole  created  in  the  target  plate  by  the  initial  impact  (D)  can  be 
calculated  using  any  one  of  a  number  of  onpirical  equations  for  hole  diameter  in  a  thin  plate 
due  to  a  high  speed  impact.  Four  such  equations  were  considered  and  implemented  in  the 
ddms  cloud  materials  diaracterization  scheme  being  developed  herein.  These  equations  are 
9vai  below. 

#nKAPP-IimSS01 

D/dp  =  exp(app){l+bvo[l-exp(-ctg/dp)]}  (57) 

vdiere  a,b,  and  c  are  empirical  constants  [10]. 


#21KAPP-II/HSS02 

D-F«(rcfl>)[T(2P-T)]>« 
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(58) 


vrfiere 


T  -  ktt,  +  Pll-tHdp/arcJ^ll/^}  (59) 

rc-P^d(pp/rt)  +  fll^  (60) 

f  S(Vo)(pp/pt)1^2(Lp/<j  .  1  +  Djn^p)  . long  rods 

P/dp-^  (61a,b) 

lS(Vo)(Pp/rt)'^irf/2ds)(Lp/dp)l/3  . disks 

S(vo)  «=  1  -  ocp{-Vo2(pj/2Yt)[u/(l-Hi)]}  (62) 

u  ®  (Pp/Pt)*^ 

Diirf/dp  -  a(pp/pt)b(3iy2dp)C[(ppVo2y(2eBt)]l/3  (64) 


where  dj  is  the  diameter  of  an  eqtnvalem  sphere,  is  the  target  material  ten^e  yield 
strength,  Bt  is  the  target  material  Brinell  Hardness  Number,  and  a,b,c,d,e,Ckt,  and  are 
enq;urical  constants  [10]. 

#3)KAPJ»-I1/HSAQ1 

D/dp  =  1  +  (Diri/dp-l){  l-expf-KW^J) 
where  Djnf  is  as  defined  in  equation  (64),  and 

h  =  a(Bt/pp)6  (66) 

and  a,b  are  onpirical  constants  [10]. 

#4)EB34.yL0 

Dftj-  n.(tt(l  .exp{-(dpft,XppVo^/Vt)®  '’'’(Pp^Pt)'®  (67) 

While  the  empirical  nature  of  these  equations  mandates  their  use  only  within  the 
impact  vdodty  r^imes  for  wUdi  they  were  designed,  the  results  obtained  fisr  vdocities  out¬ 
side  tile  presoibed  r^imes  are  in  general  not  unreasonable.  One  of  the  implications  of  these 
equations  is  that  the  amount  of  taiget  mass  in  the  debris  cloud  and  will  continue  to  grow  as 
the  impact  vdodty  is  increased.  This  is  because  the  velodty  terms  in  equations  (57,58,65,67) 
have  a  pontive  real  number  exponent;  three  of  them  state  that  hole  diameter  is  proportional  to 
¥^2/3.  Howevo*,  this  is  not  necessarily  the  case,  espedally  in  the  case  of  a  thin  target.  For  a 
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thm  target,  one  would  expect  the  hole  diameter  to  increase  until  a  certain  critical  in^pact 
velocity  (which  dq)ends  on  relative  target  and  projectile  material  and  geometric  properties) 
and  then  levd  off  as  vdodty  ouitinues  to  increase. 

Up  until  the  critical  impact  vdodty,  there  would  be  substantial  interaction  between  the 
inojectfle  and  the  target  as  the  projectile  moves  through  the  target;  above  the  oitical  impact 
speed,  the  projectile  would  move  through  the  target  so  first  (because  of  the  rdative  thinness  of 
the  target)  that  there  is  only  a  nuiumal  amount  of  projeaile/target  interaction.  Hence,  one 
would  expect  impact  velocity  to  have  a  minimal  effect  on  hole  diameter  in  a  thin  target  beyond 
a  certain  critical  value.  Unfortunately,  equations  (57,58,65,67)  do  not  have  this  charactoistic. 

A  brief  study  was  nude  using  equations  (57,58,65,67)  for  aluminum  projectiles 
impactii^  thin  aluminum  targets  at  speeds  between  2  and  25  km/sec.  The  results  are 
presented  in  Kgures  7>10;  each  Figure  corresponds  to  a  different  relative  geometric 
configuration  that  was  conadered.  In  Figure  7,  the  projectile  length-to-diameter  ratio  (Lp/dp) 
was  2  v^e  the  ratio  of  the  target  thickness  to  the  projectile  diameter  (Vdp)  was  0. 1;  in 
Rgure  8,  Lp/dp^  w^e  ts/dp^.5;  in  Figure  9,  Lp/dp=0.1  and  ts/dp=0.1;  and  in  Figure  10, 
Lp/dp»0.1  and  Vdp>K).5. 

Thus,  in  Figures  7,8  a  relatively  long  rod  impacted  a  relatively  thin  and  thick  plate, 
respectively,  while  in  Figures  9,10  a  relativdy  thin  disk  impacted  and  relatively  thin  and  thick 
plate,  respectively.  A  common  feature  of  all  four  figures  is  than  only  the  PEN4.vlO  equation 
possessed  the  ability  to  level  off  in  hole  cUameter  beyond  a  certain  impact  vdocity.  However, 
the  PEN4.V10  equation  is  for  iq>heres  only;  the  projectile  diameter  used  in  the  equation  was 
taken  to  be  equal  to  dp,  and  not  some  'equivaloit  diameter*  that  would  be  larger  than  dp  and 
confijse  the  issue.  Thus,  the  predictions  of  the  PEN4.vlO  equation  are  affected  only  by  target 
thidcness  and  not  projectile  length. 

Another  common  feature  of  all  fisur  figures  is  that  the  predictions  of  all  three  KAPP-n 
equations  continue  to  grow  as  impact  velocity  increases.  Of  these  three  equations,  the  one 
denoted  by  'KII/HSS02'  spears  to  have  some  tendency  to  flatten  out  as  the  impact  vdocity 
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increues.  Thus,  it  would  appear  that  KII/HSS02  offers  some  promise  in  being  able  to  be 
modified  to  refiect  what  would  be  expected  of  hole  diameter  as  a  function  of  inq>act  velocity. 

A  fifth  hole  diameter  option  was  added  to  the  dd>ris  cloud  characterization  scheme 
fin*  thin  disk  projectiles.  In  this  option,  the  (tiameter  of  the  hole  in  a  taiget  plate  in^mcted  by  a 
thin  hi^  speed  prcjectile  is  merdy  set  equal  to  the  diameter  of  the  impacting  disk.  This 
apinropriateness  of  this  approximation  has  been  demonstrated  in  numo-ous  expoimental 
studies  of  hi^  speed  impacts  (see,  e.g.  [33]). 

4.2.2  Calculation  of  Shocked  and  Released  Material  Masses 

To  calculate  the  masses  of  the  various  states  of  the  projectile  and  target  materials  in 
the  debris  cloud,  the  amounts  of  shocked  and  rdeased  target  and  projectile  material  had  to  be 
determined.  These  quantities  were  obtained  by  determining  the  locations  in  the  target  plate 
and  in  the  projectile  where  the  rarefiuxion  waves  had  overtaken  the  correspondii^  shock  wave 
[34].  It  was  the  material  thrc’i.^gh  which  both  the  shock  wave  and  the  rdease  wave  had 
traveled  that  was  shodced  and  rdeased  and  wdiich  was  therefore  dther  mdted  or  vaporized, 
dq)ending  on  the  particulars  of  the  impact  event.  Any  material  beyond  the  point  at  which  the 
rarefiiction  wave  had  overtaken  the  shock  vmve  was  assumed,  for  the  purposes  of  this  study, 
not  to  have  been  shodced  and  to  have  remained  in  a  solid  matter  state.  If  the  point  at  which 
the  rdease  wave  had  overtaken  the  shock  wa\«  was  beyond  the  thickness  of  the  target  plate 
or  the  length  of  the  projectile,  then  all  of  the  target  and/or  projectile  matoial  had  been 
shocked  and  rdeased. 

For  single-materid  projectiles,  refming  to  Figures  1  la,b  and  12  and  utilizing  the 
results  in  [34],  rar^iction  wave  Rj  overtakes  the  shock  wave  Sj  on  the  axis  of  symmetry  at  a 
point  in  the  projectile  given  by 

Li*0.72dp  (68) 

where  Lj  is  measured  fix>m  the  fi-ont  &ce  of  the  initially  uncompressed  projectile. 

Furthermore,  rarefiu:tion  wave  R4  will  ovotake  the  shock  wave  Si  at  a  point  in  the  projectile 
given  by 
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L4  -  ^I(cst■Hlsr^»pty(c•p•^■^^K V 

iiidMfe  t,  is  the  target  thickness,  and  c,l,c^  are  the  speeds  of  sound  in  the  shocked  target, 
prcjectile  materials  and  are  given  by  [34] 

Ci(t,p)*  -  « + l('yt,py^'p(«.p))Nt.p)l*>  W 

respectively.  Thus,  ifLi<L4,  then  Rj  overtakes  S]  first  and  the  slxKked  and  released 

|MX|jectile  length  is  taken  to  be  equal  to  Lj;  if  Li>L4,  then  R4  is  the  first  to  overtake  S]  and 
the  shocked  and  rdeased  projectile  length  is  taken  to  be  equal  to  L4. 

For  multi-material  projectiles,  the  location  in  the  projectile  where  the  rare&ction  wave 
R4  overtakes  shodc  wave  SI  is  determined  using  a  technique  derived  fii>m  that  used  fi>r 
ui^e-material  projectiles.  Conrider  Figure  13,  whidi  is  an  extension  of  Figure  12  for  a  angle 
material  iMojectile  to  the  case  of  a  multi-material  projectile.  In  Figure  13,  the  q)eeds  of  the 
waves  R4  and  S],  which  are  denoted  by  D'  and  *£*  subscripts,  respectively,  are  seen  to  dumge 
as  they  move  through  the  projectile  material  layers.  In  addition,  the  interface  velocity,  which 
is  denoted  by  a  'C  subsoipt,  is  also  seen  to  change  fiY>m  interfiu:e  to  interfiu^e  due  to  the 
dififerent  material  layer  properties. 

As  befisre,  we  are  interested  in  calculating  the  length  L4,  which  is  the  distance  form 
the  undisturbed  leading  edge  of  the  projectile  to  the  point  within  the  projectile  where  the 
rarefiiction  wave  R4  overtakes  the  shodc  wave  S]  as  it  moves  through  the  various  projectile 
liters.  This  quantity  is  obtained  by  performing  the  following  sequence  of  calculations.  It  is 
noted  that  as  we  proceed  in  the  calculations  that  follow,  the  velodties  Vb,Vq,  etc.  are 
presumed  to  be  known  for  each  successive  layer.  They  are  fiinctions  of  the  initial  impact 
ccHiditions  an  the  impedance  misnutches  at  the  various  projectile  layers.  Thus,  the  objective 
of  the  calculations  that  follow  is  to  determine  the  various  X  and  T  quantities  for  the  movement 
of  the  waves  R4  and  S]  through  a  multi-material  projectile  in  X-T  space  as  shown  in  Figure 
13. 

The  first  quantities  that  need  to  be  d^ermined  are  the  time  T(3  and  the  porition  X(3  at 
vriuch  the  rarefiiction  wave  R4  intercepts  the  target/projectile  interface.  Referring  to  Figure 
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14  rad  foUowing  the  |Ht>cedure  in  [34],  we  have 


Xc  -  ta(Vci/VA)(VB+VAy(VB+Vci) 

(71) 

Tc  -  Xc/Vci  -  (VVA)(VB+VAy(VB+Vci) 

(72) 

We  now  consider  eadi  projectile  layer  in  sequence  and  ddermine  vdiether  or  not  R4  will 

overtake  S]  whlun  a  given  l^er  or  at  some  point  b<^nd  H.  Thus,  for  the  first  layer,  referring 

again  to  Kgure  14,  we  calculate 

14(0  »  Xex  +  Xdt 

(73) 

where 

Xet  *  Xc(VeiA^C1)(Vd1‘*‘Vci)/(Vdi-Vei) 

(74) 

Tet  “  XetA^ET  •  (Xc/VciXVDl‘*‘Vciy(VDl-VEl) 

(75) 

XdT“Xet(vo/Vei) 

(76) 

Thus,  if  L4(^)  <  Lp(^),  R4  overtakes  Sj  within  the  first  projectile  layer,  otherwise,  it 


overtakes  Sj  at  some  point  beyond  the  first  layer  and  the  calculations  proceed  as  follows. 

Before  moving  on  to  the  second  layer,  we  must  first  determine  tli»  locations  of  the 
points  in  X-T  iq>ace  where  R4  and  S}  eadi  intersect  the  moving  interfiu»  between  layers  1 
and  2.  These  points  correqxind  to  points  4  and  5,  re^)ectively,  in  Figure  IS.  The  coordinates 
of  Point  5  are  detemuned  by  calculating  the  quantities  Xg]  and  T£i .  These  are  obtained 
firom  Figures  14  and  IS  using  ample  geometric  conaderations  with  the  following  results: 


XEl-VElLp('>'(VEl-*Vo)  (T7) 

Tei  =  Xei/Vei  =  Lp(')/(VEi+Vo)  (78a,b) 

To  determine  the  X-T  coordinates  of  point  4,  we  agun  refer  to  Figure  IS  and  proceed  as 
follows.  First,  noting  that 

T4*Tc  +  Tdi  (79) 

wehave 

T4  -  Tc  *Tc  +  Tdi  -  Tc  =  Tor  QU  -  XcVVdI  (80) 

so  that 

TDl-(X4-XcyVDl  (*1) 
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Second, 

T4-T5»(X4-X5y(.Vc2)  (82) 

Sid)stituting  for  T4  tccording  to  equation  (79)  and  substituting  T5=‘T£]  and  X5-X£i  allows 
equation  (82)  to  be  solved  for  X4  as  follows: 

X4-Xei-Vc2(Tc  +  Tdi-Tei)  (83) 

Equating  this  result  to  the  expresaon  for  X4  that  is  obtained  from  equation  (81)  yidds 
Xc  +  Tdi Vdi  »  Xei  -  Vc2(Tc  -  Tei)  -Vc2Tdi  (84) 

Using  equation  (83)  to  solve  for  Tj)]  yields: 

Tdi  =  [Xei  Vc2Tei  -  (Xc  +  Vc2Tc)1/(Vdi  +  ^02)  (85) 

Since  Xq  =  -TcVci  we  have 

Tdi  =  [Xei  +  vciTei  +  Tc(Vci-vc2)]/(Vdi  +  Vc2)  (86) 

Thus,  since  the  X-coordinate  of  point  4  is  given  by 

X4  =  Xc  +  TdiVdi  (87) 

the  poation  of  point  4  is  now  also  defined.  We  are  now  ready  to  the  second  projectile  layer. 
Referring  to  Figure  16  and  proceeding  as  b^ore,  we  have: 

L4(2)  =  Xei  +  Xet  +  Xpj  (88) 

XdT  “  Xdt’(Te1+TetVTeT  (89a,b) 

Xdt'^o  *  XetA^E2  (90) 

Substituting  for  X^j*  according  to  its  definition  in  equation  (89b)  yields: 

XdT  =  Xet(v</Ve2)(Te1  TEXyTET  (91) 

But  also  since  Xet^ET  ~  ^E2>  ll^s  equation  simplifies  to 

XdT  =  VoCTei  +  Tet)  (92) 

Thus,  substituting  equation  (92)  into  equation  (88)  ^elds 

L4(2) = Xei  Xet  ^oTei  +  vqTet  (93) 

Since  VqTei  =  Xjyi,  this  equation  reduces  to 

L4(2)  =  Xei  +  Xdi  +  Xet(1  +  Vj/Ve2)  (94) 

Rnally,  since  Xei  XdI  “  Lp^^\  we  have 
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L4(2)  -  Lp(l)f  Xet(1  +  v«/Ve2)  .(95) 

To  fod  X£j,  we  note  that 

TE1+TeT“Tc  +  Td1+TdT  (96) 

Since 

Tdt  *  (XCT  +  XETyVD2  (97) 

we  have 

TeT  “  Tc  +  Tdi  -  Tei  +  (XcT  +  XETyVD2  (98) 

Butsinoe 

XcT  “  Vc2(Tc  +  Tdi  -Tei)  (99) 

equation  (98)  can  be  written  as 

Tet  *  Tc  +  Tdi  -  Tei  +  [Vc2(Tc  +  Tdi  -  Tei)  +  XETyVD2  (100) 

Usiiig  the  rdadonship  Tet  Xet/Ve2>  we  (^>tain  after  siniplification  the  following 
expression  for  Xet- 

Xet  “  Ve2[(Vd2  +  Vc2y(VD2  “  Ve2))(Tc  +  Tdi  “  Tei)  (JOl) 


Substituting  this  expression  into  equation  (95)  comply  the  derivation  of  the  expression  of 
1^(2).  Thus,  if L4(2)  <  Lp(0  +  Lp(2),  R4  overtakes  Sj  within  the  second  projectile  layer,  if 
not,  we  continue  our  calculations.  As  before,  prior  to  moving  on  the  third  layer  (assuming,  of 
course,  that  it  exists),  we  must  first  determine  the  locations  of  the  points  in  X-T  space  where 
R4  and  Si  each  intersect  the  moving  inteifoce  between  layers  2  and  3.  These  points 
ccMTe^nd  to  points  9  and  10,  respectively,  in  Rgure  17.  From  Figure  16  and  17,  we  obtain 
the  fi>llowing  expressions  for  the  coordinates  for  point  10: 

Xio  *  Xei  +  Xe2  =  Xei  +  Ve2Lp^^^/(Ve2  +  ^o)  (102a,b) 

TlO  “  Tei  +  Te2  =  Tei  ^BZ^EZ  ~  TeI  Lp^^^/(VE2  '^o)  (103a-c) 

where  Xei  ^  Tei  ^  W  equations  (77)  and  (78),  respectively.  To  detomine  the 
coOTdinates  for  point  9  in  X-T  ^ace,  we  note  that 

T9  -  T4  =  (X9  -  X4)A^d2  (104a) 

and 
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T9-Tio-(X9-X,oy(-Vc3) 


(104b) 


SKnoe 

T9-T4-TD2  (105a) 

we  caii  sid>stitute  for  T4  according  to  equation  (79)  and  ob^in 

T9  «  Tc  +  Tdi  +  Td2  (105b) 

Equations  (104a,b)  are  then  combined  to  yidd: 

X9*X4  +  Td2Vd2  =  Xio-Vc3(T9-Tio)  (106a,b) 

Thus,  to  uniquely  determine  the  poation  of  poim  9,  all  that  remains  is  to  find  an  e9q)ression 
for  Tq2-  Substituting  equations  (79),  (87),  (102a),  (103a),  and  (lOS)  into  equations  (106)  and 
solving  for  Td2  yidds 

TD2“  (Xe2  +  Vc3TE2y(VD2  Vc3) 

+  [(VC2  -  Vc3y(VD2  “  Vc3)](Tc  +  Toi  -  Tgi)  (107) 

Since  X4  is  already  known,  equation  (106a)  can  be  used  to  obtain  X9.  This  completes  the 
series  of  calculafions  required  to  define  the  position  in  X-T  space  of  points  9  and  10. 

The  series  of  calculations  presemed  for  the  first  two  projectile  layers  serves  as  the 
basis  for  the  general  forms  of  the  equations  that  can  be  used  for  determining  the  location 
where  R4  overtakes  Sj  in  a  multi-material  projectile.  These  generalized  equations,  which  are 
valid  for  imperial  laym  2  through  NPMAT-1  whore  NPMAT  is  the  numbo*  of  projectile 
layers,  are  derived  as  follows. 

Reforring  to  Figure  18  and  19,  vriiidi  are  generalizations  of  Figures  16  and  17, 
re^rectivdy,  we  begin  by  writing  the  general  form  of  L4O)  as  follows: 

i-1 

-  s  xq  +*et + Xdt  (108) 

j=l 

where  the  Xgj  are  known  for  l<j<i-l  (as  are  the  accompanymg  Tgj).  Thus,  the  unknowns  in 
equation  (108)  are  Xgj  and  Xqj.  To  find  expresrions  for  these  quantities,  we  refer  to 
figure  18  and  write; 
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Tet  “  Xdt’^o  “  Xet^Vei 


,(l<»a.b) 


so  tint 

Xof-XEiCvt/VEi) 

But  also  smce 


i-1 

Xi>x’/Tet“Xdx/(S  ’^Ej  +  TEX) 
j=l 


we  have 


(110) 


(111) 


i-1 

Xdt*Xdx'(2  TEj  +  TEx)^ET  (1*2) 

j=l 

Substhuting  for  Xdx*  according  to  equation  (1 10)  and  then  rq>Iadng  Xex^ET  ^Ei 

yidds  the  following  expression  for  Xqx: 


i-1 

Xdt“Vo(]2  TEj  +  TET>TET  (*^3) 

j=l 

It  is  noted  that  equation  (1 13)  is  a  simple  generalization  of  equation  (92).  Substituting 
equation  (1 13)  into  equation  (108)  yields,  after  rq>ladng  Tex  with  Xex^Vei  and  VqTej  with 

j-1 

Xpj'  =  Xi)jTEj/2;  TeIc  (114) 

k=l 

the  following  expression  for  L4O): 

i-1  j-1 

L4<^>*Z  (XEj  +  XpjTEj/S  TEIc)  +  XEx(l+VoA^Ei)  (115) 

j=l  k=l 

Having  diminated  Xjyx  ^om  the  e^qiression  for  L4O),  all  that  remains  is  to  find  Xex-  'I'o 
b^in,  we  write 
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(116) 


i-1  i-1 

Z  T]^  +  Tet-Tc+ 2  Tpj  +  TDT 
j=l  j=l 

so  that 


i-1 

TeT  =  Tc+Z  (TDj-TEj)  +  (XcT  +  XET)A^Di  (117) 

j=I 

Notify  that 

i-1  i-1 

XcT=vciCrc+2:  Toj-  S  ('*«) 

j=i  j=i 

we  have,  after  substituting  equation  (1 17)  into  equation  (1 16)  and  simplifying,  the  following 
esqHVSsion  for  Xet- 

i-1 

XeT  ‘  VEiKVDi  +  VciVCVK  -  VEi)][Tc  +  Z  (T^  -  TEp]  (1 19) 

j=l 

Thus,  if  L^O)  <  Lp(0  +  Lp(2)  + ...  +  LpO),  th«i  R4  overtakes  Sj  within  layer  T;  if  not,  then 
we  must  determine  the  coordinates  of  the  points  in  X-T  space  where  R4  and  S 1  intercept  the 
moving  inter&ce  between  layers  i  and  i+1 .  If  we  denote  the  coordinates  of  these  points, 
^di  are  labded  11'  and  'S'  in  Figure  19,  as  (Xe4,Te4)  and  (Xsi,Tsi),  then  we  have 

i-1 

^Sl“Z^Ej  +  XEi  (120a) 

j=l 


i-1 

Tsi^Z^Ej  +  TEi  (120b) 

j=l 

idiere 
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XEi-VEjLpCO/(VB  +  Vo) 


and 

Tgi-XEiA^Ei  021b) 

Additionally, 

i-1 

Tr4  =  Tc+2;  Tpi  +  Toi  022a) 

j=l 

i-1 

Xr4*Z  XEj  +  XR4*  (122b) 

j=l 

where  the  XEj,  Tpj,  and  Tgj  are  known  for  l<j<i-l,  and  the  quantities  T];)i  and  Xe4'  are 
obtained  by  generalizing  equations  (106)  and  (107),  respectively.  Thus,  we  have 
Toi  =  (Xk  +  Vc,i+iTEi)/(VDi  +  Vc,i+l) 

i-1 

+  [(Vc,i  -  Vc.H.iy(VDi  -  Vc,i+i][Tc  +  Z  (Toj  -  Tfij)]]  (123a) 

j=l 

and 

Xr4'  =  XEi  -  Vc,i+i[Tc  +  i  CTpi  -  TEj)]  (123b) 

j=l 

At  the  last  layer,  i.e.  \^ien  i=NPMAT,  if  L^C^PbiAT)  >  then  the  entire  projectile  is 
shodced  and  released;  if  not,  then  R4  overtakes  Sj  in  the  final  projectile  material  layer. 

AO  that  remains  now  is  to  relate  the  known  quantities  Vy^,V3,V(^, ...,  etc.  to  physical 
quantities  such  as  shock  velodty,  particle  velocity,  etc.  Referring  to  [34],  these  relationships 
are  readily  obtained  and  are  presented  below. 

VA  =  «st  VB  =  Cst-Upt 

Vci  =  «pt  V£)1  =  Csp(l)  -  Upt  (124a-e) 

VEl“«sp(l)-Vo 
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Fori>2, 

Vq  ...  partide  vdodty  at  the  interftce  between  layers  i-1  and  i 
VDi-CgpCO-Vci  (125a-c) 

Yoi ...  (shock  vdocity  at  the  interfiu:e  between  laye^  i-1  and  i)  -  Vg 
It  is  noted  that  the  quantity  c^O  is  recalculated  for  each  material  layo*  based  on  the 
partide  and  shock  vdodties  obtained  for  each  layer  udng  the  impedance  mismatch  teduiique 
described  in  Section  3.5.  Finally,  if  we  substitute  the  definitions  of  Vyi^,VB,VQi,Viyi,  and 
V£i  according  to  equations  (124a-e)  into  equations  (71-75)  and  simplify,  we  detain  equation 
(69),  which  is  the  equation  obtained  in  [34]  for  single  material  projectiles. 

Rgure  20  shows  the  results  obtained  ^en  this  technique  is  applied  to  a  3-layer 
projectile  impacting  an  aluminum  plate  at  6  km/sec.  The  projectile  materials,  their  stacking 
sequence,  and  the  geometry  of  the  impact  are  also  given  in  Figure  20.  As  can  be  seen  fi'om 
Figure  20,  the  original  rarefiu:tion  wave  emanating  from  the  target  rear  surface  overtakes  the 
shock  wave  in  the  projectile  at  a  distance  of  approximately  0.71  cm  from  the  leading  edge  of 
the  undisturbed  projectile.  This  implies  that  at  the  impact  velodty  considered,  the  first  two 
projectile  layers  (i.e.  the  aluminum  and  the  sted)  are  completely  shocked  and  released  as  is  the 
first  0.202  cm  of  the  third  projectile  layer  0.e.  the  tungsten). 

It  is  the  materid  through  which  both  the  shock  wave  and  the  release  wave  travel  that 
is  shocked  and  rdeased  and  which  is  therdbre  either  melted  or  vaporized,  depending  on  the 
impact  vdodty.  Any  nuterid  beyond  the  point  at  >duch  the  rarefaction  wave  overtakes  the 
shock  wave  is  assumed,  for  the  purposes  of  this  study,  not  to  be  shocked  and  to  remain  in  a 
solid  matter  state.  If  the  point  at  which  the  release  wave  overtake  the  shock  wave  is  beyond 
the  thickness  of  the  target  plate  or  the  length  of  the  projectile,  then  all  of  the  target  and/or 
projectile  materid  is  shocked  and  rdeased.  Thus,  according  to  the  assumptions  and 
definitions  presented  herein,  the  renuuning  0.306  cm  of  the  tungsten  layer  in  the  projectile 
corresponding  to  the  impact  depicted  in  Figure  20  is  unshocked  and  unreleased. 
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In  calculating  the  amount  of  target  material  subject  to  shock  loading  and  release,  it  is 
assumed  that  the  shocked  tai;^  material  comes  from  an  area  of  the  target  equal  to  the 
presented  area  of  the  projectile  [35];  the  renuuiung  material  gected  from  the  target  in  the 
creatkm  of  the  target  |date  hde  is  assumed  to  remain  in  a  solid,  albeit  undoubtedly 
fragmented,  state.  Tins  is  due  to  the  &ct  that  if  shear  and  viscous  forces  are  n^ected,  there 
are  rx>  net  forces  actittg  on  the  projectile  and  target  masses  immediately  after  impact.  This  in 
turn  inqrlies  tiiat  the  fixoe  exerted  by  the  projectile  on  the  target  equals  the  force  exerted  by 
the  target  on  the  projectile.  Corribining  this  resuh  with  equation  (5)  and  noting  that  force  is 
the  product  of  pressure  and  area,  the  effisctive  area  of  the  target  on  which  the  inq)act  pressure 
acts  must,  to  an  first-order  approximation,  equal  the  presented  area  of  the  projectile.  This  in 
turn  inqrlies  that  the  shocked  target  nuoerial  comes  fiom  an  area  of  the  target  iq>proximately 
equal  to  the  presented  area  of  the  projectile. 

Furthermore,  it  is  also  assumed  that  the  depth  of  the  shocked  target  material  extends 
completdy  through  the  target  thickness.  Were  this  not  the  case,  then  other  target  failure 
modes,  such  as  plugging,  for  example,  might  come  into  play.  This  in  turn  would  seriously 
compromise  the  validity  of  the  assumptions  made  in  the  development  of  this  dd)ris  cloud 
model.  A  direct  consequence  of  this  assumption  is  that  the  model  developed  herein  is  not 
valid  for  "thick"  target  plates. 

Once  the  projectile  and  target  mass  contributions  to  the  dd)ris  cloud  and  the  fractions 
of  these  masses  that  were  shocked  and  released  were  obtained,  the  masses  of  the  target  and 
projectile  materials  in  each  of  the  three  states  of  matter  were  computed  by  multiplying  each 
nutter  state  percentage  by  the  appropriate  total  shocked  and  released  mass.  The  mass  of  the 
solid  shocked  and  released  nuterial  (if  aity)  was  then  added  to  the  mass  of  the  unshocked 
nuterial  0f  any)  to  obtain  the  total  mass  of  the  solid  component  of  the  dd)ris  cloud  material. 
4.2.3  Sununarv  and  Comments 

Thus,  if  we  let  Lq  denote  the  length  of  the  shocked  and  released  portion  of  the 
projectile  (origirul  length  Lp),  then  the  mass  distribution  among  the  three  matter  states  is 
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givaiby: 


lugfi 

Projectile 

Mgt-Mt-Mtsr+Mst' 

M^“Mp-Mpsr'^h4sp' 

(126a,b) 

Mgt^PstMtsr 

M^'*Pq)Mpsr 

(127a,b) 

M]t-PitMtsr 

Mip=PipMpsr 

(128a,b) 

Wvt“Pvt^tsr 

Mvp“PvpMpsr 

(129a,b) 

Mtsi^dp2tspt/4 

hdpsr^(Lo^Lp)Mp 

(130a,b) 

Mt=iiiAsp|/4 

Mp=Jtdp2Lppp/4 

(131a,b) 

indhere:  Mgt^Mgp,  and  Myt^vp  ^  masses  of  the  solid,  liquid,  and  viq>or 

components  of  the  taf:get  and  projectile  contributions  to  the  dd>ris  cloud,  req)ectively; 

Pst,Psp,  Plt>Plp>  Pvt’^vp  percottages  of  the  solid,  liquid,  and  vapor  constituents  of 

the  shocked  and  rdeased  portions  of  the  target,  and  projectile,  respectively;  and  Mpgr 
are  the  portions  of  the  target  and  projectile  that  are  shocked  and  released;  pppp  and  K^Mp 
are  the  mass  dentities  and  total  ori^nal  mass  contributions  of  the  target  and  projectile  to  the 
dd>iis  doud,  respectively;  and,  and  Mgp*  are  the  masses  of  the  diocked  and  rdeased 

portions  of  the  target  and  projectile  that  remain  in  a  solid  matter  state  upon  rdease. 

A  limitation  of  this  procedure  is  tlw  assumption  that  no  further  projectile  and/or  target 
loading  and  unloading  had  occurred  b^ond  the  point  where  the  release  waves  had  overtaken 
the  correqronding  shock  wave.  This  is  not  completely  correct  since  the  shock  wave  does  not 
simply  cease  to  exist  once  it  is  overtaken  by  a  rare&ction  wave.  Ratho-,  its  magnitude 
decreases  ova*  a  finite  amount  of  time  and  a  finite  extent  of  material.  Some  additional 
projectile  and  target  material  will  be  heated  and  possibly  mdted  until  the  strength  of  the  shock 
wave  diminishes  to  a  point  below  which  melt  due  to  plastic  deformation  no  longer  occurs. 
However,  the  procedure  set  forth  does  allow  the  calculation  of  first-order  accurate  mass 
quantities  fi>r  projectile  and  target  materials  in  the  three  states  of  matter. 
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4.3  PriffB  gpud  Vdo^cs 
4.3.1  Introductory  Comments 


The  equations  devdoped  in  the  subsequent  section  are  presented  in  their  most  general 
timn.  Th^  can  be  ^)pUed  directly  to  a  sin^e-nuiterial  projectile  and  a(Uq>ted  easily  to  iq>ply 
to  the  impact  of  a  multi-material  projectile.  In  characterizing  the  vdodties  of  the  debris  doud 
created  by  a  hypervdodty  impact  on  a  thin  plate,  there  are  two  possibilities  that  need  to  be 
considered. 

First,  all  of  the  projectile  material  is  shocked  and  rdeased.  In  this  case,  the  dd>ris 
doud  condsts  of  the  projectile  and  target  material  that  is  shocked  and  rdeased  and  the 
additiond  fiagmented  target  materid  that  is  ejected  from  the  target  plate  during  the 
perforation  process  but,  according  to  the  assumptions  made  herdn,  is  not  shocked  and 
rdeased.  In  the  debris  cloud  modd  devdoped  herdn,  all  of  this  materid  is  allowed  to  move 
axially  and  expand  radially.  The  quantities  of  inttfest  in  this  case  are  therefore  the  debris 
doud  leading  edge,  center-of-nwss,  trailing  edge,  and  expansion  velodties,  that  is,  vf,  v;,  v^ 
and  Vexp,  respective^. 

Second,  some  of  the  projectile  materid  remains,  according  to  the  assumptions 
employed  herein,  unshocked.  While  it  would  not  be  appropriate  to  call  this  unshodeed 
projectile  nuUerid  a  "residud  projectile  mass”,  it  is  reasonable  to  presume  that  this  materid  is 
less  severdy  stressed  than  that  which  is  folly  shodeed  and  then  rdeaseo.  Hence,  it  is  dso 
reasonable  to  presume  that  if  there  is  any  unshocked  projectile  materid,  then  it  does  not 
significantly  expand  radially  as  it  moves  axially.  In  this  case,  the  dd>ris  cloud  consists  of 
shodeed  and  rdeased  target  and  projectile  materids  and  the  additiond  unshocked  fi-agmented 
target  materid.  The  quantities  of  interest  are  the  dd)ris  doud  leading  edge,  center-of-mass, 
and  expansion  vdodties,  that  is,  vf,  and  Vg^p,  respectively,  and  the  velocity  of  the 

remaining  unshocked  projectile  materid,  Vp|..  Note  that  due  to  the  presence  of  the  unshocked 
projectile  mass,  there  is  no  dd)ris  cloud  trdling  edge  for  which  to  cdculate  a  velodty. 
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4.3.2 


TX*nL?lM«C*]iK: 


C-pfiyirfy  the  impact  of  a  projectile  on  a  tlun  target  and  the  dd>ris  doud  created  by  it  as 
shown  in  Figure  21.  As  indicated  in  the  Figure,  the  velocities  of  interest  are  vf,  v,,  v^xp.  and 
Vf .  As  the  initial  shodc  wave  created  by  the  impact  strikes  the  rear  surfiice  of  the  targ^  it 
creates  a  rare&ction  wave  that  travds  bade  into  the  target  and  evoitually  in  some  form  into 
the  projectile.  This  action  and  interaction  of  the  shock  wave  and  the  free  surfoce  impacts  a 
velocity  uf^  to  the  target  rear  sur&ce  equal  to  the  sum  of  the  particle  vdodty  in  the  target 
material  due  to  the  shock  wave  Upt  and  the  partide  vdodty  due  to  the  rarefaction  wave  U|^ 
foatis, 

«ftt  “  “pt  +  «rt  *  “pt  }  V(-dV/dP)isendP  (132) 

0 

where  the  P^V  curve  used  in  the  integration  is  the  isentrope  for  the  target  material.  Since 
Urt^Upt  [25],  an  alternative  form  for  equation  (132)  is 

uftt  “  2upt  (133) 

In  both  of  the  cases  described  in  Section  4.3.1,  the  vdodty  of  the  leading  edge  of  the 
dd>ris  doud  Vf  is  approximated  with  uf^  (see  also  [33]): 

Vf  =  uftt  =  «pt  J  '^(-<lV/dP)isen<lP  (1^^) 

0 

Also  common  to  both  cases  is  that  the  half-angle  measuring  the  spread  of  the  dd>ris 
doud  materials  is  ^ven  by 

e  *  tan-l(vcxpAr,)  (135) 

What  distinguidies  the  two  types  of  dd)ris  clouds  mathematically  is  the  manner  in 
which  v,,  Vexp,  and  v^  or  Vp^  are  ck;ailatcd.  When  all  of  the  projectile  material  is  shocked 
and  rdeased,  then: 

Vr  =  Vo-ufep;  (136) 
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\i  is  obtained  from  momentum  conservation  b^re  and  after  the  impact  event,  that  is, 

Vi  -  mpVo/mdc;  (137) 

and,  Vexp  is  obtained  from  the  ^jplication  of  energy  conservation  before  and  after  the  im^Mict 
event,  that  is, 

mpVo2/2-Epr  +  Etr  +  nycVi2/2  +  mdcVexp2/2  (138) 

ediere  m^«  mp-Hn^  is  the  total  dd)ri8  doud  nuss,  nip  is  the  projectile  mass,  mt  is  the  total 
target  holcKHit  nuiss,  Epr  and  Et|-  are  the  internal  projectile  and  target  energies,  respecdvdy, 
that  have  gone  into  heating  the  projectile  and  target  matmials,  and  u^«i^p-Ki|p  is  the 
vdodly  oftherearfreesurfriceoftheprcyectile.  As  in  the  case  of  is  taken  tobe 

equal  to  the  sum  of  the  partide  vdodty  in  the  projectile  material  due  to  the  passage  of  the 
shodc  wave,  Upp,  and  the  partide  vdodty  due  to  the  passage  of  the  rarefiiction  wave  in  the 
projectile  material,  Ufp,  created  by  the  reflection  of  the  shock  wave  from  the  projectile  rear 
fieesurfoce. 

In  the  event  vdien  not  all  of  the  projectile  material  is  shocked  and  rdeased,  that 
Vqxp»  and  Vp^  are  obtained  through  the  solution  of  the  following  three  timultaneous  equations: 

Vcxp'Vf-vi;  (139) 

mpVo  *  niprVpr  +  mdcVi;  (140) 

m^Vfp’12  =  Epr  +  Etr  +n*prVpr2/2  +  nycV,2/2  +  (141) 

where  in  this  case  m^c'  nip''int-mp|-  and  mpr  is  the  mass  of  the  unshocked  projectile  material. 
In  this  particular  case,  substitutiiig  for  v^q,  and  Vpf  into  equation  (141)  using  appropriate 
expressions  obtained  fit>m  equations  (139)  and  (140)  yields  a  quadratic  equation  for  v,.  This 
equation  is  then  solved  to  yidd  the  following  exprestion  for  \r,: 

v,-b/a-[(Wa)2-(c/a)]l^  (142) 

vidiere 
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t-2  +  ind</iiipr  b»vf+(inp/inpr)vo  X143a,b) 

c  -  +  (nip/mr  -  lXiV"*<ic)vo^  +  2(Epr+Etr)/nyc  (143c) 

The  quantities  and  Vpr  are  then  easily  obtained  from  equations  (139)  and  (140). 

If  the  sdution  of  the  above  system  of  equations  results  in  a  situation  vdiere  there  is  insufSdent 
energy  available  fin*  debris  cloud  expansion  or  motion  of  the  unshocked  portion  of  the 
projectile  material,  then  the  leading  edge  vdodty  is  reduced  until  some  kinetic  energy  does 
become  availaUe. 

4.3.3  Conunents 

It  is  inqx>rtant  to  note  that  equation  (136)  can  occasionally  yidd  rear  surfiu^e  vdodties 
that  nuiy  be  questionable.  For  example,  for  like-into-like  impacts,  ugp»2upp=2(vo/2>=Vo  so 
drat  equaticMi  (136)  yidds  Vf^.  However,  this  nuy  in  &ct  be  an  acceptable  result  of  one 
recalls  the  dd>ris  douds  in  the  x-ray  photographs  of  lead-on-Iead  impacts,  for  example  [24], 

In  these  photographs,  the  debris  doud  appears  to  remain  attached  to  the  target  plate,  therd}y 
giving  the  impression  that  the  rear  end  of  the  doud  does  not  move,  i.e.  that  yf=0.  In  the 
cq>per-on-aluminum  impacts  in  [33],  the  rear  end  of  the  copper  projectile  does  in  fret  move 
through  the  aluminum  target  plate  so  that  the  rear  end  of  the  dd>ris  doud  does  have  a  rather 
dear  forward  vdodty  component. 

In  addition,  equation  (136)  may  yield  negative  vdues  in  some  cases  where  a  less  dense 
projectile  impacted  a  more  dense  target  plate.  But  even  in  this  case,  pertiaps  the  negative 
vdodty  is  that  of  the  badc^lash  that  would  undoubtedly  occur  and  which  may  be  significant 
in  such  as  case.  Thus,  caution  should  be  exercised  when  using  equation  (136)  to  calculate  the 
vdodty  of  the  rear  sutfoce  of  the  debris  doud. 
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Figure  9.  Target  Hole  Diameter  Prediction  Comparison,  AluminunM>n>Aluminum  L/D^.  1,  T/D®0  I 


Hole  Diameter  Comparison,  Normal  Impact,  AL-on-AL, 
L«0.254  cm,  D»2.54  cm,  T*1.27  cm 
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Figure  10.  Target  Hole  Diameter  Prediction  Comparison,  Alumimim*on>Ahiminum  Impact.  L/DK).l,  T/D=^.5 


PROJiCTiie 


(•)  la  a  Projtetila  and  Shield  Soon  After  Impact 


PtOJiCTiLE 


(b)  After  the  Shock  In  the  Shield  Has  Reflected  From 
the  Bottom  Face  of  the  Shield 


Figure  11.  Wave  Patterns  in  a  Projectile  and  an  Impacted  Target  [34] 
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Figure  12.  X-T  Diagram  Showing  Where  a  Rarefaction  Wave  firom  the  Rear  Surftce  of  a  Target  Plate 
the  Shock  Wave  in  a  Single-Material  Projectile  [34] 
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Figure  15.  Final  First  Layer  Calculations 
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Figure  17.  Final  Second  Layer  Calculations 


Figure  19.  Hnal  i>th  Layer  Calculations 
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Figure  21 .  Debris  Cloud  Velocities 


S.0  DEBRIS  CLOUD  CHARACTERIZATION  SCHEME  VERmCATION 


5.1  Introductory  Comments 

A  FORTRAN  program  called  DEBRIS3  was  written  to  implement  the  various  pro¬ 
cedures  described  in  Sections  3  and  4.  The  source  code  is  given  in  Appendix  A,  with  input 
and  ouQxit  files  in  Appendix  B  and  C,  respectivdy.  DEBRIS3  is  an  interactive  program  that 
(nompts  the  user  fi^r  the  fisUowing  information;  1)  number  of  projectile  layers;  2)  projectile 
material;  3)  target  material;  4)  impact  velocity;  S)  target  thickness;  6)  projectile  diameter,  7) 
lei^lths  of  projectile  layers;  8)  TiUotson  EOS  option;  and,  9)  hole  diameter  option.  DEBRIS3 
also  requires  the  input  file  INDATA,  which  is  a  material  library.  INDATA  also  contains  the 
dioice  of  the  dE’^PdV  ^roxunation,  the  TUlotson  EOS  parametos  a  and  3,  and  the 
HUotson  EOS  parameter  8  vdiich  tdls  the  program  when  to  stop  a  rdease  process  in  which 
the  isentrope  is  asymptotic  to  the  V-axis.  The  units  for  the  data  in  the  file  INDATA  are 
presented  at  the  end  of  the  sample  file  in  ^pendix  B. 

DEBRIS3  genoates  the  mitput  file  IMPOUT,  ^ch  contains  a  detailed  summary  of 
the  fiallowing  information:  1)  projectile  and  target  geometric  and  material  properties;  2) 
impact  conditions;  3)  projectile  and  target  material  EOS  parameters;  4)  projectile  and  target 
matoial  end-state  calculation  results,  including  the  waste  heat  generated,  the  resulting 
temperature  inCTease,  the  percent  of  solid,  liquid,  and  vaporous  material,  and  the  masses  of  the 
solid,  liquid,  and  vaporous  components;  and,  5)  dd>ris  cloud  velocities  vf,v,,  and  v^,  and  Vgxp, 
as  applicable.  A  sample  of  the  output  file  IMPOUT  generated  by  DEBRIS3  is  given  in 
^)pendix  C.  A  word  of  caution;  while  the  TiUotson  EOS  is  relatively  straightforward  to 
inqilemmit,  its  use  requires  a  fiur  amount  of  fiunUiarity  with  its  pecuUarities. 
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<  1  SiitgWM«tenal  ProiectUes 

5.2.1  Compmaon  with  Experimental  Reauto  and  L-D  Hvdrocode  Predictions 


k 


Ddms  doud  vdodty  values  were  calculated  using  DEBRIS3  and  compared  against 
experimental  results  and  one-Klimenrional  hydrocode  predictions  obtained  from  a  prerious 
study  of  ddnis  cloud  formation  and  growth  using  thin  copper  disks  0V1>K).3)  inqMcting  thin 
aluminum  plates  [33].  As  can  be  seen  in  Table  5,  the  predictions  of  DEBRIS3  for  vf,  \r^,  and 
V|>  were  in  excellent  agreement  with  those  of  the  1*D  hydrocode  and  the  experimental  results. 


Table  S.  Comparison  of  DEBRIS3  with  Empirical  Results  and  1-D  Hydrocode  Predictions 
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(1)  E3q)aimental  Results  [33]  (2)  1-D  Hydrocode  Predictions  [33]  (3)  DEBRIS3  Predictions 
Over  all  the  cases  conridered,  the  average  difference  betwem  the  predictions  of 
DEBRIS3  for  vf,  and  Vf  and  the  corre^nding  ^cperimental  results  was  approximately  4% 
with  a  standard  deviation  of  q)proximatriy  3%.  However,  the  predictions  of  DEBRIS3  for 
^exp  exceeded  the  experimental  results  by  an  average  of  iq>proximately  40%  with  a  standard 
deviation  of  approxinutdy  1S%.  This  discrepancy  may  have  been  due  to  the  foct  that  the 
expansion  vdodty  measured  in  [33]  was  that  of  the  heario'  copper  component  of  the  dd>ris 
doud  while  the  expansion  vdodty  calculated  by  DEBRIS3  was  based  on  both  ddrris  cloud 
materials. 
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Figure  22  presents  a  comparison  of  the  predictions  of  DEBRIS3,  the  hydrocode  CTH, 
and  the  semi-empirical  code  FATEPEN2  for  dd)ris  cloud  leading  edge  velocity  vf  for  steel 
c^inders  (L/I>*1)  normally  impacting  thin  aluminum  target  plates  (T/IX).  12S).  As  is  evident 
in  Figure  22,  the  predictions  of  DEBRIS3  compare  &vorably  with  those  of  FATEPEN2  in  the 
vdodty  regime  for  v^ch  FATEPEN2  was  designed  to  be  used  (i.e.  less  than  approximatdy  S 
to  6  km/sec).  A  quick  calculation  reveals  that  the  difference  between  the  DEBRIS3 
predictions  of  leading  edge  vdochy  and  those  of  FATEPEN2  was  iq>proximately  26%  of  the 
DEBRIS3  values  with  a  standard  deviation  of  approximately  4%  for  the  impact  vdocides 
considered.  One  reason  for  this  difference  could  be  the  fact  that  the  mass  of  the  impacting 
proiecdle  conadered  (approximatdy  1S5S  grains  =  100  gms)  exceeded  the  maximum  value  of 
prqiectQe  masses  used  to  devdop  the  FATEPEN2  equations. 

The  CTH  values  plotted  in  Figure  22  are  average  values  of  the  vdodties  of  three 
Lagrangian  station  points  along  the  impact  centerline  vnthin  the  aluminum  target  plate.  These 
average  values  differed  from  the  connesponding  minimum  and  maximum  values  by 
iq^roximatdy  0.5  km/sec  at  an  impact  speed  of  2  km/sec  and  3.0  km/sec  at  an  impact  speed 
of  14  km/sec.  Inspection  of  Figure  22  also  reveals  that  there  is  excellent  agreement  between 
the  predictions  of  DEBRIS3  and  CTH  for  dd)ris  cloud  leading  edge  velocity.  The  average 
difforence  between  the  DEBRIS3  and  CTH  values  was  approximately  4%  of  the  DEBRIS3 
values  with  a  standard  deviation  of  approximately  3%. 

Figure  23  presents  a  comparison  of  the  predictions  of  DEBRIS3,  CTH,  FATEPEN2, 
PEN4,  and  KAPPII  for  dd)ris  cloud  half-angle  for  steel  cylinders  (L/D=1.0)  normally 
inq)acting  thin  aluminum  target  plates  (T/D=0. 125).  In  Figure  23,  the  average  difference 
between  the  predictions  of  KAPP-II  and  DEBRIS3  was  approximatdy  18%  of  the  DEBRIS3 
value  with  a  standard  deviation  of  approximately  10%;  the  average  difference  between  PEN4 
and  DEBRIS3  was  {qjproximately  6%  with  a  standard  deviation  of  nearly  7%.  Based  on 
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time  ranilts,  it  nuy  be  argued  that  the  ivedictions  of  DEBRIS3  agree  fiurly  well  with  those  of 
KAPP>n  and  PEN4.  However,  comparing  the  dififerences  between  DEBRIS3  and 
FATEPEN2  was  somev^uU  more  difiBcult  because  FATEPEN2  distinguishes  between  target 
debris  qnead  and  proje^e  ddms  spread  wdiile  DEBR1S3  does  not.  In  FATEPEN2,  the 
target  dd)ris  half-ap^e  is  fixed  at  25^  v^e  the  fxoje^e  debris  half-apgle  is  based  on 
material  properties,  impact  conditions,  etc. 

It  is  interestmg  to  note  that  unlike  the  smooth  curve  predictions  of  KAPP-II, 
FAIEPEN2,  and  PEN4,  the  curve  representing  the  growth  of  the  d^ris  cloud  spread 
generated  by  DEBRIS3  contains  numerous  kinks.  In  particular,  the  impact  velocities 
corresponding  to  the  vertical  lines  in  Figure  23  also  correspond  to  impact  velocities  at  which 
significant  changes  occur  in  the  way  the  initial  kinetic  energy  of  the  projectile  is  distributed  to 
various  competing  meduuiical  and  thermal  processes  during  the  impact  event.  These  features 
of  the  curve  predicted  by  DEBRIS3  are  discussed  in  the  following  paragraph. 

For  the  impact  considered  in  Figure  23,  between  2  and  S  km/sec,  increasing  the  impact 
vdodty  resulted  in  a  steady  increase  in  dd)iis  cloud  spread.  However,  at  5  km/sec,  the  target 
nnaterial  b^an  to  melt.  As  a  result,  some  of  the  additional  kinetic  energy  of  the  initial  impact 
provided  as  impact  velodty  increased  beyond  5  km/sec  was  used  up  by  the  target  material 
state  change  and  was  not  available  for  dd>ris  cloud  expanaon.  Thus,  the  rate  of  dd>ris  cloud 
expanaon  slowed,  and  the  slope  of  the  curve  decreased  as  impact  velocity  increased  beyond  S 
km/sec.  Between  8  and  9  km/sec,  the  projectile  material  began  to  melt  and  the  target  material 
b^an  to  viqiorize.  This  fiuther  decreased  the  rate  of  dd>ris  cloud  expansion.  However,  once 
the  projectile  was  completdy  melted,  the  rate  of  debris  cloud  expansion  increased.  Near  12 
km/sec,  the  projectile  tiuderial  b^an  to  experience  vaporization.  The  rate  of  d^ris  cloud 
expansion  slowed  down  only  slightly  because  by  now  the  dd)ris  cloud  consisted  of  a 
significant  amount  of  hot  vrqwrous  material.  By  IS  km/sec,  the  debris  cloud  was  nearly  all 
vq)or  caunng  its  rate  of  expansion  to  increase  dramatically. 
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S.3  Muhi-Miteriil  Projectiles 
S.3.1  Introductory  Comments 


The  vafidhy  of  the  mutti-material  modding  capability  of  DEBRIS3  was  assessed  by 
conqiaring  the  predictitMis  of  DEBRIS3  against  experimental  and  numerical  data.  The 
experimental  data  and  the  results  of  one  series  of  hydrocode  simulations  were  obtained  from  a 
study  that  analyzed  the  ^fectiveness  of  layered  projectiles  against  re-entry  vdiicle-type 
targets.  The  results  of  a  second  series  of  hydrocode  runs  were  obtained  using  the  CTH 
hydrocode  spedfically  for  the  present  investigation.  The  results  of  this  validation  exercise  are 
presented  in  the  next  three  sections. 

5.3.2  Comparison  with  Experintental  Results 

Three  high  qreed  inq>act  tests  were  performed  at  4  km/sec  using  three  different  equal- 
weight  projectiles  [36].  The  first  was  a  solid  7.S  gm  TAIOW  (i.e.  a  tantalum  alloy  with  10% 
tungsten)  sphere,  while  the  second  and  third  projeoiles  were  7.5  gm  layered  spheres  with  a 
solid  TAIOW  core  surrounded  by  a  steel  shell.  The  outer  shell  of  the  second  projectile  was 
1018  steel  (i.e.  nuld  strength  steel)  while  that  of  the  third  projectile  was  4340  sted  (i.e.  a  high 
strength  steel). 

In  simulating  these  three  impact  tests  with  DEBRIS3,  the  layered  spheres  were 
modeled  as  cylindrical  projectiles  with  three  layers.  The  middle  layer  corresponded  to  the 
^herical  core  i^e  the  first  and  third  layers  represented  the  outer  shell  material.  The 
thicknesses  of  the  first  and  third  layer  were  set  equal  to  the  outer  shell  thickness.  The 
thidcness  of  the  inner  layer  and  the  diameter  of  the  cylindrical  projectile  were  calculated  by 
settii^  the  inner  layor  thickness  equal  to  the  t^lindrical  projectile  diameter  and  then  solving 
for  the  diameter  by  equating  the  mass  of  the  cylindrical  projectile  to  the  mass  of  the  original 
layered  sphop. 
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In  addition  to  ad^ting  the  geometry  of  the  ori^nal  projectiles  used  in  the  test  ^es  to 
a  projectile  geometry  that  was  compatible  with  DEBRIS3,  some  compromises  were  also  made 
r^arding  the  projectile  and  target  materials.  In  the  original  test  series,  the  target  was  a  2-D 
flat  plate  representation  of  a  half-scale  re-entry  vducle,  i.e.  a  layer  of  silica  phenolic  bonded  to 
a  thin  layer  of  aluminum.  Since  the  current  version  of  DEBRIS3  does  not  allow  for  multi¬ 
material  targets,  the  targets  used  in  the  DEBRIS3  impact  simulations  did  not  have  the  outer 
layer  of  rilica  phenolic.  In  addition,  ^lile  a  witness  block  was  placed  behind  the  initial  multi¬ 
layer  target  plate  in  the  experimental  tests  to  record  the  damage  of  tht  perforating  projectile 
and  target  dd>ris,  DEBRIS3  was  not  developed  to  have  a  predictive  crq>ability  for  damage  to 
subsequent  witness  blocks  or  plates.  Finally,  whereas  one  of  the  original  projectile  materials 
was  a  tantalum  alloy  with  10%  tungsten,  the  corresponding  material  in  the  DEBRIS3  impact 
rimularions  was  pure  tantalum. 

As  expected,  the  simplifications  described  in  the  previous  two  paragraphs  precluded 
any  direct  comparison  of  the  predictions  of  DEBRIS3  and  the  experimental  results.  However, 
it  was  possible  to  make  qualitative  comparisons  of  the  DEBRIS3  predictions  and  the  actual 
test  results  because  the  simplifications  maintained  some  similarity  between  the  original  test 
materials  and  configurations  and  the  materials  and  geometries  of  the  DEBRIS3  impact 
rimulations.  These  qualitative  comparisons  became  possible  after  the  DEBRIS3  predictions 
were  analyzed  to  infer  the  relative  severity  of  the  damage  levels  that  could  have  been  expected 
on  subsequoit  witness  plates  had  they  been  placed  behind  the  initial  target  plate. 

First,  DEBRIS3  predicted  that  a  significant  portion  he  target  material  would  be 
melted  when  impacted  by  the  solid  tantalum  sphere.  Alternatively,  when  impacted  by  the 
layered  projectiles,  DEBRIS3  predicted  that  the  target  material  would  be  shocked  and 
released  but  would  return  to  a  solid  state  of  matter.  This  indicates  that  the  target  material 
would  probably  be  firagmented  but  not  melted.  Second,  DEBRIS3  pr^icted  that  the  kinetic 
eneigy  of  the  remaining  unshocked  projectile  material  would  be  greatest  for  the  layered 
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projectile  with  a  tantalum  core  and  a  high-strength  steel  shell  and  would  be  least  for  the  solid 
tantalum  projectile.  Taken  together,  these  two  features  indicate  that  the  cylindrical  projectiles 
simulating  the  layered  sphere  projectiles  would  probably  inflict  more  severe  damage  on  a 
witness  plate  behind  the  target  than  would  the  cylindrical  projectile  simulating  the  solid 
TAIOW  projectile.  This  agrees  with  the  actual  test  results,  i^ch  state  that  among  the  three 
impact  tests,  the  crater  dqjth  and  volume  in  the  witness  block  behind  the  target  impacted  by 
the  projectile  with  a  TAIOW  core  and  a  4340  sted  sheU  were  greatest  and  those  in  the  block 
bdiind  the  target  impacted  by  the  solid  TAIOW  projectile  were  least. 

5.3.3  Comparison  vwth  Hvdrocode  Predictions  —  First  Series 

The  first  series  of  numerical  runs  consisted  of  two  sets  of  three  high  speed  impacts  at 
1 1  km/sec  using  the  SOIL  hydrocode  [36].  The  projectiles  used  were  similar  in  construction 
to  those  in  the  previously  discussed  experimental  tests  Q  e.  one  solid  and  two  layered  spheres 
in  each  test  set).  The  m^or  distinguishing  feature  between  the  two  sets  of  impact  simulations 
in  this  series  is  the  mass  of  the  projectiles:  45  gm  projectiles  were  considered  in  the  first  set, 
while  5  gm  projectiles  were  used  in  the  second  set.  In  both  sets  of  simulations,  the  solid 
sphere  was  made  out  of  tungsten  as  was  the  core  in  the  layered  spheres;  the  shells  of  the 
layered  spheres  were  made  out  of  different  strength  steels.  In  modding  the  SOIL  impact 
Emulations  vEth  DEBRIS3,  simplifications  in  the  projectile  and  target  geometries  were  made 
similar  to  those  in  the  previous  section.  As  a  result,  the  following  comparisons  are  again  only 
qualitative  in  nature. 

As  in  the  DEBRIS3  simulations  of  the  experimental  tests,  the  DEBRIS3  simulations  of 
the  SOIL  runs  indicated  that  the  solid  projectiles  would  melt  some  of  the  target  plate  material 
wdiereas  the  layered  projectiles  would  not.  In  addition,  the  kinetic  energies  of  the  unshocked 
projectile  nmterials  firom  the  layered  projectiles  greatly  exceeded  those  of  the  unshocked 
projectile  materials  fi'om  the  solid  projectiles.  These  two  features  again  indicate  that  the 
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l^«red  |HX>jectiles  would  inflict  more  severe  damage  on  a  witness  plate  behind  the  targiet  than 
would  an  equal-weight  solid  projectile. 

Interestingly  enough,  while  the  general  trends  observed  in  the  DEBRIS3  inqract 
simulations  agreed  with  the  hypothesis  that  motivated  the  layered  projectile  investigation,  they 

* 

disi^ree  vinth  the  actual  numerical  results  obtained  as  part  of  that  investigation.  The 
corresponding  SOIL  runs  predicted  that  the  witness  block  damage  due  to  the  impacts  of  the 
sdid  projectiles  would  be  approximatdy  the  same  as  the  damage  caused  by  the  layered 
projectiles.  Apparently,  either  the  impact  and/or  geometric  parameto^  used  in  the  SOIL  runs 
masked  subtle  diffoences  in  damage  levels  resulting  from  the  solid  and  layered  projectile 
imqMCts  and  prevented  them  from  being  discernible,  or  the  DEBRIS3  modeling  of  the 
projectile  and  target  geometries  over-emphasized  some  impact  phenomenology  that  produced 
some  differences  in  response  that  would  otherwise  have  been  negligible. 

In  any  event,  it  is  q)parent  that  additional  testing  of  multi-material  projectile  that  are 

*  compatible  with  the  modeling  crqiabilhies  of  DEBRIS3  are  required  to  iiilly  validate  the 
predictive  capabilities  of  DEBRIS3.  As  an  intermediate  step,  several  CTH  runs  were 
performed  using  projectile  and  target  geometries  that  were  ideally  suited  for  and  matched  to 
the  aq)abilities  of  DEBRIS3.  The  results  of  these  runs  and  how  they  compared  with  the 
predictions  of  DEBRIS3  are  discussed  in  the  next  section. 

5.3.4  Comparison  with  Hvdrocode  Predictions  -  Second  Series 

In  the  second  series  of  hydrocode  luns,  four  high  speed  impact  simulations  were 
performed  at  10  km/sec  using  CTH  with  multi-material  cylindrical  projectiles.  The  projectile 

*  diameter  and  target  plate  thickness  were  kept  constant  at  2.54  cm  and  0.3 175  cm, 
respectively.  In  the  first  two  runs,  the  layers  w^e  relatively  "thin"  (i.e.  L/D=0. 1  each),  while 

*  in  the  second  two  runs,  the  projectile  layers  wore  relatively  "thick"  (i.e.  L/D=l .0  each).  In  the 
first  and  third  runs,  an  aluminum  target  plate  was  impacted  by  a  projectile  with  an  aluminum 
leading  layer,  a  4340  steel  middle  layer,  and  a  tungsten  rear  layer.  In  the  second  and  fourth 
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funs,  the  order  of  the  projectile  mateials  was  reversed.  A  detailed  descript'on  of  the  iippact 
and  geonwtric  parameters  are  given  in  Table  6;  the  results  of  the  DEBRIS3  impact  simulations 
and  the  corre^nding  CTH  results  are  given  in  Table  7.  In  Tables  6  and  7,  a  '1'  in  the  first 
cdumn  rrfers  to  the  leading  layer  of  the  projectile  while  a  *3'  refers  to  the  rear-most  projectile 
layer. 


Table  6.  Geometric  and  Impact  Parameters  for  DEBRIS3  and  CTH  Comparison  Runs 


Run  No. 

1 

(1) 

(2) 

(3) 

(4) 

10.0 

10.0 

10.0 

10.0 

D(cm) 

2.54 

2  54 

2.54 

2.54 

T(cm) 

3.175 

3.175 

3.175 

3.175 

Target  Material 

Aluminum 

Aluminum 

Aluminum 

Aluminum 

Layer  1  Material 

Aluminum 

Tungsten 

Aluminum 

Tungsten 

Layer  2  Material 

4340  Steel 

4340  Steel 

4340  Steel 

4340  Steel 

Layer  3  Material 

Tungsten 

Aluminum 

Tungsten 

Aluminum 

Lj  (cm) 

0.254 

0.254 

2.54 

2.54 

L2(cm) 

0.254 

0.254 

2.54 

2.54 

L3  (cm) 

0.254 

0.254 

2.54 

2.54 

L/D 

0.3 

0.3 

3.0 

3.0 

Proj.  Mass  (gms) 

38.24 

38.24 

382  40 

382.40 

The  predictions  of  CTH  and  DEBRIS3  regarding  the  state  of  the  target  and  projectile 


layer  materials  were  compared  quantitatively  and  qualitatively.  To  facilitate  quantitative 
comparisons  of  material  end-states,  average  densities  were  computed  for  each  material  layer 
using  the  DEBRIS3  and  CTH  results.  The  DEBRIS3  values  were  obtained  by  multipl^g  the 
mass  of  shocked  and  released  nuterial  by  its  final  dendty,  adding  to  it  the  product  of  the 
den^  of  the  unshocked  material  and  its  mass,  and  then  dmding  by  the  total  mass  of  the 
material  layer  under  consideration.  The  CTH  v^ues  are  simply  average  values  through  the 
particular  layer  thickness  and  were  obtained  fi^om  density  history  plots  along  the  centerline. 
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A  feature  common  to  all  four  impact  simulations  and  evident  in  Table  7  is  that  tfie 
average  taiget  material  densities  predicted  by  DEBRIS3  were  significantly  higher  than  those 
preActed  by  CTH.  However,  the  reason  for  this  is  that  they  include  the  solid  component  of 
target  material  not  considered  to  be  shocked  and  released  by  the  impact  (i.e.  the  remainder  of 
the  qected  target  material  not  swept  out  by  the  projectile).  The  contributions  of  the  solid 
material  component  to  the  average  density  of  the  taiget  material  are  agnificant  considoing 
that  th^  constitute  approximatdy  90%  of  the  target  material  in  the  dd>ris  cloud  created  by 
the  impact.  If  the  target  hole  diameter  had  been  set  equal  to  the  projectile  diameter  (which  is 
not  an  unreasonable  assumption  for  the  impact  vdocity  and  geometries  considered),  thoi  there 
would  not  have  been  any  unshocked  target  nuaoial  and  it  is  reasonable  to  presume  that  the 
average  densities  of  the  taiget  material  would  have  been  much  doser  to  the  CTH  values. 


Table  7.  Comparison  of  DEBRIS3  and  CTK  Impact  Response  Predictions 


Run  No. 

(1) 

m 

(3) 

(4) 

DEBRIS3 

CTH 

DEBRIS3 

era 

DEBR1S3 

era 

DEBR1S3 

era 

Vf 

(km/sec) 

10.63 

14.27 

17.36 

14.20 

10.65 

11.08 

16.26 

14.28 

9(d^ 

36 

32 

34 

27 

48 

37 

21 

22 

Plarg- 

(gm/cm-*) 

2.56 

■i 

2.36 

-0.0 

2.63 

1.33 

2.52 

-0.0 

(gm/cm-*) 

2.02 

-0.0 

18.52 

0.67 

2.22 

2.25 

18.98 

9.20 

P2 

(gm/cm^) 

umiii 

m 

7.25 

5.50 

7.83 

8.53 

7.83 

6.94 

P3  - 
(gm/cm-’) 

18.42 

m 

2.48 

1.41 

19.17 

17.19 

2.71 

3.39 

The  differences  between  the  DEBRISl.  predictions  of  dd)ris  cloud  leading  edge 


vdodty  and  the  corresponding  CTH  values  in  lUins  No.  1-4  are  25.4%,  22.3%,  3.9%,  and 
13.9%,  respectivdy,  of  the  CTH  values.  The  somewhat  large  differences  in  Runs  No.  1  &  2 
may  be  explained  by  the  following  considerations.  In  the  characterization  scheme  employed 


77 


fay  DEBRIS3,  the  target  shock  loading  and  release  analysis  used  to  obtain  the  debris  cloud 
leading  edge  velocity  is  truly  one-dimoisional.  That  is,  it  is  performed  using  only  the  leading 
prq^ctile  layer  and  the  target  nuUerial;  anything  bdiind  the  first  projectile  material  layer  is 
ignored.  In  the  case  ofthick  projectile  layers,  the  use  of  one-dimensional  equations  is 
q^ropriate  because  the  rear  laym  of  the  projectile  are  sufficiently  far  firom  the  impact  ^  so 
as  not  to  affect  the  magnitude  of  the  veiodty  of  the  target  rear  firee  surfiu:e.  However,  in  the 
case  of  the  thin  projectile  layers,  the  second  and  third  projectile  layers  are  close  enough  to  the 
projectile-target  interfiu:e  to  influence  the  shock  and  release  process  in  the  target  material  and 
the  resultii^  velocity  of  the  target  rear  fi^ee  surfiice.  CTH,  bdng  a  3-D  hydrocode,  is 
apparently  sensitive  to  these  effects  while  DEBRIS3,  being  a  first  principles  code,  is  not.  As  a 
result,  the  CTH  and  the  DEBRIS3  predictions  differ  somewhat  more  in  Runs  No.  1  and  2  and 
are  more  in  agreement  in  Runs  No.  3  and  4. 

The  differences  between  the  DEBRIS3  predictions  of  debris  cloud  leading  edge 
vdodty  and  the  corresponding  CTH  values  in  Runs  No.  1-4  are  1 1%,  26%,  23%,  and  5%, 
req>ectivdy,  of  the  DEBRIS3  values.  The  CTH  predictions  of  debris  cloud  half-angle  were 
obtained  indirectly  firom  d^ris  cloud  output  plots.  In  some  cases,  the  precise  angles  were 
(UfBcult  to  d^ermine  fix>m  the  CTH  plots  because  not  all  of  the  dd>ris  cloud  material  was 
retained  by  CTH  and  subsequently  plotted.  If  tha-e  is  a  very  small  fiaction  of  a  material  in  a 
cdl  in  which  more  than  one  material  is  present,  then  it  is  possible  for  that  small  fiaction  of 
material  to  goierate  n^ative  internal  iAieigies  in  that  cell.  CTH  allows  the  user  to  set  a  flag 
that  forces  CTH  to  drop  the  cell  fiom  subsequent  calculations  in  such  cases.  If  this  is  not 
done,  then  in  such  cases  the  time-stq}  becomes  so  small  that  the  impact  simulation  will  be 
fi)rced  to  termiiude  prematurely,  ^parently,  in  Runs  No.  2  and  3,  CTH  dropped  a  fair 
amount  of  cdls  as  the  calculations  proceeded  which  in  turn  produced  rather  sparse  d^ris 
clouds.  Whfle  the  agreement  between  the  DEBRIS3  predictions  and  the  CTH  values  was  in 
general  fiurly  reasonable,  this  may  explain  in  part  why  in  Runs  No.  2  and  3  the  DEBR1S3 
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predictions  of  the  ddms  cloud  half-angle  vahMS  were  significantly  higher  than  the  Cn| 
vahies. 

For  projectiles  with  thin  layers,  ahhou^  Table  7  indicates  that  there  were  significant 
dififerences  between  the  material  densities  predicted  by  DEBRIS3  and  those  obtained  with 
CTH,  closer  examination  of  the  DEBRIS3  and  CTH  predictions  of  material  state  did  in  fiu:t 
reveal  a  qualitative  agreement  in  the  results.  For  example,  the  extremely  low  material 
densities  for  Run  No.  1  predicted  by  CTH  indicate  that  the  material  fi-om  the  three  projectile 
kr/ers  in  both  cases  are  in  highly  expanded  states.  However,  the  denaty  of  the  rear-most 
portion  of  the  third  material  layer  in  Run  No.  1  (approximately  the  last  33%)  was  nearly  3.S 
times  that  of  the  fi^rward  portion  of  that  layer,  indicating  that  the  rear  third  of  the  final 
incjectile  layer  vras  significantly  more  dense  than  the  rest  of  the  projectile  material. 
Interestingly  enou^  for  Run  No.  1,  DEBRIS3  predicted  that  the  first  two  nuderial  layers 
would  be  in  a  liquid  state,  while  the  last  25%  of  the  third  layer  would  not  be  fiilly  diocked. 
Thus,  while  the  actual  den^  values  may  have  been  different  (vdiich  was  not  totally 
unexpected  given  the  relatively  simple  nature  of  the  phyacs  employed  by  DEBRIS3),  there 
was  some  agreement  between  CTH  and  DEBR1S3  with  regard  to  the  state  of  the  projectile 
material  foOowing  the  initial  impact. 

With  r^ard  to  the  target  material,  CTH  predicted  that  the  target  material  would  be  in 
a  highly  expanded  state  in  Run  No.  1  and  probd)ly  vaporized  in  Run  No.  2;  DEBRIS3 
predicted  that  the  target  material  would  be  compl^y  melted  in  Run  No.  1  wdiile  in  Run  No. 
2  it  would  be  partially  vaporized  as  wdl.  Thus,  there  was  again  some  general  agreement 
between  CTH  and  DEBRIS3  r^arding  the  state  of  the  target  material  following  a 
hypervdodty  impact  of  a  projectile  whh  thin  nutorial  layers. 

For  projectiles  with  thick  material  layers,  projectile  nuterial  characterizations 
predicted  by  DEBRIS3  were  again  found  to  agree  in  a  general  sense  with  the  post-impact 
material  states  predicted  by  CTH  (Table  7).  For  example,  DEBRIS3  predicted  that  in  both 
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Run  No.  3  and  Run  No.  4  that  the  second  and  third  projectile  layers  would  remain  unsl^ked 
«4ule  part  of  the  first  layer  would  be  shodced  and  released.  The  CTH  results  for  Runs  No.  3 
ft  4  clearly  showed  the  third  material  layers  to  be  relatively  undisturbed  and  the  second 
material  layers  to  be  only  slightly  deformed.  These  characteristics  are  also  evident  in  Table  3 
vdiere  the  density  values  predicted  by  CTH  for  the  second  and  third  projectile  material  layers 
were  near  ambient  values;  the  densities  predicted  by  DEBRIS3  for  the  second  and  third  laym 
were  naturally  exactly  equal  to  the  respective  ambient  values  due  to  the  assumptions  within 
the  DEBRIS3  model. 

With  regard  to  the  first  material  layer,  DEBRIS3  predicted  that  in  Run  No.  3  the  entire 
shocked  and  rdeased  portion  would  be  all  liquid,  whereas  the  shocked  and  released  portion  in 
Run  No.  4  would  be  a  mixture  of  liquid  and  solid  material.  The  CTH  results  for  both  cases 
showed  that  the  density  of  the  leading  edge  of  the  first  layer  was  approximately  30%  of  the 
ambient  value,  while  the  density  of  the  rear  portion  of  the  leading  layer  approached  the 
ambient  value  of  the  second  layer  nuterial,  indicating  a  significantly  more  compressed  state 
than  that  of  the  leading  edge. 

Some  mteresting  features  are  also  evident  in  the  CTH  and  DEBRIS3  predictions  of  the 
state  of  the  target  nuiterial.  In  Run  No.  3,  the  average  target  material  density  as  predicted  by 
CTH  is  qjproximately  12%  of  ambient.  This  indicates  a  significant  liquid,  if  not  vaporous, 
component  of  the  target  material  in  the  d^ris  cloud.  For  Run  No.  3,  DEBRIS3  predicted  that 
100%  of  the  shocked  and  rdeased  target  material  would  be  liquid  and  that  the  density  of  the 
shocked  and  rdeased  target  nuterial  would  be  approximately  75%  of  ambient.  The  near-zero 
value  of  the  target  nuterial  density  as  predicted  by  CTH  in  Run  No.  4  indicates  a  material 
state  near  complete  v^rarization  for  the  ejected  target  material  while  DEBRIS3  predicted  that 
q}proximatdy  24%  of  the  shocked  and  released  target  nuterial  would  be  in  a  vapor  state  and 
that  76%  would  be  liquid.  The  density  of  the  shocked  and  rdeased  target  material  predicted 
by  DEBRIS3  was  40%  of  the  ambient  value  in<hcating  a  highly  expanded  nuterial  state. 
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Figure  22.  Debris  Cloud  Leading  Edge  Velocity  Comparisons 


Debris  Cloud  Half-Angle  Comparisons 
Steel  Projectile  (L/D«  1.0),  Aluminum  Mte  (T/D*0.125) 


« 
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Figure  23.  Debris  Cloud  Half-An^e  Comparisons 


6.0  ADDITIONAL  RESULTS  AND  DISCUSSION 


6.1  End  State  Calculations 

Rg^res  24-27  compare  the  results  of  the  release  process  for  alumimiin-on-aluininum 
tnq>acts  at  three  different  energy  ievds  using  Mie-Gnindsen,  Tillotson,  Tillotson/SJC,  and 
HUotson/MPF  equations-of-state.  Figure  28  shows  the  differences  in  final  qiedfic  volume 
obtained  using  the  hfie-Gruneisen,  Tillotson,  and  Tillotson/SJC  equations-of-state. 

In  Rgure  24,  the  rdease  process  as  desoibed  by  the  Mie-Grundsen  EOS  and  the 
tillotson  EOS  are  nearly  identical.  This  is  to  be  expected  for  relatively  low  eneigy  impact 
(i.e.  those  impacts  in  \diich  the  materials  return  to  a  solid  matter  state  after  rdease).  Figure 
2S  shows  the  dramafic  difference  between  udng  the  Mie-Grundsen  EOS  and  the  Tillotson 
EOS  fisr  very  high  energy  impacts  (i.e.  those  impacts  in  which  the  materials  vaporize).  The 
hfie-Grundsen  EOS  caimot  account  for  the  expandon  of  the  gaseous  state  and  terminates  the 
rdease  process  at  a  much  lower  specific  volume  than  the  Tillotson  EOS. 

Figure  26  highlights  one  of  the  difficulties  in  using  the  Tillotson  EOS  in  its  original 
fi^rmulation.  This  difficulty  occurs  under  impact  conditions  that  are  not  violent  enough  to 
viqiorize  the  material,  yet  are  strong  enough  to  cause  the  material  to  melt  and  be  in  an  energy 
state  that  is  near  indpient  vaporization.  Under  these  conditions,  the  jump  in  the  release 
isentrope  at  V==Vs  generated  by  the  original  Tillotson  EOS  and  the  implementation  of  the 
Mixed  Phase  Formulation  both  result  in  a  final  volume  that  is  artifidally  high.  As  stated  previ¬ 
ously,  the  final  volume  was  considered  to  be  artifidally  high  because  the  jump  at  forced 

the  rdease  isentrope  to  follow  a  path  as  if  complete  vaporization  of  the  material  had  occurred. 
Some  viq>orization  will  indeed  occur  if  the  internal  energy  at  V=Vs  is  greater  than  that 
required  to  initiate  viqiorization  of  the  material.  However,  there  is  no  need  for  the  release 
isentrope  to  follow  the  path  of  complete  vaporization  unless  the  internal  energy  is  greater  than 
that  required  for  complete  vaporization. 

Implementation  of  the  jump  correction  given  by  equation  (28)  in  this  impact  energy 
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r^iine  caused  the  rdease  processes  to  terminate  at  specific  volume  values  that  w«e  much 
more  reasonable.  It  is  noted  that  this  correction  had  no  effect  when  the  impact  energy  was 
rdativdy  low  or  veiy  high.  Figure  27  shows  the  resuh  of  implementing  the  jump  correction 
given  by  equation  (28)  for  a  10  km/sec  alumimim>on-aluminum  impact.  In  a  such  a  scmario, 
a  fiur  amount  of  mdting  and  expansion  would  be  expected  to  occur.  The  Tillotson  EOS 
release  isentrope  shown  in  Figure  27  after  implementing  the  correction  is  more  reasonable 
because  it  terminates  at  a  spedfic  volume  that  is  greater  than  that  predicted  by  the  Mie- 
Gfundsen  EOS  which  cannot  accoum  for  greatly  expanded  states,  yet  is  substantially  less  than 
that  w^Jch  would  be  obtained  following  the  path  of  complete  vaporization 

Figure  28  illustrates  the  differences  in  the  final  q)ecific  volumes  obtained  in  aluminum- 
on-ahiminum  impacts  uang  the  TiUotson,  Tillotson/SJC,  and  Tillotson/MPF  EOS 
formulations  and  contrasts  these  with  the  results  obtained  using  the  Mie-Gruneisen  EOS.  For 
impact  vdodties  bdow  approximately  9  km/sec,  the  results  were,  as  expected,  nearly 
identicd.  For  impact  velocities  above  approximatdy  24  km/sec,  the  final  values  predicted  by 
the  Tillotson  EOS  and  the  two  dtemative  formulations  of  the  Tillotson  EOS  overiiq>  and 
significantly  exceeded  those  predicted  by  the  Mie-Grundsen  EOS  due  to  the  gaseous  expan¬ 
sion  of  the  released  materid  at  those  impact  vdodties. 

The  odd  behavior  in  the  find  vdues  of  specific  volume  due  to  the  jump  in  the 
unmodified  Tillotson  EOS  began  for  duminum-on-duminum  impacts  at  an  impact  velocity  of 
approximatdy  9  km/sec.  However,  the  Tillotson/SJC  formulation  produced  a  smooth 
transition  as  the  nuUerid  changes  fi’om  a  solid  state  (below  approximately  6  km/sec)  to  a 
liquid  state  (between  approximately  6  and  1 1  km/sec)  to  a  gaseous  state  (above  approximately 
1 1  km/sec).  The  spedfic  volumes  cdculated  by  the  Tillotson/MPF  formulation  closely 
fijUowed  those  of  the  hfie-Grundsen  EOS  until  an  impact  velodty  of  approximately  18 
km/sec  beyond  which  they  began  to  diverge  rapidly.  Apparently,  the  Tillotson/SJC 
formulation  predicted  a  more  ^anded  materid  end-state  than  did  the  Tillotson/MPF 
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fiomnilatioii  for  impact  vdodties  between  10  and  24  km/sec. 

6  2  Dehria  rioud  Materi^  TharacteriTtinn 


Figures  29-32  compare  the  effects  of  using  the  hfie-Gnineisen,  Tillotson/SJC,  and  the 
TUlotson/MPF  formulations,  respectively,  to  calculate  the  percentages  of  the  three  nuUter 
states  in  dd>ris  clouds  created  by  aluminum-on-aluminum  inq>acts.  While  the  results 
presented  and  discussed  apply  only  to  aluminum-on-aluminum  impacts,  the  equations 
devdoped  herdn  may  be  used  to  estimate  the  state  of  the  material  within  a  dd>ris  cloud 
created  by  the  high  speed  impact  of  virtually  any  two  materials  for  which  the  required  matoial 
properties  are  available. 

As  can  be  seen  in  Figure  29,  the  Mie-Grundsen  equation-of-state  predicted  only  a 
small  amourn  of  v^)orized  materid  at  an  impact  vdodty  of  25  km/sec.  However,  Figures  30 
thrcHigh  32  reveal  that  the  original  formulation  and  both  modified  versions  of  the  Tillotson 
equation-of-state  predicted  that  aluminum  was  completely  vaporized  at  an  impact  vdocity 
between  20  and  25  km/sec.  This  difference  is  due  to  the  &ct  that  the  Mie-Grundsen 
equation-of-state  did  not  account  for  the  expandon  of  the  material  it  neared  vaporization  and 
completed  the  rdease  process  with  the  material  in  a  much  lower  energy  state  than  did  either  of 
the  two  modified  versions  of  the  Tillotson  equation-of-state. 

Comparing  Figures  30-32  reveals  that  the  Tillotson,  Tillotsmi/SJC,  and  the 
Tillotson/MPF  formulations  agreed  in  the  percentages  of  the  various  states  of  matter  at 
qteeds  bdow  ^prox.  1 1  km/sec  and  above  approximately  24  km/sec.  Within  the  1 1-24 
km/sec  impact  vdodty  regime,  the  original  formulation  of  the  Tillotson  EOS  predicted  a 
steady  growth  in  the  amount  of  vaporized  nuiterial.  Within  the  Lame  impact  velodty  regime, 
the  lillotson/MPF  formulation  predicted  viqmrization  to  develop  more  rapidly  than  did  the 
Ullotson/SJC  formulation  which  predicted  a  more  gradual  transition  to  vaporized  material. 
This  appears  to  contradict  the  results  shown  in  Figure  28  in  which  the  Tillotson/MPF 
formulation  initially  predicted  a  more  dense  debris  cloud  than  did  the  Tillotson/SJC 
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foniHjlatton.  In  fact,  the  question  is  raised  as  to  how  a  more  expanded  material  sta^e  can  have 
less  vapor  than  a  less  expanded  state,  especially  ance  both  Tillotson  EOS  formulations 
predicted  approximately  the  same  radial  expansion  of  the  dd>ris  cloud  nuUerial! 

The  resolution  of  this  tq>parent  dilemma  lies  in  the  values  of  the  leading  and  trailing 
edge  vdodties  as  predicted  by  the  two  alternative  versions  of  the  Tillotson  EOS.  In  the  12  to 
24  km/sec  impact  vdodty  r^ime,  the  Tillotson/SJC  formulation  predicted  values  of  v^  and  vf 
that  were  smaller  and  larger,  respectivdy,  than  the  corresponding  values  predicted  by  the 
Tlllotson/MPF  formulation.  Hence,  in  this  impact  velocity  regime,  the  Tillotson/SJC 
formulation  predicted  the  dd)ris  cloud  to  have  a  larger  axial  dimension  than  did  the 
Tlllotson/MPF  EOS  while  the  radial  dimension  in  both  cases  was  approximately  the  same. 

This  naturally  resulted  in  larger  ddiris  doud  volumes  with  the  Tillotson/SJC  formulation  than 
with  the  Tlllotson/MPF  formulation,  even  though  the  vapor  content  predicted  by  the 
TUotson/SJC  formulation  was  less  than  that  predicted  by  the  Tillotson/MPF  formulation. 

6.3  Distribution  of  Solid.  Liquid,  and  Gaseous  Projectile  and  Target  Materials 

In  Figure  33,  the  total  projectile  mass  remained  constant  because  the  projectile  length 
and  diameter  were  fixed  in  all  of  the  impact  scenarios  conddered.  The  solid  dark  region 
represents  the  mass  of  the  projectile  that  was  undiocked  and  therefore  was  not  subjected  to 
melting  and/or  vaporization.  This  quantity  increased  with  impact  velocity  because  the  speed 
of  the  rarefaction  wave  in  the  projectile  increased  at  a  fiister  rate  than  did  the  speed  of  the 
shock  wave  in  the  projectile.  As  the  impact  velodty  increased,  the  rarefaction  wave  caught  up 
with  the  shock  wave  within  a  shorter  period  of  time.  This  in  turn  increased  the  amount  of  the 
projectile  material  that  was  not  subject  to  melting  and/or  vaporization.  The  remaining  shaded 
areas  in  Figure  33  show  the  amounts  of  the  diocked  and  released  projectile  material  in  each  of 
the  three  matter  states  as  the  impact  velocity  increased  fi'om  4  to  25  km/sec. 

Figure  34  shows  that  the  amount  of  target  material  in  the  debris  cloud  increased  as 
impact  velocity  increased  due  to  the  growth  in  target  hole  size  as  impact  velocity  increased. 
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For  the  projectile  and  target  geometries  con^dered,  all  of  the  target  matoial  was  shocked  and 
released.  Hence,  th«%  is  no  solid  daik  area,  only  the  three  lighter-shaded  areas  ^ch 
the  amounts  of  shocked  and  released  material  in  each  of  the  three  states  of  matter. 

6.4  DAris  Cloud  Vdodties 

Figures  35  and  36  compare  the  difio'ences  in  calculated  debris  cloud  velocities  for  thin 
didc  and  long  rod  impacts,  respectively.  In  the  event  of  the  disk  impact,  the  quantities  plotted 
leading  edge,  trailing  edge,  center-of-mass,  and  expansion  velocities  relative  to  the  initial 
impact  velocity.  In  the  case  of  the  rod  impact,  the  quantities  plotted  were  leading  edge, 
center-of-mass,  expansion,  and  unshocked  projectile  mass  velocities  also  relative  to  the  initial 
impact  velocity. 

In  the  thin  disk  impact  (Figure  35),  the  length  of  the  disk  was  equal  to  the  thickness  of 
the  target  plate.  In  this  case,  all  four  normalized  velocity  components  remained  relatively 
constant,  with  minor  increases  and  decreases,  respectively,  in  the  normalized  leading  edge  and 
trailing  edge  velodties,  respectively.  This  implies  that  the  changes  in  the  various  components 
of  the  dd)ris  cloud  velocity  field  in  the  event  of  a  straight-on  thin  disk  impact  are  directly 
proportional  to  changes  in  the  initial  impact  velocity.  Taken  together,  the  slight  increase  in 
the  leading  edge  velocity  and  the  slight  decrease  in  the  trailing  edge  velocity  (both  relative  to 
the  impact  velocity)  indicate  that  the  dongation  of  the  dd>ris  cloud  becomes  more  and  more 
pronounced  as  the  impact  velodty  is  increased. 

In  the  long  rod  impact,  (Figure  36),  the  length  of  the  rod  was  equal  to  four  times  the 
thickness  of  the  target  plate.  In  this  case,  all  four  normalized  velocity  components  changed 
dramatically  as  the  impact  velodty  was  increased.  As  discussed  in  Section  5.2.2  for  the  same 
taiget/projectile  combination,  many  of  the  changes  evident  in  Figure  34  coindded  with  impact 
vdodties  at  which  changes  occurred  in  the  way  in  which  the  kinetic  energy  of  the  impacting 
projectile  was  distributed  among  the  various  competing  mechanical  and  thermal  processes 
during  the  impact  event.  Closer  examination  of  Figure  36  reveals  some  additional  features  of 
interest.  These  are  discussed  in  the  following  paragraph. 
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First,  tlwre  is  a  significant  decrease  in  the  normalized  velocity  of  the  unsho<^ed 
projectile  material  while  the  normalized  leading  edge  velocity  increases  significantly  as  the 
impact  vdodty  is  increases  beyond  S  km/sec.  Taken  together,  this  implies  that  the  leading 
edge  of  the  dd>ris  cloud,  vsdiich  contains  the  more  molten  tuui  vaporous  material,  becomes 
increasingly  separated  fi'om  the  trailing  unshocked  projectile  material  as  the  impact  vdodty  is 
increased  into  the  hydrodyiuunic  r^ime.  The  accompanying  rise  in  the  normalized  expanaon 
vdodty  indicates  that  this  leading  non-solid  material  is  also  being  spread  out  to  a  greater  and 
greater  extent  as  the  impact  vdodty  increases.  Second,  the  normalized  cent^-of-mass 
vdodty  remained  relatively  constant,  with  a  value  i^>proximately  equal  to  that  in  Figure  35 
O.e.  in  the  case  of  a  thin  disk  projectile).  This  was  expected  dnce  the  cento’-of-nuss  velocity 
is  based  on  a  conservation  of  momentum  calculation  which  is  relatively  insensitive  to  changes 
in  energy  distribution  among  the  competing  processes  du  'ng  the  impact  event. 
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Figure  24.  Low  Energy  Impact  Shock  Loading  and  Rdea%  Curves  (H ...  Hugoniot,  MG ...  hfie*Gruneisen  EOS 
Release  Isentrope,  T ...  Tillotson  EOS  Release  Isentrope 
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Figure  25.  High  Energy  Impact  Shock  Loading  and  Rdease  Curves  (H ...  Hugoniot,  MG ...  Mie-Gruneisen  EOS 
Release  Isentrope,  T  ...  Tillotson  EOS  Release  Isentrope) 
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Figure  26.  Moderate  Energy  Impact  Shock  Loading  and  Release  Curves  (H ...  Hugoniot,  MG ...  Mie-G 
EOS  Release  Isentrope,  T ...  Tillotson  EOS  Release  Isentrope,  MT ...  Modified  TiUotson  EOS  Release 
Isentrope) 
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Figure  28.  Final  Specific  Volume  vs.  Impact  Velocity  (MG ...  ^fie-Gruneisen  EOS  Values,  T ...  Tillotson  EOS 
Values,  MT ...  Modified  Tillotson  EOS  Values) 
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Debris  Cloud  Composition,  AL-on-AL  Impact 
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Figure  30.  Debris  Cloud  Material  Composition  Using  the  Tillotson  EOS,  Aluminum-^n- Aluminum  Impact 
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Figure  3 1 .  Debris  Cloud  Material  Composition  Using  the  Tillotson/SJC  EOS,  Aluminum-on-Aluminum  Impact 
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Figure  32.  Debris  Cloud  Material  Composition  Using  the  Tillotson/MPF  EOS,  Aluminum-on-Aluminum  Impact 
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Figurj  33.  Projectile  Material  Mass  Distribution,  Aluminum-on-Aluminum  Impact 


Target  Material  Mass  Distribution,  Aluniinunw>n- Aluminum  Impact 
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Figure  35.  Dd>ris  Cloud  Velocities,  Ttun  Disk  Impact,  Steel-on-Aluminum 
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Figure  36.  Debris  Cloud  Velocities,  Long  Rod  Impact,  Steel-on-Aluminum 


7.0  SUMMARY  AND  RECOMMENDATIONS 

7.1  Summary 

A  robust  lethality  assessment  nwthodology  nuist  indude  the  effects  of  discrete  particle 
in^Mcts  as  wdl  as  the  re^nse  of  the  target  to  in^Mildve  debris  doud  loadings.  A  first-order 
accurate  sdieme  has  been  implemented  to  determine  the  amount  of  material  in  each  of  the 
three  states  of  matter  in  a  dd)ris  doud  created  by  a  hypervelodty  impact  on  a  thin  target.  A 
modified  version  of  the  Ullotson  EOS  was  used  to  calculate  the  residual  taergy  in  the  projec¬ 
tile  and  target  materials  upon  rdease  fi’om  their  respective  shocked  states.  Elemoitary 
thermorfynamic  prindples  were  use  to  deternune  the  percentages  of  shocked  and  rdeased 
projectile  and  target  materials  that  were  mdted  and/or  vaporized  during  die  rdease  process. 
Using  assumed  projectile  and  target  geometries,  these  percentages  were  then  used  to  calculate 
the  mass  of  the  projectile  and  target  matoials  in  solid,  liquid,  and  gaseous  form.  Dd>ris  cloud 
velocities  were  calculated  using  the  prindples  of  momentum  and  energy  conservation;  the 
spread  of  the  debris  doud  nuterial  was  thoi  readily  obtained. 

The  predictions  of  the  ddiris  cloud  modd  were  compared  against  experimratal  data, 
the  predictions  of  three  different  empirically-based  codes,  and  against  the  predictions  of  1-D 
and  3-D  hydrocodes.  In  general,  the  predictions  of  the  characterization  scheme  devdoped 
herdn  compared  fitvorably  with  the  experimoital  results,  the  lethality  assessment  schemes' 
predictions,  and  the  predictions  of  the  hydrocodes.  While  some  of  the  details  in  the  dd>ris 
doud  modd  dififered  firom  empirical  evidence,  H  is  noted  that  the  dd>ris  doud  model 
presented  herein  was  devdoped  solely  through  the  application  of  fimdamental  physical 
prindples  without  any  empirical  'adjustmem'  factors.  In  this  light,  the  agreement  between  the 
dementary  theory  predictions  and  the  experimental  results  is  highly  encouraging. 

7.2  Recommendations 

Based  on  the  work  completed  thus  far,  the  following  reconunendations  are  offered  for 
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ocm^Miing  the  devdopment  of  a  lethality  assessment  modd  that  would  be  i^licab\e  in  in^Mct 
scenarios  vdiere  materid  mdt  and/or  vaporization  can  be  expected  to  occur. 

The  next  stqp  in  the  first-order  charactmzation  of  the  ddnis  douds  oeated  in  a 
hypervdodty  impact  would  be  to  determine  the  nature  of  the  dd>ris  cloud  solid  fi:agment 
population.  This  indudes  cdculating  the  numbor  of  projectile  and  target  materid  fiag^  ^ts, 
as  well  as  their  rizes,  speeds,  and  tnyectories.  In  addition  to  the  fi'agmentation  model. 
FATEPEN,  PEN-4,  and  KAPP-II,  tlw  fi'agmentation  modds  developed  by  Grady,  et  d.  [37- 
39]  can  be  used  to  predict  the  number  of  fi'agments  that  would  result  firom  a  KEW  impact. 

The  predictions  of  the  various  fi'agmentation  modds  can  be  compared  against  one  anotlwr  and 
against  available  experimentd  data  to  detomine  vriiidi  fi'agmentation  modd  is  best  suited  for 
use  in  a  lethality  assessment  methodology.  Hypervdodty  impact  test  results  for  a  variety  of 
target  tystems  are  avdlable  fi’om  a  number  of  sources,  including  NASA  [40],  NSWC  [41], 
NRL  [42],  BRL  [43],  and  others  [44-50]. 

After  a  sads&ctory  first-order  accurate  procedure  that  characterizes  debris  doud 
composition  is  completed,  the  accuracy  of  the  procedure  needs  to  be  improved.  This  indudes 
modifying  the  methods  presented  herdn  to  include  a  more  appropriate  hole  diameto-  predictor 
equation,  the  impact  of  non-monolithic  projectiles  that  are  more  representative  of  actud  KEW 
geometries,  and  the  impact  of  yawed  and/or  obliqudy  inddent  projectiles.  Additional 
modifications  to  improve  the  accuracy  of  the  dd>ris  cloud  calculations  are  as  follows. 

First,  the  method  of  calculating  the  percentages  of  projectile  and  target  matmal  in  the 
three  states  of  matto*  should  also  be  rqrlaced  vrith  a  more  rigorous  thermodynamic  procedure. 
One  method  (see,  e.g.  [24])  would  require  calculating  the  entropy  of  the  shocked  state,  that  is, 
the  entropy  imparted  to  the  material  by  shocking  it  to  a  pven  pressure.  The  material  will 
retain  that  entropy  during  isoitropic  release  to  the  final  release  pressure  and  specific  volume. 
The  calculation  is  completed  by  identifying  the  material  state  with  that  entropy  at  the  final 
rdease  pressure  by  consuhipg  classical  thomodynamic  tables  (see,  e.g.  [51,52]). 

Second,  a  diode  wave  attenuation  procedure  [34,53,54]  should  be  implemented  to 
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cbUm  more  accurete  mass  values  for  the  material  that  is  mdted  and/or  vaporized  ip  a  high 
qteedKEW  impact  Such  a  procedure  will  result  in  a  residual  energy  profile  akmg  the  lei^ 
of  the  prcyectile  and  tlvougfi  the  thidoiess  of  the  target.  Energy  levds  at  various  portions 
can  then  be  compared  to  energy  levds  necessary  to  b^jn  material  meh  or  vaporization.  In 
addition,  the  assumption  that  the  inqMCt  pressure  acts  on  an  area  equal  to  the  area  of  the  hole 
created  in  the  target  plate  needs  to  be  reconridered. 

Third,  in  its  present  formulation,  it  is  oitirdy  possible  that  the  value  of  the  parameter 
Eq  in  foe  TUlotson  EOS  can  be  different  for  different  impact  velocities  even  when  tlw 
progectOe  and  target  nuterials  are  hdd  constant.  Since  Eq  is  part  of  an  EOS  and  an  EOS  is  a 
material  property,  the  value  of  Eg  should  be  constant  and  should  not  depend  on  impact 
conditions.  If  Eg  were  to  diange  vnfo  a  change  in  impact  conditions,  this  would  imply  the 
existence  of  an  EOS  airfime  that  also  changes  with  impact  conditions,  which  is  not  posrible 
[55].  Thus,  it  is  imperafive  to  address  the  manner  in  >duch  foe  value  of  Eg  is  chosen  in  foe 
appficafion  of  the  Tillotson  EOS. 

Subsequent  to  the  devdopment  of  a  satisfoctory  ddrris  cloud  characterization  scheme, 
an  impulsive  loading  algorithm  for  foe  target  should  be  developed  to  account  for  the  effects  of 
the  non-solid  ddrris  doud  constituents  as  wdl  as  foe  solid  non-perforating  dd>ris  doud  fig¬ 
ments.  This  effort  requires  as  input  foe  masses  and  velodties  of  foe  non-solid  dd>ris  cloud 
materials,  the  area  of  foe  inner  wall  over  which  the  impulsive  loading  is  applied,  and  tl^ 
geometric  and  nnateriai  properties  of  the  inno’  wall,  including  foe  spacing  between  the  outer 
and  inner  walls  and  the  orientation  of  the  inner  wall  with  respect  to  that  of  the  outer  wall. 
Issues  to  be  addressed  include  vriiefoer  foe  impacts  of  the  target  and  projectile  debris  doud 
materials  need  to  be  considered  separately  or  can  be  considered  simultaneously,  whether  foe 
effects  of  the  molten  and  vaporous  debris  doud  components  need  to  be  considered  separatdy 
or  can  be  combined,  and  how  to  account  for  foe  decrearing  time  of  foe  load  qrplication  and 
the  increaang  area  over  which  it  is  applied  as  the  initial  impact  velodty  increases. 

The  impulsive  loading  algorithm  can  be  validated  at  velocities  attainable  using  existing 
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hypervdocity  launchers  by  comparing  the  predictions  of  the  algorithm  with  availal^e  impact 
test  data.  The  algorithm  can  be  modified  if  necessary  until  a  satisfiictoiy  levd  of  accuracy  is 
reached.  It  can  then  be  combined  with  the  dd}ris  doud  charactoization  sdieme  and  a  suitable 
finagmentation  modd  to  yidd  an  improved,  robust  lethality  usessment  method  for  high  speed 
KEW  impacts. 
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SOIBDO 

PBOORMI  DBBRI83 

THIS  PROGRAM  PERFORMS  THB  FOLLOWING  TASKS: 


1.  XT  CALCOLATBS  THB  RBLBASB  OF  TARGET  AMD  PROJECTILE 
MATERIALS  FROM  SHOCKED  CCMIDITICMIS  DUE  TO  A  HYPERVELOCITY  IMPACT 
OF  A  MULTI-MATERIAL  PROJECTILE  ON  A  FLAT  THIN  TARGET  PLATE 
USING  THB  TILLOTSON  EQUATION  OF  STATE  TO  CALCULATE  THE  RELEASE 
XSENTROPE; 


2.  IT  CALCULATES  THB  RESIDUAL  MATERIAL  TEMPERATURES  FOR 
THB  TARGET  AND  PROJECTILE  MATERIALS; 


3.  IT  ESTIMATES  THB  PERCENTAGES  OF  THB  TARGET  AND  PRO¬ 
JECTILE  MATERIALS  IN  BACH  OF  THE  THREE  MATTER  STATES  BASED  0» 
THE  WASTE  HEAT  GENERATED  BY  THE  RELEASE  PROCESS; 


4.  IT  CALCULATES  THB  AMOUNT  OF  SOLID,  LIQUID,  AMD  GAS¬ 
EOUS  MASS  IN  THB  PROJECTILE  AMD  TARGET  MATERIAL  C(»ITRIBUTI<»fS 
TO  THE  DEBRIS  CLOUD  CREATED  IN  A  HYPERVELOCITY  IMPACT;  AMD, 


5.  IT  CALCULATES  THE  DEBRIS  CLOUD  VELOCITIES  AMD  THB 
SPREAD  OF  THE  DEBRIS  CLOUD  MATERIAL. 


.  THE  TILLOTSCm  EQUATICBl  OF  STATE  USED  BY  THIS  PROGRAM  INCLUDES 
.  THB  MIXED  PHASE  EQUATICXIS,  THE  CHECK  AT  V-VS,  A  CHOICE  OF  TWO 
.  ADJUSTMENTS  TO  THB  EQUATX(»(  OF  STATE  TO  ELIMINATE  THE  DISOCNI- 
.  TXMUITY  AT  V-VS,  A  CHOICE  OF  WHICH  HOLE  DIAMETER  EQUATKM  TO 
.  USE,  AND  A  CHOICE  OF  WHICH  dE— PdV  APPROXIMATION  TO  USE. 

IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z) 

DOUBLE  PRECISION  IMEP,IVEP,IMBT,IVBT,XSP,XST,XP,XT,LO,LP,LPTOT, 
$HUP, NUT, MPROJ, MTAR6, MUSTM, MTSR 
DOUBLE  PRECISIcm  OOPA(IO) ,RPA(10) ,KPA(10) ,BPA(10) ,ALPHPA(10) , 
$CPSPA(10) ,CPLPA(10) ,TMPA(10) ,TVPA(10) ,GPIA(10) ,HFPA(10) ,HVPA(10) , 
$BHNPA(10) ,ALFPA(10) ,LPA(10) ,BBTPA(10) ,EPSPA(10) ,SYPA(10) ,SUPA(10) , 
$NUPA(10) ,NUSPM(10) ,MPSR(10) ,FSRP(10) ,BXP(10) ,PMS(10) ,PML(10) , 

$PMV ( 10 ) , MPLYR ( 10 ) , UFSPl ( 10 ) , UFSP2 < 10 ) , PMSSR ( 10 ) , PMSSNR ( 10 ) , 
$UPPA(10) ,UPSA(10) ,CSP(10) ,VB(10) ,VD(10) ,VC(10) ,TDSUM(10) ,XE(10) , 
$TB(10),TD(10) 

INTEGER  ROPT,ROPTPA(10),ROPTT,HCOPT 
CHARACTER*!  PGMSTP 

CHARACTER*!  PIDA(IO) ,TID,PXDCHK,TIDCHK 
CBARACTBR*10  PNATA(IO) ,TMAT 

COMMON/TDATA/A, B, AA, BB, ALF, BBT, EO, EOM, EOI , ROPT, JCOPT 
OPEN (1, FILE- 'INDATA'  ) 

OPEN  ( 2 ,  FILE- '  IMPOUT ' ) 

.  READ  PROJECTILE  AND  TARGET  MATERIAL  PROPERTIES.  THE  PARAMETERS 
.  MUST  BE  IN  THE  FOLLOWING  UNITS: 


PID,TID  .  PROJECTILE  AND  TARGET  MATERIAL  ID  CODES 

PMAT,TMAT  .  PROJECTILE  AND  TARGET  MATERIALS 

COP, COT  .  ADIABATIC  BULK  SOUND  SPEED,  KN/S 

RP,RT  .  AMBIENT  MATERIAL  DENSITY,  6M/CU.CM. 

KP,KT  .  SLOPE  OF  US-UP  LINE,  DIMENSIONLESS 

BP,ET . ELASTIC  MODULUS,  LBS/SQ.IN. 

NUP,MUT  .  POISSON'S  RATIO,  DIMENSIONLESS 

GP,GT  .  AMBIENT  GRUNEISEN  COEFF. ,  DIMENSIONLESS 

ALFAP,ALFAT  ...  LINEAR  COEFF  OF  TERMAL  EXP,  1/DEG-C 

CPSP,CPST  .  SPECIFIC  HEAT  (SOLID),  CAL/GM-DEG-C 

CPLP,CPLT  .  SPECIFIC  HEAT  (LIQUID),  CAL/GM/DEG-C 
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TNPfTKT  . 

TVPfTVT  ....... 

H7P  f BVT  ....... 

RVP  t UVT  ....... 

BRNPfBHMT  . 

ALrPfALFT  . 

BKTPfBBTT  . 

XPSPfKPST  . 

SYPfSYT  . 

SUP, SOT  . 

IfP  ............ 

DP  ............ 

TS  ............ 

KOPTP,ROPTT  ... 
ROPT-1  .... 

ROPT«2  .... 

ROPT-3  .... 

HOOPT  . 

HCOPT-1  ... 
HCOPT-2  ... 
HOOPT«3  ... 
HOOPT-4  ... 
HC!OPT-S  ... 


MELT  TINPSRATUIUe,  DBG-C 
VAPORIZATION  TXMPXRATURZ,  0X6>C 
LATXMT  HBAT  OF  PUSI«I,  CAL/GN 
lATSNT  BXAT  OF  VAPCmiZATION,  CAL/6M 
BRIHXU.  HARDHBSS  NUMBER,  KG/SQ.MM 
TILLOTSCBI  EOS  OOHSTAMTS 
TILLOTSCBI  SOS  CXMSTAMTS 
TILLOTSCM  EOS  CCBfSTAMTS 
TENSILE  YIELD  STRENGTH,  MPA 
ULTIMATE  TENSILE  STRXMGTO,  MPA 
PROJECTILE  LAYER  LENGTH,  IN 
PROJECTILE  DIAMEnR,  IN 
TARGET  PLATE  THICKNESS,  IN 
TILLOTSON  SOS  RELEASE  OPTICM 
BACKNARD  PRESSURE  APPROXIMATION 
AVERAGE  PRESSURE  APPROXIMATION 
CURRENT  PRESSURE  APPROXIMATICNI 
HOLE  DIAMETER  EQUATION  OPTKXI 
KAPPII/HSSOl 
KAPPII/HSS02 
KAPPII/HSAOl 
PSN4/V10 

HOLE  DIA  -  PROJ  DIA 


IIRITI(*,1) 

1  FORMATC  ENTER  NUMBER  OF  PROJECv  .LS  MATERIAL  LAYERS  (12)  AND  HIT  R 
$XTURN*) 

RSAD(*,6)  NPMAT 

6  FCHU(AT(I2) 

DO  79  I«l, NPMAT 
NRZTS(*,3)  I 

3  F(ffiMAT('  ENTER  PROJ  MATL  ID  aX)B  FOR  LAYER  NO.  *,12,'  (A2)  AND  HIT 
$  RETURN') 

RXAO(*,5)  PIDA(I) 

5  FORNAT(A2) 

79  CONTINUX 
NRITS(*,7) 

7  FORMATC  ENTER  TARGET  MATERIAL  ID  CODE  (A2)  AMD  HIT  RETURN') 
REAO(*,9)  TID 

9  FORMAT(A2) 

DO  88  I«l, NPMAT 
REWIND  1 
RXAD(1,4) 

4  FORMAT(/////) 

99  READ(1,8)  PIDCHK 

8  FORMAT (A2) 

IF  ( PIDA( I ).8Q. PIDCHK)  THEN 

RSAD(1,10)  PMATA(I),COPA(I),KPA(I),RPA(I),GPIA(I),BHNPA(I) 

10  FORMAT(A10,5F10.5) 

RXAD(1,100)  EPA(I),NUPA(I),ALPHPA(I),CPSPA(I),CPLPA(I), 

$  SPSPA(I) 

100  FORMAT(2(E10.3,F10.5),2(F10.5)) 

RXAD(1,102)  TMPA(I),TVPA(I),HFPA(I),HVPA(I),ALFPA(1),BETPA(I) 

102  FORMAT (6F10. 5) 

RXAD(1,104)  SYPA(I),8UPA(I),ROPTPA(l) 

104  FORMAT(2F10.5,I1) 

XNDIF 

IF  (PIDA(I). ME. PIDCHK)  THEN 
IF  (PIDCHK.EQ. 'XX')  THEN 
WRITE  (*,17)  I 

17  FORMATC  PROJ  MATL  FOR  LAYER  NO.  ',12,'  NOT  FOUND  IN  MATERIAL  LIBR 
$ARY.',A'  please  check  PROJ  NATL  ID  CODES  AND  BEGIN  AGAIN.') 

STOP 


IMDIF 

ZF  (FiocaK.m. 'XX*)  thin 
HIM)  (1,2) 

2  FOm(HT(////) 
qozo  99 

IMDIF  • 

IMDIF 

•8  OOMTZMDI 
C 

RINIMD  1 
RIM>(1,4) 

999  RIM)(1,8)  TIDCHK 

IF  (TZD.IQ.TIDCHX)  THIN 
HIM)(1,10)  TMHT,00T,1CT,RT,GTI,BHMT 
RIM)(1,100)  IT,MDT,M.PHHT,CPST,CPLT,IPST 
RIM>(1,102)  TMT,TVT,HFT,HVT,HLFT,BITT 
RIM)(1,104)  8TT,8aT,HOPTT 
IMDIF 

ZF  (TID.MI.TIDCHR)  THIN 
IF  (TIDCHX.IQ. *XX')  THIN 
WRITI  (*,117) 

117  FORMHTC  TAROIT  MXTIRZAL  MOT  FOUND  IN  MHTIRIAL  LIBRARY.',/, 

$'  FLIASI  CHICR  TARGIT  MATIRZAL  ZD  CODI  AMD  BIGIN  AGAIN. ' ) 

STOP 

IMDIF 

IF  (TIDCaX.MI. 'XX' )  THIN 
HAD  (1,2) 

GOTO  999 

IMDIF 

INDIF 

C 

C .  RXAD  IMPACT  VILOCITY  IN  XM/S 

C 

1IRITI(*,29) 

29  FCHGfATC  INPUT  IMPACT  VILOCITY  IN  XM/SBC  (F5.2)  AMD  HIT  RETURN') 

RIAD(*,30)  V 

30  FORMAT(F5.2) 

C 

C .  RIAD  TARGIT  THICKNISS  AND  PROJECTILE  DIAMETER 

C 

1IRITE(*,13) 

13  FORMATC  ENTER  TS  AND  DP  VALUES  IN  INCHES  (F10.5,/,F10.5)  AND  HIT 
$RITURM') 

RIAD(*,11)  TS,DP 
11  FORMAT(F10.5,/,F10.5) 

C 

C .  RIAD  PROJICTILI  LAYER  THICKNESSES 

C 

DO  66  I-1,MPMAT 
1fRITE(*,44)  I 

44  FORMATC  ENTER  LP  VALUES  IN  INCHES  (F10.5)  FOR  LAYER  NO.  ',12, 

$*  AMD  HIT  RETURN') 

READ(*,33)  LPA(I) 

33  FORMAT(F10.5) 

66  CXHITZMUE 
C 

C .  READ  TARGET  HOLE  DIAMETER  EQUATION  OPTI(»l 

C 

flRITE(*,16) 

16  FORMATC  IMPUT  HOLE  DIAMETER  EQUATI(»(  OPTION  (II)  AND  HIT  RETURN', 

$/,5X, 'HOOPT-l  ...  KAPPII/HSS01',/,5X, 'HCOPT-2  ...  KAPPII/HSS02 ' , /, 

$5X, 'H00PT«3  ...  KAPPII/HSA01',/,5X, 'HCOPT-4  ...  PEN4/V10' , /,5X, 

$'H00PT-5  ...  DH-DP') 

READ(*,18)  HCOPT 
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18  rOIU<8T(Zl) 


....  KBM)  TILLOT80N  lOS  DZSOOMTZNUITY  ADJUSTMBNT  OPTION 
imZTB(*,21) 

21  POMATC  INPUT  lOS  DZSOONTINUITT  ADJUSTMKNT  OPTZ(»l  (II)  AND  HIT  RB 
$TUiaf',/,5X, 'JOOPT-1  ...  8CHCHIBXR0  JUMP  CX>ItRBCTI(MI' , /,5X,  ' JOOPT-2  . 
$..  MIXXO-PHASX  FORMULATION') 

RIAD(*,23)  JOOPT 
23  FORNAT(Il) 

....  80NX  PRSLIMIMART  CALCULATIONS  ... 

PI«4.0*ATAN(1.0) 

RPAVO«0.0 

LPTOT-0.0 

MPROJ-0.0 

DO  242  I-1,NPMAT 

RPAV6«RPAVG4LPA( I ) *RPA( I ) 

LPTOT-LPTOT-«>LPA(  I ) 

MPROJ«MPROJ-»-LPA(  I )  *RPA(  I ) 

242  CXMITINUB 

RPAVO-RPAVG/LPTOT 

MPROJ-PI*MPROJ* (DP/2 . 0 )* (DP/2 . 0) * (2 . 54*2 . 54*2 . 54 ) 

RBNIND  2 

IIRITS(2,40)  MPROJ,TMAT,V 

40  FORMAT (‘HTPBRVBLOCITY  IMPACT  OF  A  ‘,F8.3,‘  GM  MULTI -MATERIAL  PROJE 
$CTILX  (Ml  A*,/,A10, '  TARGET  AT  A  ',F5.2,'  XM/SEC  IMPACT  VELOCITY') 

WRITE(2,50)  TMAT,COT,KT,RT,TS*2.54,DP*2.54 

50  FORMAT( A 'TARGET  MATERIAL  PROPERTIES  . . . ' ,//,3X, 'MAT  -  ',A10,/,3?:, 

$'CO  «  ',F7.3,'  iai/S’,/,3X,'K  »  ' ,F7.3,/,3X, 'RHO  •  ',¥1.2, 

$'  GM/CU.CM. ',/,3X, 'TS  -  ',F7.3,'  CM' ,//, 'PROJECTILE  MATERIAL  PROP 
$ERTIES  (DP  -  ',F7.3,'  CM)  ...') 

DO  53  I«1,NPMAT 

NRITE(2,51)  I,PMATA(I),OOPA(I),KPA(I),RPA(I),LPA(I)*2.54 

51  FORMAT(A3X, 'MAT  ',12,'  ■  '  ,AiO,/,3X, '00  »  ',¥1.2,'  XM/S',/, 

S3X, 'X  -  ',F7.3,/,3X, 'RHO  -  ',F7.3,'  GM/CU.CM. ', /,3X, 

$'LP  »  ',¥1.2,'  CM') 

53  C»NTINUE 

....  CALCULATE  TARGET  MATERIAL  AND  1ST  PROJECTILE  LAYER  MATERIAL 
....  PARTICLE  AND  SHOCK  NAVE  VELOCITIES  AND  INTERFACE  HUGONIOT  IMPACT 
....  PRESSURE.  THIS  INFORMATION  IS  REQUIRED  TO  DETERMINE  THE  DEPTH 
....  TO  NHICH  THE  TARGET  MATERIAL  IS  SHOCXED  AND  TO  RELEASE  THE 
....  TARGET  MATERIAL  FROM  ITS  SHOCXED  STATE. 

IF  (TMAT.EQ.PMATA(l))  GOTO  35 
A-XPA( 1 ) -XT* (RT/RPA( 1) ) 

B-2 . 0*XPA( 1 ) *V+OOPA( 1 ) +COT* (RT/RPA( 1 ) ) 

C>COPA( 1 ) *V+XPA( 1 ) *V*V 

D«B*B-4.0*A*C 

UTP- (B-SQRT (D ) ) / ( 2 . 0*A) 

GOTO  38 
35  DTP-V/2.0 
38  UPPA(1)-V-UTP 
UTS«COT+XT*OTP 

UPSA  ( 1 ) -COPA  ( 1 )  •••XPA  ( 1 )  *UPPA  ( 1 ) 

PT«RT*0TS*0TP 

....  CALCULATE  SHOCX  NAVE  AND  PARTICLE  VELOCITIES  IN  SUBSEQUENT 
....  PROJECTILE  LAYER  MATERIALS 
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If  (MFIiAT.Gai.2)  THBM 
DO  85  X-2,MPN8T 

CDU.  UPrCM.(RPA(X>l)  ,RPA(X)  «KPA(X>1)  ,KPA(X)  ,CX»A(Z-1)  fCXXPAd) , 

$  UPfA(I-l),UPPA(I)} 

IIMA(  X  ) -OODAC  X  ) -t-KPAt  X  )  •OPPA(  X  ) 

88  OOMTXIIDB 
IMDXP 
C 

C .  CAX.COZATI  LOCATXON  XN  PROJBCTXLI  WBEAM  (mXOXNAI.  TAAOtT 

C .  RAMtPACTX<W  NAVI  OVIiaAXBS  PR0JICTXX.8  SHOCK  NAVI 

C 

CAZX  R8XMT(MPNAT,OT8,aTP,OP8A,tlPPA,I.PA,DP,T8,V,LO) 

C 

C .  TAIKSIT  MATIRXAL  RILIASI  CALCULATXON  PHASI 

C 

NRXTI(2,599) 

599  PORNAT( A  ******  TARGIT  MATIRXAL  RILIASI  CALCOLATION  *«***•) 

C 

VTO-l.O/RT 
VTl«RT*OTS/  (OTS-tJTP) 

VTl-l.O/VTl 

PR-PT*1.0I09 

IBT*0 . 5*PH* ( VTO-VTl ) / 1000 . 0 
PHMB>PH/101 . 3I<H)9 
C 

NRXTI(2,60)  aTP,aTS,PT,PHNB,IHT,VTO,VTl 
60  TC»U(AT(/, 'XMXTXAL  0(lfDXTX<»l8  FOR  TARGIT  MATIRXAL  . . .  * , /,3X, 'PARTXC 


$LI  VILOCXTT . OP  «  *,18.3,'  XM/S*  r/,3X, 'SHOCK  NAVI  8PIID 

$  .  08  •  ',F8.3, *  XM/S*,/,3X, ‘BOGONXOT  XMPACT  PRI8S0RI  . 


$...  PR  >  *,18.3,*  OPA  •  *,P6.3,'  M8AR’,/,3X, 'HOGONXOT  XMPACT  IMXRO 

$y  .  IH  •  *,110.4,*  JOOLI8/XO*,/,3X, 'SPICXPXC  VOLOMI  AT  RI8T  . 

$....  VO  •  *,P10.3,*  CO. CM. /GM*,/,3X, 'SPICXPXC  VOLOMI  AT  XMPACT  ... 
$  VI  •  *,P10.3,*  CO.CM./GM*) 

C 

C .  CALCOLATI  AMBXIMT  GROMIISIN  OOlPFICIINT  AND  GAMMA/SP.VOL.  RATIO 

C .  FOR  TARGIT  MATIRXAL 

C 

IT-IT*68947.0 
8ITAT-3 . 0*ALPHAT 
XF  (NOT.LT.0.5)  THIN 
K8T«IT/3 . 0/ ( 1 . 0-2 . 0*NOT) 

OOTC«DSQRT( (KST/10.0)/(RT*1000.0) )/1000.0 
INDXF 

XF  (NOT.IQ.0.5)  THIN 
KST— 1.0 
OOTO-1.0 
INDXF 

XF  (NOT.LT.O.S)  0T«2 . 388SI-08*K8T*BITAT/CPST/RT 
XF  (NOT.IQ.O.S)  GT«GTX 
C 

NRXTI(2,7S)  IT/10. 0,N0T,KST/10.0,ALPKAT,CPST,CPLT 
75  FORMAT (/, 'PARANITIRS  RIQOIRXO  FOR  CALCOLATING  TARGIT  MATIRXAL  RIUC 
$ASI  FROM  SH0CKID*,/,*8TATI  OSXNG  THI  TZLLOTSCHI  IQDATI(»(  OF  STATS t ' 

$,/,3X, 'TARO  MATL  ILASTXC  MODOLOS  .  I  -',110.4, 

$*  N/SQ.M.  *,/,3X, 'TARO  MATL  POXSSON  RATIO . NO 

$F10.3,/,3X, 'TARO  MATL  BOLK  NCDOLOS  .  K  -',110.4, 

$'  M/8Q.M. *,/,3X, 'TARO  MATL  LXN.  COIF.  OF  THBRM.  IXP.  ...  ALFA  -', 

$110.4,'  /OIO-C',/,3X, 'TARO  MATL  SP  BRAT  (SOLID)  .  CPS 

$  »',F10.3,'  CAL/GM/OIO-C*,/,3X, 'TARO  MATL  SP  HIAT  (LIQOID)  . 

$ .  CPL  -*,F10.3,'  CAL/(m/DIO-C') 

NRXTI(2,80)  6T,6TX,SyT,S0T,BHNT,TMT,TVT,RFT,HVT 
80  F(»MAT(3X, 'TARO  MATL  AMB  M-ORON  COIF  (CAL,INP)  ...  GAMO  -',F10.3, 

$',',F6.3,/,3X, 'TARG  MATL  TXILD  STRBN6TH  . .  SY  -', 

$F10.3,'  MPA',  A3X, 'TARG  MATL  OLT  STRENGTH .  SO 
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$,ri0.3, '  MPA *,/,3X, 'TARO  MATL  BRM  HDMS  MO .  BHM  - 

$*,P10.3,/,3X, 'TARO  MATL  MILT  TDIPIRATURI  .  TM  -',P10. 

$2«'  Oie-C',/,3X, 'TARO  MATL  VAP(Ml  TIMPIRATORI . TV' 

IPIO.2,'  I»0-C',/,3X,*TAR0  MATL  RIAT  OP  PU8Z0H . BP 


$«'«P10.2,'  CAL/QM',/,3X,'TARO  MATL  HEAT  OP  VAPORIZATION . H 

fV  >',P10.2,'  CAL/OM') 

C 

SBST>CP8T*4186. 

SBLTi«PLT*4ia6. 

Bf«"HPI*4ie6. 

HVTi«VT*41S6. 

IMIT>TMT*8B8T 

IVIT«XMBr»-HPT4-  (TVT-TMT)  *8HLT 
C 

MRITB  <2,76)  IMIT,IVIT 

76  PafyiAT(3X, 'TARO  MATL  IRICPIIMT  MILT  INUGY . IMI 

$110.4,  '  JOOLI8/XO',/,3X, 'TARO  MATL  IMCIPIIMT  VAPOR  IMIROy  . 

$IVI  ■',110.4,'  JOOLU/KO') 

C 

MRZTI  (*,230) 

230  P0RMAT(/,1X, 'IMTIR  lO  MOLTZPLZIR  VALUE  FOR  TARGET  MATERIAL  (F4.2) 
$AMD  HIT  RETURN' ) 

READ  (*,240)  lOM 
240  F01IMAT(F4.2) 

ALF-ALFT 
EET»BETT 

ROPT-RCWTT 

GALL  T00N8T(VT0,00T,XT,0T,TNT,TVT,BFT,BVT,IVET) 

CALL  TRIL81(VTO,VT1,PH,IXT,UTP,ORT,UF8T1,UF8T2,IVET,HVT,OOT,XT, 

$  EP8T) 

CALL  TINC(8B8T,8HLT,TMT,TVT,HFT,HVT,IXT,IMIT,IVET,PS,PL,PV,TRT) 
NRITI(2,87)  0F8Tl,UF;'iT2 

87  F0RMAT(A'FREI  8URF  VEL  (UP-t-UR)  . ',F7.3,’  KN/SEC ' , / , ' FREE  S 

SURF  VIL  (2.0*UP)  . ',F7.3,'  KN/8EC’ ) 

CALL  TNCALC(V,P8,PL,PV,RPAVO,RT,TRT,TNT,TVT,TS,LPTOT,OP,BHMT, 
$T80LT,TM8,TML,TNV,NTAR6,H00PT,SYT,8UT,UT8,UPSA(1),UTP,UPPA(1), 
$MU8TN,MT8R) 

...  PROJECTILE  MATERIALS  RELEASE  CALCULATIONS  PHASE 
CHKL*0.0 

DO  9999  I«1,NPMAT 
NRITI(2,89)  I 

89  FOIOIAT( A '*****  PROJECTILE  MATERIAL  RELEASE  CALCULATICNfS,  LAYER  NO 

$.',12,'  ***♦*') 

...  READ  MATERIAL  PARAMETERS  FOR  CURRENT  PROJECTILE  LAYER  MATERIAL 

OOP^COPA(I) 

RP«RPA(I) 

KP«XPA(I) 

IP«IPA(I) 

ALPHAP«ALPHPA(I) 

CPSP^CPSPA(I) 

CPLP*CPLPA(I) 

TMP^TMPA(I) 

TVP-TVPA(I) 

OPZ-GPIA(I) 

HFP-BFPA(I) 

BVP*HVPA(I) 


ooo  u  oouu 


MOff-mVAd) 

M.vr<4a.rrA(X) 

x.»-Lm(X) 

■ItfWCTi^d) 

i»*r-iMra(X) 

n»««TPA(X) 

■IIP-tllVA(X) 

MOP-IRirA(X) 

IK»TP-1I0VTPA(X) 

niAT«fllA1A(X) 

IIWa4IPPA(X) 

UM>UMA(X) 

CHXLKSXL't-X.PA  ( X )  *2 . 54 

....  CA[«aXATI  880000  8TATB  0IA1ITXTXB8  FOR  PROJBCTXLK  LAYIR  HATRRXAL 

PP>RP*UPS*OPF 
ra-PF*l. 01409 
FOMI-ra/lOl .  31409 
VFO-l.O/RP 
V»1«RP*1IP8/  (0P8-0PP) 

VP1>1.0/VP1 

IHP>0. S*FH* (VPO-VPl) /lOOO.O 

wtxn(2,93)  x,app,aM,pp,raMB,iHP,vpo,vpi 
93  fOIOIAT(/,'XMXTXAL  GOMOXTXOM8  FOR  PROJICTXLB  LAXRR  NO.  ',X2,'  MATIR 

$XAL  ... ',/,3X, 'PARTXCXA  VRLOCXTY  .  OP  -  ',98.3, •  XM/8',/ 

$,3X, ‘8H0aC  NAVI  XPXID  .  08  >  ',F8.3,'  XM/8* ,/,3X, '8B0CK 

$  PRX88aRS . PB  -  •,F8.3,*  GPA  •  ‘F6.3,'  NBRR’,/,3X, 

$*8B0aC  IMBROy . IB  -  M10.4,  *  J00LX8/K0',/,3X, 

|*8PICXPXC  VOAOMB  (AT  RX8T)  ...  VO  >  SP10.3,'  CO.OI./GN’ ,/,3X, 
I'SPBCXFXC  VOtOMB  (BBOOCXD)  ...  VI  >  ’,fl0.3,’  CD.OI./OM’) 

....  CAIXOLATI  AMBXBRT  GR0NIX8XN  CXHtPFXCXXNT  AND  GRNNA/8P.VOL.  RATXO 
....  FOR  PROiJBCTXXO  NATXRXAL. 

BP«BP*68947.0 
nTAP>3.0*ALPBRP 
XF  (NOP.X.T.O.S)  TBXB 
K8P-BP/3 . 0/ ( 1 . 0-2 . 0*N0P  ) 

OOPOD80tT(  (K8P/10.0)/(RP*1000.0)  )/1000.0 
BNDXF 

XF  (MUP.BQ.0.5)  TBBN 
K8P— 1.0 
OOPC>-1.0 

BNDXF 

XF  (NOP.LT.O.S)  0P*2.388SB-08*K8P*BBTAP/CP8P/RP 
XF  (NDP.BQ.O.S)  0P>0PX 

1IRXTB(2,10S)  X, BP/10. 0,NUP,K8P/10.0,ALPBAP,CPSP,CPLP 
105  FaRNAT(A*PARRMBTBR8  RBQDXRBD  FOt  CALCULATING  RXLBASB  OF  PROJ  LAYS 
$R  NO.  ',12,'  NATBRXAL  FROM* ’8B0CXBD  8TATB  USING  TBB  TILLOTSON  B 


SOUATXON  OF  8TATBt’,/,3X, ’NATL  BLASTXC  NOTULUS  .  B 

$',B10.4, '  N/SQ.N. ',/,3X,'MATL  POISSON  RATIO  .  NO  « 

$',F10.3,/,3X, 'NATL  BULK  MODULUS  .  K  -',B10.4, 

$'  N/8Q.N. ',/,3X, 'NATL  LXN.  COTF.  OF  TBBRM.  BXP.  ...  ALFA  -',B10.4, 

$*  /DBG-C',/,3X, 'NATL  8P  BBAT  (SOLID)  . CPS  -',F10.3, 

$'  CAL/GM/DB8-C',/,3X, 'NATL  SP  BBAT  (LIQUID)  . CPL 

$F10.3,'  CAL/GM/DBO-C') 


1IRXTB(2,110)  GP,GPX,8YP,8UP,BBNP,TNP,TVP,BFP,BVP 
no  F0IINAT(3X, 'NATL  AMB  M-GRUN  COBF  (CAL,XNP)  ...  GAMO  ,F10.3, ' , ' , 

$F6.3,/,3X,'IIATL  TXBLD  STRBNGTB  .  SY  -',F10.3,'  MPA' 

$,/,3X, 'NATL  ULT  STRBNGTB  .  SU  »',F10.3,'  MPA',/, 

$3X, 'NATL  mw  BDN8  MO .  BBM  -' ,P10.3,/,3X, 'MATL  MB 


c 


c 


c 


c 


$LT  TINPSMTOKI . TN  -’,ri0.2,'  DKG-C  ,/,3X, ‘NATL  VAPO 

$»  TBmKATOXI . TV  DIO-C' ,/,3X, 'MAn  HEAT  O 

$T  fl»ZOH . HP  **,P10.2,‘  CAL/GN',/,3X,  *MATL  HEAT  OP 

$  VAPORISATION  .  HV  •• ,T10.2,'  CAL/GN') 


EBEP-CPSP*41E«. 

m.P-CTLP*4ie6. 

HPP-RPPM186. 

BVP>HVP*4186. 

IiapaTIIP*8B8P 

IVEP-IMEP4>BPP+(TVP-TNP)  *8HLP 


NEin  (2,77)  IMEPflVEP 

77  P0R11AT(3X, 'NATL  IMICPIEMT  MELT  EMEEGY .  XME  -‘,E10.4,  '  JOG 

SLE8/XO‘,/,3X, 'NATL  IMCZPIEMT  VAP(»  ENERGY  .  IVE  -',E10.4,'  J 

lOOLES/XO*) 


UNITE  (*,231)  I 

231  POraiAT(/, IX, 'ENTER  EO  MOLTIPLZER  VALUE  POR  PROJECTILE  LAYER  NO.  ', 
$12,'  MATERIAL  (P4.2) ' ,/,lX, 'AMD  BIT  RETURN') 

READ  (*,241)  EON 
241  P<NUIAT(P4.2) 


ALP-ALPP 

BET«BETP 

ROPT-ROPTP 

CALL  TOORST(VPO,COP,KP,OP,TMP,TVP,HPP,HVP,ZVEP) 

CALL  TRBL81(VPO,VPl,PH,EXP(I),UPP,URP,UPSPl(I),UPSP2(X),IVEP, 

$  BVP,OOP,XP,EPSP) 

CALL  TIMC(8H8P,SHLP,TMP,TVP,BPP,BVP,EXP(I),ZMEP,IVEP,PS,PL,PV, 

$  TRP) 

CALL  PNCALC(Z,UP8,UTS,UPP,UTP,RP,PS,PL,PV,T8,LP,DP,PMS(Z) , 

$  PM88R(X),miL(Z),PMV(I),MPLYR(Z),NUSPM(Z),NPSR(Z), 

$  P8RP(Z),POMSTP,ISTOP,LO,CHXL) 

PNS8RR(Z)>0.0 

C 

IP  (PGNSTP.EQ. 'Y')  GOTO  498 
C 

9999  CONTINUE 


498  NRITE  (2,499) 

499  PORMAT(/^ 'NASS  DISTRIBUTKW  SUMMARY  ...’) 

DO  501  I«1,RPMAT 

IP  (PGMSTP.EQ. 'Y'.AND.I.GT.XSTOP)  THEN 
PNS(I)«PI*(OP/2.0)*(DP/2.0)*(LPA(I)*2.54)*RPA(Z) 

PNL(I)«0.0 

PNV(I)-0.0 

MUSni(I)«PMS(Z) 

PMSSR(Z)>0.0 

PNSSMR(Z)-0.0 

MPSR(I)-0.0 

PSRP(Z)-0.0 

EMDIP 

NRITE  (2,500)  I,PMS(I) ,MUSPM(I) ,PMSSNR(I) ,PMSSR(Z) ,PML(I) ,PMV(I) 

500  P(ffiNAT(3X, 'PROJECTILE  LAYER  NO.  ',12,'  ...  SOLID  -  ',F7.2, 

$'  GMS',/,32X, 'UMSH  ....  ',P7.2,'  GMS' ,/,32X, 'SNR  .  ',P7.2, 

$'  OMS',/,32X, 'S8R  .  ',P7.2,'  GMS' ,/,31X,' LIQUID  ...  ',P7.2, 

$'  GMS ',/,31X, 'VAPOR  ....  ',P7.2,'  GMS') 

501  COMTIMUE 
TMPRAG-MTARG-MTSR 

NRITE  (2,502)  TSOLT,TMPRAG,TMS,TML,TNV 

502  PORMAT(3X, 'TARGET  MATERIAL  .  SOLID  -  *,F7.2,'  GMS*,/, 

$32X,'PRAG  ....  ',P7.2,'  GMS' ,/,32X, 'S6R  .  ',F7.2,'  GMS',/,31X, 

$'LIQUID  ...  ',P7.2,'  GMS ',/,31X, 'VAPOR  ....  •,F7.2,*  GMS*) 
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....  COMVOn  DIBRZ8  CLOOD  VnOCITXBS 

CaOL  DCVIL(0F8Tl,arSPl,V,lffAM,MPlte>J,lfPLYR,KT8R,MP8R,M(JSTM, 
$  Mn8m,PM881IR,BXT,BXP,P8RP,llPMRT) 


CX08I(1) 

CtJ0n(2) 

CL088(3) 

8tW 

BHD 


RUBROmiNB  TCiaH8T(V0,00«K,0,1M,TV,RP,HV,R8) 

IMPLICIT  DOOBLB  MaK!I8ZON  (A>H,0-S) 

OOOBLB  PRBCI8I0M  K 
THTBCIB  ROPT 

OOMNOM/TDATA/A, B, AA, BB, ALP, BBT,BO, BON, BOI , ROPT, JOOPT 

...  THZ8  8UBROOTII1B  CALC0LATB8  THB  VALUBS  OP  THB  OmSTAllTB 
...  RBOCIRBO  BT  TBB  TILLOTSON  BQOATION  OP  STATE  (8wRI  PIKAL 
...  RBPCWT  worn  PROJ.  MO.  06-4438). 

AA- ( 1000 . O/VO) * (00*1000 . 0 )* (C0*1000 . 0 ) 

BB^AA* (2 . 0*K-1 . 0-0. 5*0) 

A-O.S 

B-O-O.S 

Rl-TM/TV 

R2-HP/HV 

BOZ>BXP(-0. 199) * (K**6. 5939) * (R2**0. 5720) / (0**0.7680) 

$  /(Rl**0.0210) 

BOI«BOI*(B8*HV) 

BO-BOM*B01 

WRZTB  (2,10)  AA,BB,A,B,ALP,BBT,BOI,BON,BO 
10  POra(AT(/, 'AINDZTXOMAL  PARANBTBR8  RBQDIRBD  POR  CALCULATING  MATERIAL 
SRBLBASE  PROM*,/, 'SHOCKED  STATE  USING  THB  TZLLOTS<Nf  EQOATKXI  OP  STA 
$TBl*,/,3X, *AA  «  ',E11.4,'  M/8Q.N. ',/,3X, *BB  -  *,X11.4,*  N/SQ.M. ' 
$,/,3X,'A  •  ',P7.4,/,3X, *B  -  ' ,P7.4,/,3X, 'ALP  »  •,P7.4,/,3X, 

$*BET  •  ',P7.4,/,3X, 'BOX  «  *,E11.4,’  JOULBS/KG' ,/,3X, 'EOM  -  ',P7.4, 
$,/,3X, 'BO  »  ’,B11.4, '  JOULBS/KG’) 


RBTORM 

no 

SUBROOlIMB  TRXLS1(VO,V1,PHO,BX,UP,UR,UPS1,UPS2,ZVB,HV,CO,K,BPS) 
IMPLICIT  DOUBLE  PRECISION  (A-H,0-X) 

DOUBLE  PRBCISKNI  Q(401) ,NU(401) ,V(401) ,B(401) ,P(401) ,R(401) 

DOUBLB  PRECISION  RP(401) ,U(401) ,S(401) ,PH(401) , IVB,K 
INTEGER  ROPT 

OQMNOH/TDATA/A, B , AA , BB , ALP , BBT, EO, BOM , EOI , ROPT , JCOPT 

..  THIS  SUBROUTINE,  TOGETHER  WITH  THE  SUBROUTINE  PCALC,  CALCULATE 
..  THB  RELEASE  OP  THE  PROJECTILE  AMD  TARGET  MATERIALS  USING  THE 
..  TZLLOTSCNI  EQUATION  OP  STATE.  XT  IS  ASSUMED  THAT  FOR  MOST  METALS 
..  THB  SPECIFIC  VOLUME  VS  IS  APPROX.  13.1%  GREATER  THAN  THE  AMBIENT 
..  SPECIFIC  V(n.UMB  VO. 

ESP-ZVE-i-HV 

V8«1.131*VO 

WRITE  (2,5)  ZVE,HV,ESP,VS,EP8 

5  FORMAT ( 3X, *ES  •  ',E11.4,’  JOULBS/KG' , /,3X, 'HV  -  ',B11.4, 

$*  JOULBS/KG ',/,3X, 'BSP  •  ',B11.4,’  JOULBS/KG’ ,/,3X, 'VS  >  ', 

$P7.4, *  CU.CM./0M',/,3X, 'BPS  -  ',P7.4) 
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PR(l)>rao 

V(1)»V1 

F(i)-rao 

I(1)«0.5*P(1)*(VO-V1)/1000.0 
DBLV-(VO-V1)/50.0 
1IO(1)«VO/V(1)-1.0 
Q  ( 1 )  >AA*iai  ( 1 )  ( 1 )  *Ma  ( 1 ) 

R ( 1 ) >IO* ( VO/ 1000 . 0 ) * ( VO/ 1000 . 0 ) 

RP(1)-B*R(1) 

....  MOTBt  MO(l),Q(l),R(l)«RP(l)  ARX  INXTIALZZBD  BUT  NOT  USED 

PBQ2-0.0 

DBLP-0.0 

DB-0.0 

UR-0.0 

ZZ-0 

DO  10  Z«2,401 
V(Z)-V(Z-1)+DBLV 

PH{ Z ) -CO**2* { 1000. O/VO) * ( 1. 0-V< Z ) /VO) / ( 1 . O-K* ( 1 . 0-V( Z ) /VO) ) **2 

Pa(Z)«PH(Z)*1.0B06 

MU(Z)-VO/V(Z)-1.0 

R(Z)-B(Z-1)*(V(Z)/1000.0)*(V(Z)/1000.0) 

$  <fBO*(VO/1000.0)*(VO/1000.0) 

ZF  (V(Z).LT.VO)  THBN 

Q  (  Z  ) -AA*M0  (  Z  )  >00*110  (  Z  )  *110  (  Z  ) 

RP(Z)-A*B(Z-1)*(V(Z)/1000.0)*(V<Z)/1000.0) 

$  >B*BO*(VO/1000.0)*(VO/1000.0) 

CALL  PCALC(B(Z-l),P(Z-l),V(Z),Q(Z),R(Z),RP(Z),VO,DBLV,P(Z)) 
BNDZF 

ZF  (JOOPT.BQ.l.AND.V(Z).GE.VO)  THEN 

....  ZMPLEMBNTATZON  OF  SCHONBERC  JUMP  CORRBCTZON  (NL-TR-93-7028) 

. . . .  TOGETHER  NZTH  MZXED  PHASE  FORKULATZON 

ZF  (V(Z).LT.VS)  THEN 

ZF  (B(Z-l).LT.ZVE)  THBN 

g(Z)-AA*MU(Z)>BB*MU(Z)*MU(X) 

RP(Z)-A*B(Z-1)*(V(Z)/1000.0)*(V(Z)/1000.0) 

$  >B*EO*(VO/1000.0)*(VO/1000.0) 

CALL  PCALC(E(Z-l),P(Z-l),V(Z),Q(Z),R(Z),RP(Z),VO,DELV,P(Z)) 
ENDZF 

ZF  (E(Z-’1).GE.ZVE.AND.E(Z-1).LT.ESP)  THEN 
Q  (  Z  )  «AA*NU  (  Z  )  >00*101  (  Z  )  *NU  (  Z  ) 
RP(Z)-A*B(Z-1)*(V(Z)/1000.0)*(V(Z)/1000.0) 

$  >B*BO*{VO/1000.0)*(VO/1000.0) 

CALL  PCALC(B(Z-l),P(Z-l),V(Z),Q(Z),R(Z),RP(Z),VO,DELV,PC) 

C»V(Z)/VO-1.0 

U ( Z ) -DEEP ( -ALF*C*C ) 

S ( Z ) -AA*MU ( Z ) *DEXP ( -BET*C) 

B«B*U(Z) 

Q(Z)-U(Z)*S(Z) 

RP(Z)-A*E(Z-1)*(V(Z)/1000.0)*(V(Z)/1000.0) 

$  >B*EO*(VO/1000.0)*(VO/1000.0) 

CALL  PCALC(E(Z-l),P(Z-l),V(Z),Q(Z),R(Z),RP(Z),VO,DELV,PB) 
B-B/U(Z) 

T1-PB*(E(Z-1)-ZVB) 

T2-PC*(BSP-B(Z-1)) 

OEN-ESP-ZVE 
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oou 


rcz)«(Ti-»^T2)/Dni 

iMDir 

c 

If  (I(X-2).6I.I8P)  TBni 

C-V(I)/VO-1.0 

n  (  X  )  •DBXP  ( >JUUP*C*C  ) 

8  (  X  ) -AJk*llD  ( X )  *DBXP  ( -Bn*C) 

B«8*U(X) 

Q(X)-0(X)*8(X) 

ltP(X)-A*B(X>l)*(V(X)/1000.0)*(V(X)/1000.0) 

$  ••>B*10*(VO/1000.0)*(VO/1000.0) 

CKLL  PCBLC(B(X~l),P(X-l),V(X),Q(X),R(X),RP(X),VO,DBLV,P(X)) 
B«B/U(X) 

BMDXP 

C 

DBLV8«V8-V(X) 

XF  (DBLVS.lt. OBLV)  THBN 

C»V(X)/VO-1.0 

D ( X ) -DBXP ( -ALP*C*C) 

8 ( X ) -AB*ND ( X ) *OBXP ( -BBT*C ) 

B>B*U(X) 

Q(X)-0(X)*8(X) 

RP(X)«A*B(X-1)*(V(X)/1000.0)*(V(X)/1000.0) 

$  4B*BO*(VO/1000.0)*(VO/1000.0) 

CALL  PCMLC(B(X-l),P(X-l),V(X),Q(X),R(X),RP(X),VO,DBLV,PBQ2VS) 
B«B/U(X) 

DBLP-PBQ2V8-P(X) 

BMDXP 

C 

BMDXP 

C 

XP  (V(X).GB.VS)  THBM 

C»V(X)/VO-1.0 

0 ( X ) -DBXP ( >ALP*C*C ) 

8  ( X  )  •RR*IID  ( X )  •DBXP  ( -BBT*C) 

B«B*U(X) 

Q(X)«U(X)*S(X) 

RP(X)«A*B(X-1)*(V(X)/1000.0)*(V(X)/1000.0) 

$  ■••B*BO*(VO/1000.0)*(VO/1000.0) 

CALL  PCALC(B(X-l),P(X-l),V(X),Q(X),R(X),RP(X),VO,DBLV,PBQ2) 
B>B/D(X) 

P(X)«PBQ2>DBLP 

BMDXP 

C 

BMDXP 

C 

XP  (JOOPT.BQ.2.AMD.V(X).GB.VO)  THBM 

XMPLBNBNTATXOM  OF  PDRB  NXXBD  PHA8B  PORMDLATIC8I 


XP  (B(X>1).LT.XVB)  THBM 

Q  (  X ) -AA*MD  ( X  )  •l'BB*N(l  ( X )  *MU  (  X  ) 

RP(X)»A*8(X-1)*(V(X)/1000.0)*(V{X)/1000.0) 

$  •t-B*BO*(VO/1000.0)*(VO/1000.0) 

CALL  PCALC(B(X-l),P(X-l),V(X),Q(X),R(I),RP(I),VO,DBLV,P(I)) 
BMDXP 

XP  (B(X-1).GB.XVB.AMD.B(X-1).LT.BSP)  THBN 

Q(X)«AA*MD(X)+BB*MU(X)*MU(X) 

RP(X)>A*B(X-1)*(V(X)/1?00.0)*(V(X)/1000.0) 

$  +B*BO*(VO/1000.0)*(VO/1000.0) 

CALL  PCALC(B(X-l),P(X-l),V(X),Q(X),R(X),RP(X),VO,DBLV,PC) 

C»V{X)/VO-1.0 

U ( X ) -DBXP ( -ALP*C*C ) 
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8  (  Z ) -AA*iai  ( X  )  *DBXP  ( >UT*C ) 

B<>B*U(Z) 

Q(I)-0(I)*8(I) 

»P(I)-A*1(I-1)*(V(1)/1000.0)*(V(I)/1000.0) 

$  •»-B*IO*(VO/1000.0)*(VO/1000.0) 

CBZX  BCAZ.C(B(Z-l),P(Z-l),V(Z),Q(Z),R(X),RP(Z),VO,DBLV,PB) 

B-B/U(X) 

T1-PB*(B(X-1)>XVB) 

T2«K!*  (B8P-B  ( X-1 ) ) 

DBR«B8P*XVB 

P(X)-(Tl+T2)/Dm 

BMDXF 

C 

XP  (B(X-l).GB.BSP)  THBN 

C-V(X)/VO-1.0 

U  (  X  ) -DBXP  (  >AZJr*C*C ) 

8 ( X ) -AA*MU ( X ) *DBXP ( -BBT*C ) 

B-B*U(X) 

Q(X)-0(X)*S(X) 

RP(X)>A*B(X-1)*(V(X)/1000.0)*(V(I)/1000.0) 

$  •t>B*BO*(VO/1000.0)*(VO/1000.0) 

CALL  PCALC(B(X-1),P(X-1),V(Z),Q(Z),R<Z),AP(Z),VO,OBLV,P(Z)) 
B>B/a(X) 

BHDXP 

c 

BHDXP 

C 

C .  CALCDLATB  BNBRGXBS  BASED  <Xi  RBLBASB  APPROXXMATI(»l  OPTXON 

C 

XP  <ROPT.BQ.l)  THEN 

E  ( X  )  <-E  ( X-1 ) -P  ( I>1 )  *DELV/ 1000 . 0 

EHDXP 

XP  (ROPT.EQ.2)  THBH 

E(X)«E(X-l)-0.5*(P(X-l)-fP(X))*DELV/1000.0 

EHDXP 

XP  (ROFT.BQ.3)  THEH 
E(X)-E(X-1)-P(X)*DELV/1000.0 
EHDXP 
C 

DP-P(X)-P(X-1) 

XP  (DP.6E.0.0)  THEH 
HRXTE  (2,11)  X 

11  PORMATC***  AH  ERROR  HAS  OCCURRED  ZH  RELEASE  PROCESS  AT  THE  ',X3, 
$'-TH  XTERATIOH  •••') 

STOP 

BHDXP 

DUR-DSQRT(-DP* (DELV/1000. 0) ) 

UR-UR-t-DUR/ 1000 . 0 
XX-XX+l 

XP  (P(X).6E.0.0)  THEH 

XP  (ROPT.EQ.l)  DE«DE+P(Z-1)*DBLV/ 1000.0 
XP  (ROPT.EQ.2)  DE-DB-l-0.5*(P(Z)-l'P(X-l))*OELV/1000.0 
XP  (ROPT.EQ.3)  DE«OE-«>P(X)*DELV/1000.0 
ADP»OABS(DP) 

DPR-ADP/  ( P  (  X-1 )  <«-DELP  ) 

XP  (DPR.LT.EPS)  GOTO  15 
EHDXP 

XP  (P(X).LT.O.O)  GOTO  15 
10  0(»1TXHUE 
C 

15  EX«E(1)-DB 
VP»V(XX) 

UPS1-UP4-UR 

UPS2-2.0*UP 
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c 

imXTB(2,20)  V7,I(1),DB,BX 

20  r0niBT(/, 'BMD-STATB  CALCDIATX<»I  lUe8UI.TS  DSXNG  THB  TXIXOTSON  BOS  .. 

$.',/,*IIATBItXJa.  FXM  8P  VOL  (VP)  .  ’,P10.3,'  CU.CM. /ON' ,/, 

S’NATBRIAL  8HOCK  BNIBOY  .  *,110.4, ‘  JOULX8/KG* , /, 'MATBRXAL  B 

SMBBOY  RBCOVBRBD  ....  *,B10.4,*  JOOLBS/K6* , /, 'HASTB  HBAT  OBMBRATBO 

$ .  ',B10.4,'  JOOLBS/KO*) 

C 

RBTURM 

BHD 

C 

SUBBOOTINB  PCALC(B,P,V,Q,R,RP,VO,DBLV,PX) 

XMPLXCIT  DOUBLB  PRBCXSXOM  (A-B,0-Z) 

IBTBOBR  ROPT 

OOMMCm/TDATA/ A,  B  ,  AA,  BB  ,  ALP ,  BBT  ,  K>,  BOM,  BOX  ,  ROPT ,  JOOPT 
C 

C.....  CALCULATB  PRBSSURBS  BA8B0  ON  RBLBASB  APPROXXMATXON  OPTXON 
C 

XP  (ROPT.BQ.l)  THBN 
T1-B-P*DBLV/ 1000 . 0 
OTl-Tl/BO 

DBR«0T1* (V/VO) * (V/VO)+l 
Pl«  (  A-t-B/OBM )  *  ( Tl/  ( V/ 1000 . 0 ) ) 

BMDXP 

C 

C2P>0.0 

C3P-0.0 

XP  (ROPT.BQ.2)  THBN 
DBLV-OBLV/2.0 

C2P-P* (OBLV/1000 . 0 ) * (V/1000 .0) * ( V/1000. 0 ) * ( 1 , O+A* ( 2 . 0*DBLV/V) ) 
C3P1> ( 1 . 0+A) *B* (V/ 1000 . 0 ) * (V/1000 . 0 ) 

$  •«>(1.04-B)*BO*(VO/1000.0)*(VO/1000.0) 

$  ♦Q*(V/1000.0)* (V/1000. 0)*(V/1000.0) 

C3P-P* ( DBLV/ 1000 . 0 ) *C3P1 

$  - ( P* (DBLV/ 1000 ) ) • ( P* (DBLV/1000 . 0 ) ) • (V/1000 . 0 ) * (V/1000 . 0 ) 

BNDIP 
C 

XP  (R0PT.BQ.2.0R.R0PT.BQ.3)  THBN 

Cl- ( V/ 1000 . 0 )*( DBLV/ 1000 . 0 )*( 1 . 0+A* ( DBLV/ V ) ) 

C2-Cl*R/( (V/1000. 0)*(DBLV/1000.0) )-f(DELV/V)*RP 
$  ♦g* (V/1000. 0)* (V/1000. 0)*(DELV/1000.0)-C2P 

C3- ( A*B+Q* (V/ 1000 . 0 ) ) *R+B*B*K)* (VO/ 1000 . 0 ) * ( VO/1000 . 0 ) -C3P 
DISC«C2*C2-4 . 0*C1*C3 

PXl-  ( C2-»-DSQRT  ( DISC ) )  /  ( 2 . 0*  (V/1000 . 0 )  *C1 ) 

PI2- ( C2-0SQRT ( D tSC ) ) / ( 2 . 0* ( V/ 1000 . 0 ) *C1 ) 

PX-PZ2 

BNDXP 

C 

IP  (ROPT.BQ.2)  DELV-2.0*0ELV 
C 

RETURN 

END 

C 

SUB1M>UTXNB  TXNC(SHS,8HL,TN,TV,HP,HV,BXH,XMB,XVB,PS,PL,PV,TR) 
IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z) 

DOUBLB  PRECISION  IMB,IVB 


C 

C .  THIS  SUBROUTINB  CALCULATES  THB  RESIDUAL  TEMPERATURE  INCREASE 

C .  IN  A  MATERIAL  THAT  HAS  BEEN  RELEASED  PROM  THE  SHOCKED  STATE 

C .  AND  ESTIMATES  THB  PERCENTAGES  OP  VAPORIZED,  MELTED,  AND  SOLID 

C .  MATERIAL  DUE  TO  THE  RELEASE  PROCESS. 

C 

C .  IP  HASTE  HEAT  IS  LESS  THAN  THE  ENERGY  REQ'D  TO  START  MELT, 

C .  CALCULATB  TEMPERATURE  RISE  USING  W.H.-S.H. * (TEMP. INCR. ) 
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IF  (KXH.lt. IMB)  THBN 

DT-XXH/SHS 

TK-DT 

DKL-0.0 

1IRITK(2,S0)  ZMK,DKL,BXH 

50  FORMAT ( 'XNXROT  RXQ,  INCIPIENT  MELT  ...  ',E10.4,'  JOOLES/KG ' , / , 

$'KNBROy  AVAILABLE  FOR  MELT  -  SEIO.A,'  JOOLES/KG ' , / , 

$'BXCESS  ENERGY  AVAILABLE  .  ‘,E10.4,'  JOULES/KG') 

PV-0.0 
PL-0.0 
PS-100.0 
GOTO  100 
BNDIF 

...  IF  HASTE  HEAT  EXCEEDS  THE  ENERGY  REQ’D  TO  START  MELT,  BUT  IS 
...  LESS  THAN  THAT  RKQ'O  TO  COMPLETE  MELT,  RESET  THE  VALUE  OF  THE 
...  ENERGY  AVAILABLE  FROM  THE  HASTE  HEAT  VALUE  TO  THE  VALUE  RSQ*D 
...  TO  START  MELT.  THIS  IMPLIES  THAT  SOME  ENERGY  IS  AVAILABLE  FOR 
...  MELTING  A  PORTION  OF  THE  MATERIAL.  NOTE:  THE  TEMPERATURE  RISE 
...  EQUALS  THE  MELT  TEMPERATURE  OF  THE  MATERIAL. 

IF  (EXH.GE. IME.AND.EXH.lt. IME-l-HF)  THEN 
TR-TM 

DBL-BXH-7NE 

REQM-IME^HF 

NRITE(2,60)  IME,RB<H(«DEL 

60  FORMAT ( 'ENERGY  REQ,  INCIPIENT  MELT  ...  ',E10.4, '  JOULES/KG ' , / , 

$ 'ENERGY  REQ,  CCBIPLETE  MELT  ....  ',E10.4,'  JOULES/KG',/, 

S'EMERGY  AVAILABLE  FOR  MELT  ....  ',E10.4,'  JOULES/KG') 

PV-0.0 

PL-100. 0*DEL/HF 
PS-100. 0-PL 
GOTO  100 
ENDIF 

...  IF  THE  HASTE  HEAT  EXCEEDS  THE  ENERGY  REQ'D  TO  COMPLETELY  MELT 
...  THE  MATERIAL,  BUT  IS  LESS  THAN  THAT  REQ'D  TO  START  VAPORIZA- 
...  TION,  COMPUTE  THE  TEMPERATURE  INCREASE  CAUSED  BY  THE  EXCESS 
...  ENERGY  AND  ADD  IT  TO  THE  MELT  TEMPERATURE  OF  THE  MATERIAL. 

IF  (EXH.GE. IME+HF. AND.EXH.lt. IVE)  THEN 

DEL-EXH-IME-HF 

DT-DEL/SHL 

TR-TM+DT 

REQM-IME-I-HF 

NRITE(2,70)  IME,REQN,DEL 

70  FORMAT ( 'ENERGY  REQ,  INCIPIENT  MELT  ...  ',E10.4,'  JOULES/KG',/, 
S'BNERGY  REQ,  COMPLETE  MELT  ....  ',E10.4,'  JOULES/KG',/, 

$ 'EXCESS  ENERGY  AVAILABLE  .  ',E10.4,'  JOULES/KG') 

PV-0.0 
PL-100.0 
PS-0.0 
GOTO  100 
BNDIF 

...  IF  HASTE  HEAT  EXCEEDS  THE  ENERGY  REQ'D  TO  START  VAPORIZATION, 
...  BUT  IS  LESS  THAN  THAT  REQ'D  T'  COMPLETE  VAPORIZATION,  RESET  THE 
...  VALUE  OF  THE  ENERGY  AVAILABLE  FROM  THE  HASTE  HEAT  VALUE  TO  THE 
...  VALUE  REQ'D  TO  START  VAPORIZATION.  THIS  IMPLIES  THAT  SOME 
...  ENERGY  IS  AVAILABLE  FOR  VAPORIZING  A  PORTION  OF  THE  MATERIAL. 
...  NOTE:  THE  TEMPERATURE  RISE  EQUALS  THE  VAPORIZATION  TEMPERATURE 
...  OF  THE  MATERIAL. 
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zr  (BXH.OS.IVB.Ma).IXH.I.T.ZVS-t>HV)  THBN 

OIL-BXH-IVB 

IISgV«ZVI<t>HV 

TR-TV 

1IRZTS(2,80)  IVB,RKgiV,DIL 

80  rmaukTCiMBiioy  riq,  zncipiimt  vap -  ',8io.4,'  jodlbs/xg',/, 

S'KIBlIOy  RIQ,  OOHPLSn  VAP .  JOTLBS/KG' , 

S’milnoy  availabli  por  vap  .  *,bio.4,‘  jooues/xg*) 

PV>100.0*1»L/BV 
PL-100. 0-PV 
P8-0.0 
GOTO  100 
BMDZP 

c 

IP  (BXH.gr. IVR-t-HV)  THBN 

BCVAP-ZVB'I'BV 

PV-100.0 

PL-0.0 

P8-0.0 

NRITB  (2,90)  BCVAP 

90  FORMAT ( 'BMBROy  RBQ,  GOMPLBTB  VAP  .  ',B10.4, '  JOULBS/KG * , / , 

$•***  tHB  MATBRIAL  is  COMPLBTBLy  VAPORIZBD  ***•) 

GOTO  120 
BNDIF 
C 

100  WRITB(2,110)  TR,P8,PL,PV 

no  FORMAT CRBSZDOAL  MATBRIAL  TBMP .  *,F10.3,’  DB6-C PBRCBM 

$T  OF  SHOCKBD  AMD  RBLBA8B0  MATBRIAL  ,3X, • IH  SOLID  STATB  ... 

$, F6. 2, '%*,/, 3X, *IM  MOLTBH  FORM  ...  ’ ,F6.2, •%• ,/,3X, ’ IN  VAPOR  FORM 
$....  •,F6.2,*%’) 

C 

120  RBTURN 
BND 
C 

8UBROUTINB  TMCALC(V,PS,PL,PV,RP,RT,TR,TM,TV,TS,LP,DP,BHN,TSOL, 

$  MS,ML,MV,MTARG,HC»PT,Sy,SU,UTS,DPS,1ITP,UPP, 

$  MUSM,MSR) 

IMPLICIT  DOUBLB  PRBCISION  (A-H,0-Z) 

DOUBLB  PRBCISION  LP,MTARO,MS,ML,MV,MUSN,MSR,MTSR 
INTBGBR  HCOPT 
C 


C .  THIS  SOBROUTINB  CALCOLATBS  THB  MASSBS  OF  SOLID,  LIQUID,  AMD 

C .  GASBOUS  TAR6BT  MATBRIAL  IN  THB  DBBRIS  CLOUD. 

C 


LP-LP*2.54 

DP-DP*2.54 

TS-TS*2.S4 

T22-0.72*DP 

CST-C8(UT8,UTP) 

CSP-CS(UPS,UPP) 

T1N-CSP4-UPS-UPP 

TID-CST-UTS+UTP 

Tl-TIH/TID 

T2-CST/CSP 

T3-UTS/UPS 

T42-LP*T1*T2*T3 

PZ-4.0*ATAM(1.0) 

D8- ( 1 . 5*DP*DP*LP ) **0 . 3333333333 


C 

C .  MOTBt  THB  PROJBCTILB  LBMGTH  AND  DBNSITT  PASSED  TO  THE  ROLE 

C .  DIAMETER  CALCULATOR  SUBROUTINES  IS  THE  TOTAL  LENGTH  OF 

C .  THB  PROJECTILE  AND  ITS  AVERAGE  MATERIAL  DENSITT, 

C .  RESPBCTIVBLy 


u  u  u  u  u 


ir  (Hoc^.BQ.i)  cao:j.  dholii(ts,dp,rp,v,oh) 

IP  (HOOPT.1Q.2)  CKLL  DHOLB2(T8,DP,RP,RT,I.P, V,BHN,8Y,I>S,DH)  , 

IP  (BOOPT.IQ.3)  C8IJ.  DHOLI3(T8,DP,RP,RT,LP,V,BaN,DH) 

IP  (HOOPT.IQ.4)  GALL  OHOLB4(T8,OP,RP,RT,V,8y,OH) 

IP  (BOOPT.IQ.5)  DH«DP 
NTARSaPI * (DH/2 . 0 ) * ( DH/2 . 0 ) *T8*RT 
HT8ROPI* ( DP/2 . 0 )*( DP/2 . 0 ) *TS*RT 
C 

CALL  80RT(T8,T22.T42,T0) 

P8R«TO/TS 

•  C 

IP  (P8R.lt. 1.0)  THBN 
NRITB  (2,50)  TO,TS 

50  P(»HRT(/, 'DEPTH  OP  TARG  MATL  SUBJ  TO  S6R  -  ',P9.4, '  CM  <  TARO  THIC 
$101188  >  *,P9.4,'  CM. ',/,' PROGRAM  HALTED  IN  8UBROOTINE  TMCALC. ' ) 
NRITE  (*,50)  TO,TS 
8TOP 
EHDIP 
C 

MSR>P8R*MTSR 

iai8M-MT8R-M8R 

M8«(P8/100.0)*MSR 

ML-(PL/100.0)*MSR 

MV-(PV/100.0)*MSR 

T80L-MTAR0-M8R4-MS 

TN8-MTARG-T80L 

C 

WRITE  (2,20)  RP,LP,OP,T8,DH,MTARG,TO,MSR,MDSM,MSR,MS,ML,MV, 

$  T80L,TN8 

{?  20  PORMAT(/, 'AVO  PROJ  DEM8ITT  .  ',P9.4, ’  GM/CD. CM. ’,/, 'TOTAL  PR 

$OJ  LENGTH  .  •,P9.4,*  CM* ,/, 'PROJECTILE  DIAMETER  ....  *,P9.4, 

$•  CM*,/, 'TARO  PLATE  THICKNESS  ...  *,P9.4,'  CM',/,'TARG  PLATE  ROLE 

$DIA  ....  ',P9.4,'  CM',//, 'MASS  OP  REMOVED  TARG  MATL .  ', 

^  $P9.4,'  GMS',/, 'DEPTH  OP  TARG  MATL  SUBJ  TO  SfiR  .  ’,P9.4,'  CM', 

$/,'TOT  MASS  OP  TARG  MATL  SUBJ  TO  S6R  ...  ',P9.4,'  GMS',/, 'MASS  OP 

$UMSH  TARGET  MATL .  ',F9.4,'  GMS',/, 'MASS  OF  SH  AND  REL 

$TARO  MATL  .  ',F9.4,'  GMS* ,/, 3X, 'MASS  OP  S6R  SOLID  MATL  .... 

$ .  ',F9.4,  '  GMS ',/,3X, 'MASS  OF  S6R  LIQUID  MATL .  ', 

$F9.4,'  GMS ',/,3X, 'MASS  OF  S6R  VAPOR  MATL  .  ',F9.4,'  GMS' 

$,/, 'TOTAL  SOLID  MASS  COMPONENT  .  ',F9.4,'  GMS ' , / , ' TOT AL  N 

SON-SOLID  COMPONENT  .  ',F9.4,'  GMS') 

C 

RETURN 

END 

C 

SUBROUTINE  DHOLEl (TS,DP,RP, V,DH) 

IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z) 

THIS  SUBROUTINE  CALCULATES  THE  HOLE  IN  A  THIN  PLATE  DUE  TO 
THE  NORMAL  IMPACT  OF  A  SOLID  RIGHT  CIRCULAR  CYLINDER  USING 
KAPPII  EQUATION  HSSOl  (KAPPII  USER'S  MANUAL). 

A« 

B- 
C- 

T1«DEXP(A*RP) 

T2-1 . 0-DEXP ( -C*TS/DP ) 

DHDP»T1* ( 1 . 0+B*V*T2 ) 

DH-DP*DHDP 
IF  (DH.LT.DP)  DH-DP 
RETURN 
END 

SUBROUTINE  DHOLE2(TS,OP,RP,RT,LP,V,BHN,Sy,DS,DH) 
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nnonnn  nonno 


IMPLICIT  OOOBLI  PRKCISION  (A>H,0->Z) 

DOQBLI  PMBCI8I0H  1CT,LP 

THIS  SUBROOTIMI  CALC0IATI8  TOX  HOLI  IM  A  THIN  PLATE  DUX  TO 
THX  NORNAL  IMPACT  OP  A  SOLID  RIGHT  CIRCULAR  CYLINDER  USING 
XAPPII  EQUATION  H8S02  (KAPII  USER'S  MANUAL). 


A- 

B> 

EL 

XT« 

PH- 

T1*RT/RP 

T2-RP*V*V/ ( 2 . 0*X*BHN ) 

R1«A* (T1**B) * (T2**0. 3333333333 ) 

U>D8QRT(T1) 

pf|R«V*V*  (RP/2 . 0/SY)  *U/  ( 1 . 0+U) 

SV>1 . 0>0XXP  ( -1000 . 0*PtfR ) 

AR-LP/DP 

IF  (AR.GX.1.0)  PDP«(8V/U)*(AR-1.0+Rl/2.0) 

IF  (AR.LX.1.0)  PDP-(SV/U)*(R1*DP/DS/2.0)*(AR**0. 3333333333) 
p«PDP*DP 

DRC-0 .533*  (RP/RT)  <i>0 . 467 
RC-P/DSQRT(DRC) 

Q-1.0-<OP/2.0/RC)**2 

IP  (Q.LT.0.0)  THEN 

DH— 1.0 

RETURN 

ENDIP 

IP  (Q.GE.0.0)  THEN  '' 

T-KT*TS+P* ( 1 . 0-DSQRT (Q) ) 

DH«2 .0*PN* (RC/P) *DSQRT(T* (2 .0*P-T) ) 

IP  (DH.LT.DP)  DH-DP 

RETURN  m 

ENDIP 

END 

SUBROUTINE  DHOLB3(TS,DP,RP,RT,LP,V,BHN,DH) 

IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z) 

DOUBLE  PRECISION  K,LP 

THIS  SUBROUTINE  CALCULATES  THE  HOLE  IM  A  THIN  PLATE  DUE  TO 
THE  NORMAL  IMPACT  OP  A  SOLID  RIGHT  CIRCULAR  CYLINDER  USING 
KAPPII  EQUATION  HSAOl  (KAPPII  USER’S  MANUAL). 


A> 

B- 

C- 

B» 

P» 

0« 

K«P*((BHN/RP)*’*C) 

Rl-RP/RT 

R2«(3.0*LP)/(2.0*DP) 

R3«(RP*V*V) / (2 . 0*B*BHN) 

DR«A* (R1**B) * (R2**C) * (R3**0. 3333333333 ) 

R4«(TS/DP) **0.6666666666 

DHDP-1 . 0+ ( DR- 1 . 0 ) * ( 1 . 0-DEXP ( -K*R4 ) ) 

DH-DP*DHDP 

IP  (DH.LT.DP)  DH>DP 

RETURN 

END 

SUBROUTINE  DHOLE4(TS,DP,RP,RT,V,SY,DH) 
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IMPLICIT  DOUBLE  PRBCISKHf  (A-R,0>S) 


C 

C .  THIS  SUBROUTINB  CALC0LATB8  THB  HOLE  IN  A  THIN  PLAT!  DUB  TO 

C .  THB  NORMAL  IMPACT  OP  A  SOLID  SPHBRB  USING  THB  PBN4.V10  HOLE 

C .  DIAMETER  EQUATION  (BOBING-DI80>30550-2). 

C 


R1«DP/TS 

R2-1000 . o*RP*v*v/sy 
R3-RP/RT 

Tl»Rl*(R2**0.415)/(R3**0.15)/29.9 
DHT8«11 . 02* ( 1 . 0-DEXP ( -T1 ) ) 

DB-TS*DRTS 
IF  (DH.LT.DP)  DH-DP 
RETURN 
END 
C 

SUBROUTINE  UPPCAL(RPA,RPB,KPA,XPB,OOPA,COPB,V1,UPP) 
IMPLICIT  DOUBLE  PRECISION  (A>B,0-8) 

DOUBLE  PRECISICNI  KPA,KPB 


C 

C .  THIS  SUBROUTINB  CALCULATES  THB  PARTICLE  VELOCITY  IN  A  MATERIAL 

C .  *B'  DUE  TO  A  SHOCK  NAVE  THAT  HAS  ENTERED  MATERIAL  'B*  FROM  AM 

C .  ADJACENT  MATERIAL  'A* 

C 


A«RPA*KPA-RPB*KPB 
B-RPA*COPA<fRPB*COPB-l-4 . 0*RPA*KPA*V1 
Pl«2 . 0*RPA*V1*  (COPA-t-2 . 0*KPA*V1 ) 

DISC«B*B-4 . 0*A*P1 
UPP>(B-SQRT<DISC) )/(2.0*A) 

r  c 

RETURN 

END 

C 

SUBROUTINB  SORT(A,B,C,SN) 

IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z) 

SM-A 

IF  (SM.GT.B)  SM>B 
IF  (SM.GT.C)  SM>C 
RETURN 
END 
C 

SUBROUTINE  DCVEL (UFSTl , UFSPl , V, MTARG,MPROJ, MPLYR, MTSR, MPSR, MUSTM 
$  ,MUSPN,PNSSNR,BXT,EXP,FSRP,NPNAT) 

IMPLICIT  DOUBLE  PRECISKM  (A-H,0-Z) 

DOUBLE  PRECISION  MDC,NPROJ,NTARG,MUSTN,MTSR 

DOUBLE  PRECISKHI  MUSPN(IO) ,MPSR(10) ,EXP(10) ,FSRP( 10) ,MPLyR(10) , 

$  UFSP1(10),PNSSNR(10) 

CHARACTER*!  ALLSR 
C 

C. ...  THIS  SUBROUTINE  ONIPUTES  DEBRIS  CLOUD  VELOCITIES  AND  THE  DEBRIS 
C....  CLOUD  SEMI-CONE  ANGLE.  FIRST,  A  CHECK  IS  PERFORMED  TO  SEE  IF 
C. ...  AMY  UNSHOCKED  PROJECTILE  MATERIAL  REMAINS. 

C 

ALLSR>'Y' 

*  DO  10  1-1,NPMAT 

IF  (FSRP(I).NS.l)  THEN 
ALLSR- 'N* 

GOTO  100 

«  EMDIF 

10  CXMITINUE 
C 

C. ...  IF  NO  UNSHOCKED  PROJECTILE  MATERIAL  REMAINS,  THEN  ALL  OF  THE 
C....  PROJECTILE  MATERIAL  COMBINES  WITH  THE  TARGET  MATERIAL  TO  FORM 
C....  THE  DEBRIS  CLOUD. 


131 


c 

VL^OrSTl 

vKPV-urspi  (MmiiT) 
vrnii»iiTiinT*Yf 

■XMP-CIIPIIOJ/IOOO.O)  •  (V*1000.0)*  (V*1000.0) /2 .0 

ll<«XT*lCT8m/ 1000 . 0 

B2«0.0 

DO  20  i-ifUnar 
B2>B2-t>BXP  ( Z  )  *MP8R  ( I )  / 1000 . 0 
20  wQTTW^B 

B3«  (  MntOJ/ 1000 .  Ot-HTMU/ 1000 . 0 )*(  VCXM* 1000 . 0 )*(  VOOM* 1000 . 0 )/ 2 . 0 
DBI.B-BI1IP>B1>X2-B3 
DBUI«(IIPllOJ-«-lCTARO)  /lOOO.O 
VBXP-8QRT<2  .O^DBU/DBLM) /lOOO.O 
VBXPl-(VL-VR)/2.0 
C 

1IRZTB(2,30) 

30  PantAT(/,'IIO  0M8H0CKB0  RBSIDOBL  PROJBCTIU  PRAOMBIIT  RXMAZMS .  ' ,  / , 
S'Mnr  SOLID  PPOJBCTILB  AMD  TAMOBT  MATBRIAL  RBMAINING * , / , *  IS  LIKBLY 
$tO  BB  PRAQMBMTBD. • ) 

NRITB(2,40)  VR,VR/V,VOOM,VOOH/V,VL,VL/V,VBXP,VBXP/V,VBXP1, 

$  VBXPl/V 

40  PORMAT(/,  DBBRIS  CLOOD  VBLOCITY  SOMNARY  . . . * ,/,3X, 'RBAR  SURPACB  VE 

SLOCITY  (VR)  . 'fWl.S,'  KM/SBC  (>' ,F5.3, 'V)  •  ,/,3X, 'CBMTBR 

$-OP-MASS  VBLOCITY  <VOON)  . ‘,¥1.3,*  XM/SBC  (»•  ,F5.3, ’V) ' 

$3X, ’LBADIMO  BOOB  VBLOCITY  (VL)  . ',F7.3,'  KM/SBC  (•', 

$FS.3,'V)’,/,3X, *BXPAMSI(Xf  VBL  (VBXP  ->  BMBRGY  CCMIS)  ...',F7.3,'  KM 

$/SBC  (-‘,FS.3,'V)',/,3X, 'BXPAMSKMI  VBL  (VBXP  -  (VF-VR)/2)  . *, 

$F7.3, •  XM/SBC  (-‘/FS,!, -V)*) 

C 

OOIO  200 
C 

C. . . .  IF  SOMB  ONSBOCKBD  PROJBCTILB  MATERIAL  REMAINS,  THEM  INVCHCB  THE 
C....  ALTERNATIVE  DBBRIS  CLOOD  VELOCITY  CHARACTERISATION  SCHEME. 

C 

100  CONTINUE 
VL-UFSTl 
TPMH-0.0 
TMUSPN-0.0 
MDC-NPROJ•^MTARG 
DO  no  I«1,NPNAT 
TPMH-TPHR't'EXP  ( I )  *NPSR  ( I ) 

TNDSPM-TMUSPN+MUSPN  ( I )  -t-PNSSNR  ( I ) 

MDC-MDC-MUSPN( I ) -PNSSNR( I ) 
no  COMTIMOE 

A-2  .  Ot-MDC/TMOSPN 
115  B«VL+(NPROJ/TNUSPN)*V 
C1-VL*VL 

C2-  (MPROJ/TMDSm-l .  0 )  •  (MPROJ/NDC)  *V*V 

C3-2 . 0*  ( EXT*NTSR/ 1000 . 0<^TPNB/ 1000 . 0 )  / 1000 . 0/MDC 

C-C14C2<K:3 

DISC- (B/A)*(B/A) -C/A 

IF  (DZSC.LT.O.O)THEN 

VL-VL-O.l 

OOIO  115 

BMDIF 

VCOM-B/A-SQRT ( DISC ) 

VEXP-VL-VCOM 

VRES- (MPROJ*V-MDC*VCOM) /TNOSPM 
C 

NRITB(2,120) 

120  F(»OIAT(/,*SOMB  UMSHOCKBD  PROJ  MATL  REMAINS  ...') 

1IRITE(2,130)  TMDSPM,TMOSPM/MPROJ,VRES,VRES/V 
130  F0RMAT(A3X, *TOT  MPROJ,UNSH  ...',F7.3,'  GMS  («* ,F5.3, 'MPROJ) ' ,/, 
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$3X, ’AVG  VP,UM8H . •,F7.3,'  KN/S  (-’ ,F5.3, 'V) ' 'ALL  OTHKR  SO 

$LZO  FROJICTZLI  MATERIAL  (IF  AMY)  IS  LIKSLY  TO  SB  FRAONBHTBD. ' 

$'A1IY  SOLID  TAROBT  MATERIAL  REMAIMIHG  IS  ALSO  LIKELY  TO  BE  FRASMEMT 
fED. 'tA'TRE  DEBRIS  CLOUD  CONSISTS  OF  SMOCKED  AMD  RELEASED  FROJECTI 
$LE  MATERIAL  AMD' 'ALL  EJECTED  TARGET  MATERIAL. ' ) 

C 

1IRITE(2,140)  VOOM,VOOM/V,VL,VL/V,VEXP,VEXP/V 

140  FORMATC A 'DEBRIS  CLOUD  VELOCITY  SUMMARY  3X, 'DEB  CLD  CEMTER- 

$OF-MA88  VBL  (VOOM)  ...',F7.3,'  KM/SEC  (-' ,FS.3, 'V) ' , /,3X, 'DEB  CLD 

SLEADIMG  EDGE  VEL  (VL)  . ',F7.3,'  KM/SEC  (-' ,F5.3, 'V) ' ,/,3X, 

$*DEB  CLD  EXPAM8IOM  VEL  (VEXP)  . ',F7.3,'  KM/SEC  (»',FS.3, 

$'V)') 

C 

200  DCAMO»ATAR(VEXP/VOOM) 

DCANO- ( 180 . 0/3 . 141592 )  *DCAM6 
C 

1IRZTE(2,498)  OCAMG 

498  F(»ttlAT(//, 'DEBRIS  CLOUD  HALF-ANGLE . ',F7.3,'  DBG') 

C 

RETURM 

EMD 

C 

DOUBLE  PRECISIOM  FOMCTICMI  CS(US,DP) 


IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z) 

C 

C .  THIS  FUNCTKMI  CALCULATES  THE  SPEED  OF  A  RAREFACTION  flAVE  IN  A 

C .  SHOCKED  MEDIUM 

C 

T1-(US-0P)/0S 

f  C8Q«US*US*(0.49-»-Tl*Tl) 

CS«DSQRT(CSQ) 


RETURN 

END 

C 

SUBROUTINE  RSINT(NPMAT,UTS,UTP,UPSA,UPPA,LPA,DP,TS, V,LO) 
IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z) 

DOUBLE  PRECISION  LO,L11,L41,LPTOT,UPPA(10) ,UPSA(10) ,CSP(10) , 
$VE(10),VD(10),VC(10),TB(10),XB(10),TD(10),TDSUM(10),LPA(10) 

C 

C .  THIS  SUBROUTINE  CALCULATES  THE  LOCATKHI  WITHIN  THE  PROJECTILE 

C .  AT  WHICH  THE  TARGET  RARBFACTKXI  WAVE  OVERTAKES  THE  PROJECTILE 

C .  SHOCK  WAVE 

C 

L11-0.72*DP*2.54 

CST»CS(UTS,UTP) 

VB-CST-UTP 

VA-UTS 

TC»(TS/VA) * ( VB4VA) / (VB+UTP) 

LPTOT-LPA ( 1 ) *2 . 54 
VC(1)-UTP 
DO  10  1«2,NPMAT 
VC(I)«UPPA(I) 

LPTOT»LPTOT+LPA( I ) *2 . 54 
10  OCXITINUB 
C 

I-O 

99  CONTINUE 
l»I>fl 

CSP(I)«CS(UPSA(I) ,UPPA(I) ) 

VE(Z)>UPSA(I)-V 

C 

C .  CCmSIDER  THE  FIRST  MATERIAL  LAYER  OF  THE  PROJECTILE. 

C 

IF  (I.EQ.l)  THEN 


133 


ooouo  oouoooo  ouooo  ooo  oo 


VD(1)*«S»(1)>0T» 

nTi^  ( 1 )  *  ( VD  ( 1 ) -fVC  ( 1 ) )  *TC/ ( VD  ( 1  >  - VI  ( 1 ) ) 

nt-HY 

n*<*iiy/vi(i) 

tWTIT 

IM«XR*V/VB(1) 

L41*(m4>XDT)*2.54 

CnLl«l.rA(l)*2.S4 

CaiLL  MKT(CaCLl,Lll,L41«LO) 

...  zr  OMLT  on  nK»BCTix.i  nmtiAL  is  imvolvid,  nroitii  mxm  vium 

...  M  Tn  ZMSTIOM  ALOM  TU  ntOJBCXXLS  AXIS  OF  THI  mmSICTIOII 

...  or  Tn  RAnrACTZoM  navi  aid  Tn  saocx  mavi. 
zr  (imiAT.IQ.l)  OOTO  40 

...  zr  Tnn  is  nan  than  on  ntoaicTiLB  MATnzAL  uum  and  Tn  min 
...  VAzn  ncms  Tn  TBicxnss  or  Tn  MATniAL  LAm,  Tns  Tn 
...  SHOCK  HAVE  ZS  OVnTAXXM  AT  A  POZMT  BIYOHD  Tn  THICXnSS  OF  Tn 
...  rZMT  LAXn.  ZMZTIALXSS  AARAY  nTRZKS  RIQDXnD  FOR  MMXT  LAYSR 
...  CAXXXIZATZOMS  AIB>  NOVI  OM  TO  THR  SKOOMD  MATIRIAL  LAYKK. 

zr  (lirMAT.aK.2.AlfD.IX).OT.Z.PA(l))  Tin 
Tl  ( 1 ) -X.rA(  1 )  /  (  VI  ( 1 ) -i-V) 
n(i)«vi(i)*Ti(i) 

TD(l)>(n<l)‘fVC(2}*TI(l))/(VD(l)4^VC(2)) 

OOir-(VC(l)-VC(2) )/(VD(l)4VC(2) ) 

TC  ( 1 )  «TD  ( 1 ) -toon  *TC 
OOTO  99 

nozr 

...  zr  Tnn  zs  non  than  on  frojictzli  mathial  layir  aid  Tn  mzm 

. . .  VALOR  LZI8  WZTBXR  Tn  FIRST  LAYBR,  THn  Tn  SHOCK  NAVI  IS  OVIR- 
. . .  TAXn  lY  A  RARSFACTION  NAVI  NITHIM  Tn  FIRST  LAYSR. 

zr  (liniAT.OI.2.AllD.LO.LI.LrA(l))  OOTO  40 

nozr 

. . .  OOMSIDIR  Tn  RIMAINING  PROJICTILI  LAYERS 

zr  (1. 01. 2)  THn 
VD(Z)>csr(X)-urPA(Z) 

XIS«0.0 
TSS«0.0 
Toirr»Tc 
OO  20  J-1,Z-1 
XIS-XIS-<'XI(J) 

TS8-TIS4'TI(J) 

TDirr-TDirr+TD ( j) -Tl ( j) 

20  OOHTZHOI 

XIT^  ( I )  •  ( VO  ( X  )  ♦VC  ( I ) )  •TO  IFF /  ( VD  ( I )  - VI  ( I ) ) 

XI8>XI8^XIT 

TIT»XIT/VI(I) 

TSS-ns^TIT 

XOTr-XIT*V/VI(I) 

XDT-XDTP*TIS/TIT 
L41-(XI8^XDT)*2.54 
CBXLZaO.O 
OO  30  J-1,I 

CHXLI*cnLI^LPA(  J)  *2 . 54 
30  OONTINOI 

CALL  8(»T(CHKLZ,Lll,L41,LO) 

...  IF  Tn  HIM  VALDK  IXCIIOS  Tn  CUNDLATIVI  LnOTH  OF  Tn  PROJICTILI 
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.  MATERIAL  THROUGH  THE  CURRENT  LAYER,  THEN  THE  SHOOC  WAVE  IS  OVER- 

.  TAKEN  AT  A  POINT  BEYOND  THE  THICKNESS  OF  THE  CURRENT  LAYER  AND  LO 

.  IS  SET  EQUAL  TO  THE  CURRENT  CUMULATIVE  LENGTH.  INITIALISE  ARRAY 

.  ENTRIES  REQUIRED  F<»  NEXT  LAYER  CALCULATICEIS  AND  HOVE  ON  TO  nlE 

.  NEXT  MATERIAL  LAYER  (UNLESS  THIS  IS  THE  PINAL  MATERIAL  LAYER  IN 

.  MBICH  CASE  RETURN  THE  TOTAL  PROJECTILE  LENGTH  AS  LO) . 

IP  (LO.BQ.CHXLI)  THEN 
IP  (I.EQ.NPNAT)  THEN 
LO-CHKLI 
GOTO  40 
BMDIP 

IP  (I.LT.NPNAT)  THEM 
TE(I)-LPA(I)/(VE(I)-»^V) 

XE(I)>VE(I)*TB(I) 

TD(I)«(XE(I)-fVC(I-t-l)*TB(I))/(VD(I)-*‘VC(I-»-l)) 
TD(I)-TD(I)+(VC(I)-VC(I+1))*TDIPP/(VD(I)+VC(I+1)) 

GOTO  99 
EMDIP 
ENDIP 

.  IP  THE  MIN  VALUE  LIES  NITHIN  THE  CURRENT  LAYER,  THEN  THE  SHOCK 

.  NAVE  IS  OVERTAKEN  BY  A  RAREPACTIQN  NAVE  NITHIN  THE  CURRENT  LAYER. 

IP  (LO.LT.CHKLI)  GOTO  40 
EMDIP 

40  RETURN 
END 

SUBROUTINE  PMCALC(I,UPS,UTS,UPP,UTP,RP,PS,PL,PV,TS,LP,DP,TSOL, 
$NS,ML,NV, NPLYR, MUSM, MSR, PSR, PGMSTP, ISTOP, LO, CHXL ) 

IMPLICIT  DOUBLE  PRXCISICW  (A-H,0<-S) 

DOUBLE  PRECISION  LP,LO,LSR,MPLYR,MS,NL,MV,MSR,MUSM 
CHARACTER*!  PGMSTP 

.  THIS  SUBROUTINE  CALCULATES  THE  MASSES  OP  SOLID,  LIQUID,  AND 

.  GASEOUS  PROJECTILE  MATERIAL  IN  THE  DEBRIS  CLOUD. 

PI«4.0*ATAN(I.0) 

LP-LP*2.54 

MPLYR-PI* (DP/2 . 0) * (DP/2 . 0) *LP*RP 

IP  (CHKL.GT.LO)  THEN 
LSR«LP> ( CHKL-LO ) 

PGMSTP-'Y’ 

EMDIP 

IP  (CHKL.LE.LO)  THEN 
L8R-LP 
PGNSTP«'N' 

EMDIP 
PSP-LSR/LP 
MSR>>PSR*NPLYR 
MUSM-NPLYR-MSR 
MS« ( PS/100 . 0 ) *MSR 
ML«(PL/100.0) *MSR 
MV-(PV/100.0)*MSR 
TSOL-NUSN-fNS 
TNS-MPLYR-TSOL 
C 

NRITB  (2,100)  RP,LP, OP, NPLYR 

100  PORMAT(/, 'PROJECTILE  LAYER  DENSITY  .  ',P9.4,'  GM/CU.CM. ' , /, 

S'PROJECTILE  LAYER  LENGTH  .  ',P9.4,'  CM ',/,' PROJECTILE  LAYER  D 

$IAMETER  ....  ',P9.4,'  CM' ,/, 'PROJECTILE  LAYER  MASS  .  ',P9.4 


man  (2,200)  I,O,X.m,M08M«M8K,ll8,ia.,KV,T8OL,TM8 

200  ioMaT(/,'nioj  lihoth  iosrk  moivs  bits  smsvb  ....  ',r9.4, *  cm',/, 

I'LlMOn  09  ntOJ  LTK  NSTL  80BJ  TO  SSR  ...  ',89.4, '  CM',/,  *MSS8  OF  U 

$M8B  PBOJ  LTK  NSTBRIAL .  ',89.4,'  OHS', /, ‘MASS  OF  8H  AMD  RBL 

$  PBOJ  LTM  NATL .  •,F9.4,'  QMS' ,/,3X, 'MASS  OF  SU  801.ZD  LYR  MA 

ITZ. .  ',F9.4,  '  0M8',/,3X,  *NA88  OF  8fiR  LZQDXD  LYR  NATL . 

$...  ',F9.4,'  QMS *,/,3X, 'MASS  OF  8SR  VAP<»  LYR  MATL .  ', 

$F9.4, '  0N8',/, 'TOTAL  SOLID  LAYBR  MASS  OOMPONXMT  .  ',F9.4, 

$*  ON8',/, 'TOTAL  MOM-SOLID  LAYBR  OOMPOMXRT  .  ',F9.4,'  QMS') 


Z8TOF-I 

RBTORN 

BMD 


# 


t 


APPENDIX  B  ~  Input  Fde  INDATA 
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- CO - 

- NO - 

— T.VAP — 
-OLISTR-- 

— 'ALFA - 

— H.FUS-- 
~ROPT - 

•>-OAMO - 

— -CPL - 

—ALFA - 

•“•“•mix 

— TLDSTR-- 

- CPS - 

—H.VAP — 

- BPS - 

— BBTA - 

AL 

AI41NIiniM 

5.380 

1.340 

2.712 

2.130 

120.0 

0.103B'K>8 

0.35 

0.240B-04 

0.235 

0.255 

0.005 

660.0 

2450.0 

95.0 

2450.0 

8.0 

5.0 

290.0 

434.02 

A1 


2XXX  ALUM 
0.106B<^08 
640.0 
314.0 

5.350 
0.33  0. 
2450.0 
433.02 

1.340 

.209B-04 

85.0 

2.800 

0.212 

2450.0 

2.000 

0.242 

5.0 

120.0 

0.005 

5.0 

A2 

5XXX  ALUM 

5.310 

1.340 

2.670 

2.000 

84.0 

O.lOlB+08 

0.33  0. 

.225B-04 

0.215 

0.245 

0.005 

641.0 

2450.0 

85.0 

2450.0 

5.0 

5.0 

211.0 

312.02 

A3 

6XXZ  ALtm 

5.380 

1.340 

2.700 

2.000 

93.0 

0.100B4-08 

0.33  0. 

.233B-04 

0.212 

0.242 

0.005 

652.0 

2450.0 

85.0 

2450.0 

5.0 

5.0 

202.0 

268.02 

A4 

7XXX  ALOM 

5.290 

1.340 

2.810 

2.000 

150.0 

0.103B'«>08 

0.33  0. 

.221B-04 

0.217 

0.245 

0.005 

636.0 

2450.0 

85.0 

2450.0 

5.0 

5.0 

479.0 

531.02 

BB 

BBRYLLIUM 

7.975 

1.124 

1.820 

1.160 

120.0 

0.419B-^08 

0.08  0. 

,1408-04 

0.570 

0.832 

0.005 

1281.0 

2884.0 

260.0 

8195.0 

5.0 

5.0 

225.0 

300.02 

CD 

CADMIUM 

2.307 

1.640 

8.640 

2.270 

24.0 

0.6728+07 

0.33  0. 

.343B-04 

0.058 

0.063 

0.005 

321.0 

765.0 

13.5 

212.0 

5.0 

5.0 

34.0 

52.02 

CU 


COPPBR 

0.190B+08 

1083.0 

240.0 

3.940 
0.34  0. 
2590.0 
340.02 

1.489 

,170B-04 

49.0 

8.930 

0.097 

1150.0 

2.000 

0.114 

5.0 

37.0 

0.005 

5.0 

BP 

BPOXY 

3.020 

1.520 

1.180 

0.800 

-1.0 

0.650B+06 

0.50  0. 

.  500E-04 

0.250 

0.285 

-1.0 

350.0 

-1.0 

-1.0 

-1.0 

-1.0 

-1.0 

-1.0 

-1.02 

FB 

ZRM 

0.290B+08 

1539.0 


4.580  1.490 

0.30  0.120B-04 
3035.0  65.0 


7.870 

0.120 

1591.0 


1.570 

0.150 

5.0 


95.0 

0.005 

5.0 
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T.»>n 

0.2001^07 

327.0 

9.0 

2.030  1.470 

0.45  0.293B-04 
1740.0  6.0 

17.02 

11.340 

0.031 

210.0 

2.770 

0.036 

10.0 

7.0 

0.005 

10.0 

IX 

UtXAM 

2.750 

1.480 

1.180 

0.860 

37.0 

0.3451406 

0.50  0. 6501-04 

0.290 

0.315 

-1.0 

225.0 

-1.0 

-1.0 

-1.0 

-1.0 

-1.0 

1.0  -1.02 


NO 

NOLyBOBinm 

0.4601408 

2610.0 

350.0 

5.173 
0.31  0. 
5555.0 
450.02 

1.220 

.0611-04 

70.0 

10.200 

0.079 

1242.0 

1.520 

0.104 

5.0 

200.0 

0.005 

5.0 

NI 

mCKBL 

4.667 

1.530 

8.860 

1.800 

200.0 

0.3301408 

0.30  0. 

.  1431-04 

0.130 

0.157 

0.005 

1454.0 

2865.0 

74.0 

1523.0 

5.0 

5.0 

59.0 

317.02 

PT 

PZATimiM 

3.680 

1.500 

21.370 

2.940 

70.0 

0.2771408 

0.39  0. 

,1101-04 

0.037 

0.042 

0.005 

1769.0 

4349.0 

26.0 

632.0 

10.0 

10.0 

100.0 

200.02 

81 

304  8T11L 

4.590 

1.550 

7.910 

1.670 

237.0 

0.2841408 

0.28  0, 

.1121-04 

0.110 

0.125 

0.005 

1425.0 

3035.0 

65.0 

1590.0 

5.0 

5.0 

250.0 

500.02 

82 

430  STBIL 

4.680 

1.550 

7.830 

1.670 

251.0 

0.2991408 

0.29  0. 

.1041-04 

0.110 

0.125 

0.005 

1470.0 

3035.0 

65.0 

1590.0 

5.0 

5.0 

275.0 

480.02 

83 

4340  8T11L 

4.570 

1.550 

7.830 

1.670 

290.0 

0.2901408 

0.30  0. 

.1121-04 

0.110 

0.125 

0.005 

1510.0 

3070.0 

65.0 

1590.0 

5.0 

5.0 

469.0 

745.02 

TA 

TANTALUM 

3.374 

1.201 

16.650 

1.690 

200.0 

0.2601408 

0.35  0. 

.0651-04 

0.033 

0.039 

0.005 

2996.0 

5425.0 

38.0 

1007.0 

10.0 

10.0 

288.0 

380.02 

8N 

TIN 

2.560 

1.520 

7.280 

1.850 

4.0 

0.6031407 

0.33  0. 

.2691-04 

0.058 

0.062 

0.005 

235.0 

2450.0 

14.0 

580.0 

10.0 

10.0 

23.0 

31.02 

TI 

TITANIUM 

4.786 

1.049 

4.512 

1.100 

330.0 

0.180B-<-06 

1676.0 

810.0 

0.30  0. 
3260.0 
1013.02 

.1008-04 

99.0 

0.150 

2182.0 

0.167 

5.0 

0.005 

5.0 

V 

TONGSTBN 

4.150 

1.237 

19.170 

1.480 

400.0 

0.5908-^08 

0.30  0. 

.0408-04 

0.035 

0.046 

0.005 

3410.0 

5900.0 

53.0 

1054.0 

10.0 

10.0 

1379.0 

1517.02 

18 

tIMC 

3.042 

1.500 

7.140 

2.150 

82.0 

0.1088408 

0.33  0. 

.2748-04 

0.100 

0.115 

0.005 

420.0 

907.0 

25.0 

420. 0 

10.0 

10.0 

138.0 

183.02 

KO 

GOLD 

3.060 

1.570 

19.240 

3.100 

33.0 

0.1248-<-08 

0.42  0. 

.  1618-04 

0.034 

0.038 

0.005 

1063.0 

2960.0 

16.0 

413.0 

10.0 

10.0 

95.0 

125.02 

80 

S1LV8R 

3.230 

2.500 

10.490 

2.500 

25.0 

0.1208-1^08 

0.37  0. 

.2118-04 

0.062 

0.071 

0.005 

961.0 

2210.0 

25.0 

554.0 

10.0 

10.0 

55.0 

175.02 

MO 

M80M8SiaN 

4.490 

1.240 

1.740 

1.500 

45.0 

0.6408-»^07 

0.29  0. 

.3008-04 

0.295 

0.336 

0.005 

650.0 

1110.0 

88.0 

1326.0 

5.0 

5.0 

197.0 

278.02 

XX 
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APPENDIX  C  -  Output  FUe  IMPOUT 
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HYPIRVBLOCZTY  IMPACT  OP  A  38.240  CM  MULTI-MATERIAL  PROJECTILE  OH  A 
ALUMINUM  TARGET  AT  A  10.00  KM/SEC  IMPACT  VELOCITY 

TARGET  MATERIAL  PROPERTIES  ... 


MAT 

m 

ALUMINUM 

CO 

m 

5.380  KM/S 

K 

m 

1.340 

RBO 

m 

2.712  GM/CU.CM 

TS 

m 

.317  CM 

PROJECTILE  MATERIAL  PROPERTIES  (DP  -  2.540  CM)  ... 


NAT 

1 

« 

ALUMINUM 

CO 

m 

5.380  KN/S 

K 

m 

1.340 

RHO 

m 

2.712  GN/CU.CM. 

LP 

m 

.254  CM 

MAT 

2 

m 

4340  STEEL 

CO 

m 

4.570  KM/S 

K 

m 

1.550 

R») 

m 

7.830  GM/CU.CM. 

LP 

m 

.254  CM 

MAT 

3 

m 

TUNGSTEN 

CO 

m 

4.150  KM/S 

K 

m 

1.237 

RHO 

m 

19.170  GM/CU.CM. 

LP 

m 

.254  CM 

*****  target  material  release  CALCULATION 


> 


INITIAL  CONDITIONS  FOR  TARGET  MATERIAL 


PARTICLE  VELOCITY . UP  - 

SHOCK  NAVE  SPEED . US  - 

HUG<HIIOT  IMPACT  PRESSURE  ....  PH  » 

HUG<X(IOT  IMPACT  ENERGY . EH  - 

SPECIFIC  VOLUME  AT  REST . VO  - 


SPECIFIC  VOLUME  AT  IMPACT  ...  VI  - 


5.000  KM/S 
12.080  KM/S 

163.805  6PA  -  1.617  NBAR 

.1250E+08  JOULES/KG 
.369  CU.CM./GM 
.216  CU.CM./GM 


f 


PARAMETERS  REQUIRED  FOR  CALCULATING  TARGET  MATERIAL  RELEASE  FROM  SHOCKED 
STATE  USING  THE  TILLOTSON  EQUATION  OF  STATE: 


TARG  MATL  LIN.  COEF.  OF  THERM.  EXP. 

TARG  MATL  SP  HEAT  (SOLID)  . . 

TARG  MATL  SP  HEAT  (LIQUID)  . . 

TARG  MATL  AMB  N-GRUN  COEF  (CAL,INP) 
TARG  NATL  YIELD  STRENGTH . . 


TARG  MATL  BRN  HDNS  NO 


E 

.7102E-M1  N/SQ.M. 

NU 

.350 

K 

.7891E-M1  N/SQ.M. 

ALFA 

.2400E-04  /DEG-C 

CPS 

.235  CAL/GM/DEG 

CPL 

.255  CAL/GM/DEG 

GAMO 

2.129,  2.130 

SY 

290.000  MPA 

SU 

434.000  MPA 

BHN 

120.000 

TM 

660.00  DEG-C 

TV 

2450.00  DEG-C 

HF 

95.00  CAL/GM 

HV 

2450.00  CAL/GM 

IME 

.6492E-«-06  JOULES /KG 

IVE 

.2958E-»-07  JOULES/KG 

ADDITICniAL  PARAMETERS  REQUIRED  FOR  CALCULATING  MATERIAL  RELEASE  FROM 
SHOCKED  STATE  USING  THE  TILLOTSON  EQUATION  OF  STATE: 

AA  •  .7850E-i-ll  N/SQ.M. 

BB  -  .4831E-I-11  N/SQ.M. 
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h  -  .5000 

B  -  1.6292 

JUUr  -  5.0000 

BBT  •  5.0000 

BOX  ■  .6687B<«-07  JDOLBS/KO 

BON  >  1.0000 

BO  -  .6687B-I-07  JOOLBS/KO 

B8  «  .2958B-t-07  JOOUSS/KG 

HV  ■  .1026B-«-08  JOIILB8/KG 

BSP  >  .1321B4-08  JOOLBS/K6 

V8  •  .4170  CU.CM./GM 

BPS  •  .0050 


BMD-STATB  CALCOXATXON  RESULTS  USING  THB  TILLOTSMI  BOS  ... 


MBTBRXAL  PIN  SP  VOL  <VP)  .  .494  CU.CM./GM 

MRTBRXAL  SHOCK  BMBRGY . 1250B-«-08  JOULBS/KG 

NRTBRIAL  BMBRGY  RBCOVBRBO . 9992B-^07  JOULBS/KG 

WRSTB  BBRT  GBNBRRTBO . 2508B+07  JOULBS/KG 

BMBRGY  RBQ,  IMCIPXBNT  MELT  ...  .6492B-i-06  JOULBS/KG 

BMBRGY  RBQ,  COMPLETE  MELT . 1047B+07  JOULBS/KG 

BXCBSS  BMBRGY  AVAILABLE . 1461B-«-07  JOULBS/KG 

RBSIDUAL  MATERIAL  TEMP  .  2028.692  DBG-C 


PBRCBMT  OF  SHOCKED  AND  RELEASED  MATERIAL  ... 
IN  SOLID  STATE  ...  .00% 

IM  MOLTEN  FORM  ...  100.00% 

IN  VAPOR  FORM . 00% 


c 


FREE  SURF  VBL  (UP+UR)  .  10.652  KM/SBC 

FREE  SURF  VBL  (2.0*UP)  .  10.000  KM/SBC 


AVG  PROJ  DENSITY  . 

TOTAL  PROJ  LENGTH  . 

PROJECTILE  DIAMETER  .... 
TARO  PLATE  THICKNESS  ... 
TAR6  PLATE  HOLE  DIA  .... 


9.9040  GM/CU.CM. 

.7620  CM 
2.5400  CM 
.3175  CM 
5.3722  CM 


MASS  OF  REMOVED  TARG  MATL .  19.5174  GMS 

DEPTH  OF  TARG  MATL  SUBJ  TO  S6R .  .3175  CM 

TOT  MASS  OF  TARO  MATL  SUBJ  TO  S6R  ...  4.3631  GMS 

MASS  OF  UNSH  TARGET  MATL .  .0000  GMS 

MASS  OF  SH  AND  REL  TARG  MATL .  4.3631  GMS 

MASS  OF  S6R  SOLID  MATL .  .0000  GMS 

MASS  OF  S6R  LIQUID  MATL .  4.3631  GMS 

MASS  OF  S6R  VAPOR  MATL .  .0000  GMS 

TOTAL  SOLID  MASS  COMPCHIENT .  15.1543  GMS 

TOTAL  NON-SOLID  COMPCHIENT .  4.3631  GMS 


PROJECTILE  MATERIAL  RELEASE  CALCULATIONS,  LAYER  NO.  1  ***** 


w 


INITIAL  CXHIDITIONS  FOR  PROJECTILE  LAYER  NO.  1  MATERIAL  ... 


PARTICLE  VELOCITY  .  UP  «  5.000  KM/S 

SHOCK  NAVE  SPEED . US  -  12.080  KM/S 

SHOCK  PRESSURE . PR  >  163.805  GPA  >  1.617  MBAR 

SHOCK  ENERGY . EH  «  .1250E-«-08  JOULES/KG 


SPECIFIC  VOLUME  (AT  REST)  ...  VO  >  .369  CU.CM./GM 

SPECIFIC  VOLUME  (SHOCKED)  ...  VI  »  .216  CU.CM./GM 


*  PARAMETERS  REQUIRED  FOR  CALCULATING  RELEASE  OF  PROJ  LAYER  NO.  1  MATERIAL 

FROM  SHOCKED  STATE  USING  THE  TILLOTSON  EQUATION  OF  STATE: 

MATL  ELASTIC  MODULUS .  B  -  .7102E+11  N/SQ.M. 

MATL  POISSCHI  RATIO . NU  «  .350 

MATL  BULK  MODULUS  . K  -  .78918+11  N/SQ.M. 

MATL  LIN.  COBF.  OF  THERM.  EXP.  ...  ALFA  «  .24008-04  /DEG-C 
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MATX.  8P  HIAT  (SOLID)  .  CPS 

MSXL  SP  HIAT  (LigOID)  . CPL 

NSTL  SMB  MHmOIl  OOBP  (CSL,Z1IP)  ...  QSMO 

MSTL  YZBLD  STRUWra .  SY 

MSTL  ULT  SnumCtTH . SO 

MSYL  BUM  HDN8  MO .  BHN 

MSTL  NBLT  TBMPBRSTOBB . TM 

MSTL  VSPOR  TBMPBRSTORB . TV 

MSTL  HBST  OF  FUSION .  HF 

MSTL  HBST  OF  VSPORIXSTION . HV 

MSTL  IMICPIBMT  MBLT  BMBROY . IMB 

MSTL  IMCIPZBNT  VSPOR  BNBROY  .  IVB 


.235  CSL/GM/DBO-C 
.255  CSL/GM/DBO-C 
2.129,  2.130 
290.000  NPS 
434.000  MPS 

120.000  . 

660.00  DBO-C 
2450.00  DBG-C 
95.00  CSL/GM 
■2450.00  CSL/GM 

.6492B-»-06  J(Kn.B8/KG  # 

.2958B-*‘07  JOOLBS/KG 


SDDITIONSL  PSRSMBTBRS  RBQUZRBD  FOR  CSLCULSTING  MSTBRZSL  RBLBSSB  FROM 
SBOCKBO  8TSTB  USING  TRB  TILL0T80N  BQUATION  OF  STSTBi 
SS  *  .7850B+11  N/SQ.M. 

BB  -  .48318411  N/SQ.M. 

S  -  .5000 

B  -  1.6292 

SLF  >  5.0000 

BBT  >  5.0000 

BOZ  -  .6687B407  JOULBS/KG 

BON  -  1.0000 

BO  -  .6687B407  JOULBS/KG 

BS  «  .2958B407  JOULBS/KG 

HV  •  .1026B408  JOULBS/KG 

BSP  -  .1321B408  JOULBS/KG 

VS  «  .4170  CU.CM. /ON 

BPS  •  .0050 

BMD-STSTB  CSLCULSTKM  RBSULTS  USING  THB  TILLOTSON  BOS  ... 


MSTBRZSL  FIN  SP  VOL  (VF)  .  .494  CU.CM. /GM 

MSTBRZSL  SHOCK  BNBROY . 1250B408  JOULBS/KG 

MSTBRZSL  BNBROY  RBCOVBRBD  . 9992B407  JOULBS/KG 

NSSTB  HBST  GBNBRSTBD . 2508B407  JOULBS/KG 

BNBROY  RBQ,  IMCIPZBNT  NBLT  ...  .6492B406  JOULBS/KG 

BNBROY  RBQ,  COMPLBTB  MBLT . 1047B407  JOULBS/KG 

BXCBSS  BNBROY  SVSZLSBLB . 1461B407  JOULBS/KG 

RBSIDUSL  MSTBRZSL  TBMP  .  2028.692  DBG-C 


PBRCBNT  OF  SHOCKBD  SND  RBLBASBD  MSTBRZSL  ... 


IN  SOLID  STSTB  ...  .00% 

IN  MOLTBN  FORM  ...  100.00% 

IN  VSPOR  FORM . 00% 

PROJBCTILB  LSYBR  DBNSITY  .  2.7120  GN/CU.CM. 

PROJBCTILB  LSYBR  LBNGTH  .  .2540  CM 

PROJBCTILB  LSYBR  DISMBTBR  ....  2.5400  CM 

PROJBCTILB  LSYBR  MSSS  .  3.4904  GMS 


PROJ  LBNGTH  NHBRB  RWSVB  HITS  SNSVB  -  .7009  CM 

LBNGTH  OF  PROJ  LYR  MSTL  SUBJ  TO  S6R  ...  .2540  CM 

MSSS  OF  UNSH  PROJ  LYR  MSTBRZSL .  .0000  GMS 

MSSS  OF  SH  SND  RBL  PROJ  LYR  MSTL  .  3.4904  GMS 

MSSS  OF  S6R  SOLID  LYR  MSTL .  .0000  GMS 

MSSS  OF  S6R  LIQUID  LYR  MSTL .  3.4904  GMS 

MSSS  OF  SSR  VSPOR  LYR  MSTL .  .0000  GMS 

TOTSL  SOLID  LSYBR  MSSS  COMPONBNT  .  .0000  GMS 

TOTSL  NON-SOLID  LSYBR  COMPONBNT  .  3.4904  GMS 


PROJBCTILB  MSTBRZSL  RBLBSSB  CSLCULSTIQNS,  LSYBR  NO.  2  ***** 


ZNITZSL  CONDITIONS  FOR  PROJBCTILB  LSYBR  NO.  2  MSTBRZSL  ... 
PSRTICLB  VBLOCITY  .  UP  -  3.359  KM/S 
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SHOCK  NKVK  SmO . OS  -  9.777  KM/8 

SHOCK  ntSSSOBB . PH  -  257.156  OPA  «  2.539  NBAR 

SHOCK  IMBBCnr  . SB  -  .5642B-I-07  J00LB8/K0 


SPHCXPZC  VOLDMS  (AT  KBST)  VO  •  .128  CD. CM. /ON 

SPHCIPZC  VQLDMB  (SHOCKBD)  ...  VI  -  .084  CO.OI./GN 


PAKANStSKS  KBQDXRSD  PC«  CALCOLATZNO  RBLBASB  OP  PROJ  LAYER  NO.  2  MATERIAL 
PROM  SHOCKED  STATS  USINO  THE  TZLLOT80M  EQUATION  OP  STATSt 


MATL  ELASTIC  NOOOLOS  . 

MATL  POISSON  RATIO  . 

NATL  BULK  NOOOLOS  . 

MATL  LIN.  OOSP.  OP  THERM.  EXP.  ... 

NATL  SP  HEAT  (SOLID)  . 

NATL  SP  BEAT  (LIQUID)  . 

MATL  AMB  M-ORON  0(»P  (CAL, IMP)  ... 

MATL  YIELD  STRSM6TH  . 

MATL  OLT  STRSNOTH  . 

NATL  BRM  HDM8  MO  . 

MATL  MELT  TEMPERATURE  . 

MATL  VAPOR  TEMPERATURE  . 

MATL  BEAT  OP  PU8ION  . 

MATL  BEAT  OP  VAPORIZATKNI . 

MATL  INICPISNT  MELT  ENERGY  . 

MATL  INCIPIENT  VAPOR  ENERGY  . 


S  -  .1999S4-12  N/SQ.M. 

NO  -  .300 

K  -  .1666S+12  N/SQ.N. 

ALFA  «  .1120S-04  /DE6-C 

CPS  «  .110  CAL/GM/DEG-C 

CPL  -  .125  CAL/GM/DSO-C 

GAMO  >  1.553,  1.670 

SY  -  469.000  MPA 

SO  «  745.000  MPA 

BHN  -  290.000 

TM  -  1510.00  DSG-C 

TV  «  3070.00  DSG-C 

BP  «  65.00  CAL/GM 

HV  -  1590.00  CAL/GM 

IMB  >  .6953B<f06  JOOLBS/X6 
XVB  *  .1784B-t-07  JODLBS/KG 


ADDITIONAL  PARAMETERS  REQUIRED  FOR  CALCULATING  MATERIAL  RBLBASB  PROM 
SHOCKED  STATE  USINO  THE  TZLLOTS(N(  EQUATION  OF  STATE t 


AA 

m 

.1635B<<>12  N/SQ.M. 

BB 

m 

.2165B-»-12  N/SQ.M. 

A 

m 

.5000 

B 

m 

1.0525 

ALF 

m 

5.0000 

BBT 

m 

S.OOOO 

EOI 

m 

.1447B-<-08  JOULES/KG 

EON 

m 

1.0000 

EO 

m 

.1447E-«-08  JOULBS/KG 

E8 

m 

.1784B-I-07  JOULBS/KG 

HV 

m 

.6656B407  JOULES/KG 

ESP 

m 

.8439E-I-07  JOULBS/KG 

VS 

m 

.1444  CU.CM. /GN 

BPS 

m 

.0050 

END-STATE  CALCULATION  RESULTS  USING  THE  TILLOTSON  BOS  ... 


MATERIAL  FIN  SP  VOL  (VF)  .  .142  CU.CM. /GM 

MATERIAL  SHOCK  ENERGY . 5642E-«-07  JOULBS/KG 

MATERIAL  ENERGY  RECOVERED . 4211E-t-07  JOULES/KG 

WASTE  HEAT  GENERATED . 1431E<t-07  JODLBS/KG 

ENERGY  RBQ,  INCIPIENT  MELT  ...  .6953B-«-06  JOULBS/KG 

ENERGY  RBQ,  COMPLETE  MELT . 9674B+06  JOULBS/KG 

EXCESS  ENERGY  AVAILABLE . 4639B-t'06  JODLBS/KG 

RESIDUAL  MATERIAL  TEMP .  2396.557  DEG-C 

PERCENT  OF  SHOCKED  AMD  RELEASED  MATERIAL  ... 

IN  SOLID  STATE  ...  .00% 

¥  IN  MOLTEN  FORM  ...  100.00% 

IN  VAPOR  FORM . 00% 

PROJECTILE  LAYER  DENSITY  .  7.8300  GN/CU.CM. 

%  PROJECTILE  LAYER  LENGTH  .  .2540  CM 

PROJECTILE  LAYER  DIAMETER  _  2.5400  CM 

PROJECTILE  LAYER  MASS  .  10.0775  GMS 


PROJ  LENGTH  WHERE  RWAVE  HITS  SWAVB  _  .7009  CM 

LENGTH  OF  PROJ  LYR  MATL  SUBJ  TO  S6R  ...  .2540  CM 
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1AM  or  AIM  PROJ  LYR  lATBRIAL .  .0000  OHS 

MMS  or  M  AMD  MU.  PROJ  LYR  NATL .  10.0775  GN8 

MASS  or  SSR  SOLID  LYR  MATL .  .0000  G»f8 

MAM  AP  SSR  LIQUID  LYR  MATL .  10.0775  OMS 

MMS  or  SSR  VAPOR  LYR  MATL .  .0000  OMS 

TOTAL  SOLID  LAYBR  MASS  OOMPOMIMT .  .0000  OMS  J 

TOTAL  M0M-SCM.X0  LAYBR  OONPOMBMT .  10.0775  OMS 

*****  PROJBCTILB  NATBRIAL  RBLBASB  CALC0LATX0M8,  LAYBR  NO.  3  ***** 


XMXTXAL  OONDXTXCHIS  PCffi  PROJBCTXLB  LAYBR  NO.  3  MATBRXAL  ...  f 


PARTXCLB  VBLOCXTY . OP  »  2.559  KM/S 

SHOCK  MAVB  SPBBD . 08  -  7.315  XM/8 

SHOCK  PRBS80RB . PH  >  358.834  OPA  -  3.542  MBAR 

SHOCK  BMBROY . BH  -  .3274B<i>07  JOOLBS/XO 


SPBCXrXC  VOLOMB  (AT  RB8T)  ...  VO  «  .052  C0.CM./6M 

SPBCXPXC  VOLOMB  (SHOCKBO)  ...  VI  -  .034  CO.CM./OM 


PARAMBTBRS  RBQOXRBD  rOR  CALCOLATXNO  RBLBASB  OP  PROJ  LAYBR  NO.  3  MATBRXAL 
ntOM  SHOCKBO  8TATK  OSXNO  THB  TXLLOTSM  BQOATXOM  OP  STATBt 


MATL  BLA8TXC  MMULOS . 

MATL  POXSSON  RATXO . 

MATL  BULK  MMOLUS . 

MATL  LXM.  OOBP.  OP  TBBRM.  BXP.  ... 

MATL  SP  BBAT  (SOLXD)  . 

MATL  SP  BBAT  (LXQOXD)  . 

MATL  AMB  M-GRON  OOBP  (CAL,XMP)  ... 

MATL  YXBLD  STRBNGTH  . 

MATL  OLT  STRBNGTH  . 

MATL  BRM  HONS  NO  . 

MATL  MBLT  TBMPBRATORB  . 

MATL  VAPOR  TBMPBRATORB  . 

NATL  HBAT  OP  POSXOM  . 

MATL  BBAT  OP  VAPORXSATXON  . 

MATL  XNXCPXBMT  MBLT  BMBRGY  . 

MATL  XNCXPXBNT  VAPOR  BNBRGY  . 


B 

.40688*12  N/SQ.M. 

NO 

.300 

K 

.3390B*12  N/SQ.M. 

ALPA 

.4000B-05  /OBG-C 

CPS 

.035  CAL/GM/DBG-C 

CPL 

.046  CAL/GM/DBG-C 

GAMO 

1.448,  1.480 

SY 

1379.000  MPA 

SO 

1517.000  MPA 

BHN 

400.000 

TN 

3410.00  DBG-C 

TV 

5900.00  DBG-C 

HP 

53.00  CAL/GM 

HV 

1054.00  CAL/GM 

XMB 

.4996B*06  JOULBS/KG 

XVB 

.1201B*07  JOULBS/KG 

t 


ADOXTXONAL  PARAMBTBRS  RBQOXRBD  POR  C^LCVLATZKO  MATBRXAL  RBLBASB  PROM 
SHOCKBO  STATB  OSXNG  THB  TXLLOTSM  BQOATXON  OP  STATBt 


AA 

■ 

.33028*12  N/SQ.M. 

SB 

m 

.2476B*12  N/SQ.M. 

A 

m 

.5000 

B 

m 

.9481 

ALP 

m 

10.0000 

BBT 

m 

10.0000 

BOX 

m 

.2574B*07  JOULBS/KG 

BOM 

m 

1.0000 

BO 

m 

.2574B*07  JOULBS/KG 

BS 

m 

.12018*07  JOULBS/KG 

HV 

m 

.44128*07  JOULBS/KG 

BSP 

m 

.56138*07  JOULBS/KG 

VS 

m 

.0590  Cl  CM./GM 

BPS 

m 

.0050 

BMD-STATB  CALCOLATXON  RBSOLTS 
MATBRXAL  PXN  SP  VOL  (VP)  .... 

MATBRXAL  SHOCK  BNBRGY  . 

MATBRXAL  BMBRGY  RBCOVBRBD  ... 

MASTB  BBAT  GBMBRATBD  . 

BMBROY  RBQ,  XNCXPXBNT  MBLT  .. 
BNBRGY  RBQ,  OOMPLBTB  MBLT  ... 
BNBRGY  AVAXLABLB  POR  MBLT  . . . 
RB8X00AL  MATBRXAL  TBMP  . 


OSXNG  THB  TXLLOTSON  BOS  ... 

.055  CU.CM. /GM 
.3274B-f07  JOULBS/KG 
.2564B-f07  JOOLBS/K6 
.7095B-«'06  JOULBS/KG 
.4996B+06  JOULBS/KG 
.7215B-I-06  JOULBS/KG 
.2099B-«>06  JOULBS/KG 
3410.000  DBG>C 
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WWCIOT  or  SBOaCID  and  UUMSKD  MATntZAL  ... 
XH  aOLXO  raOB  ...  S.40% 

ZR  II0I.TRM  rom  ...  94.60% 

ZM  VAPOR  PORN . 00% 


PROJICTZZX  ZAZRR  DSMSZTZ .  19.1700  ON/CO. CM. 

PRDJICTZXJI  ZAXIR  LIMOTH . 2S40  CM 

PROJBCTZXR  Umn  DZANITSR  ....  2.S400  CM 

PROJICrZXR  LAXIR  MASS .  24.6725  QMS 

PROJ  IMmtB  NHIRI  RNAVI  HITS  SNAVI  ....  .7009  CM 

ZJMOTH  or  PROJ  LTR  MATL  SOBJ  TO  S6R  ...  .1929  CM 

MASS  or  UMSH  PROJ  LTR  MATIRZAL .  5.9316  QMS 

MASS  OP  SH  AMD  RHL  PROJ  LTR  MATL .  18.7409  QMS 

MASS  OP  S6R  SOLZD  LTR  MATL .  1.0128  GMS 

MASS  OP  S6R  LZQOZO  LTR  MATL .  17.7281  OMS 

MASS  OP  S6R  VAPOR  LTR  MATL .  .0000  OMS 

TOTAL  SOLZD  LAYKR  MASS  OOMPOMKMT .  6.9444  GMS 

TOTAL  MOH'SOLZD  LATIR  C0MP0M8MT .  17.7281  OMS 


MASS  DZSTRZBOTZOM  SDNMART  ... 

PROJRCTZL8  LAYKR  MO.  1  ...  SOLID  ....  .00  OMS 

OMSH  ....  .00  OMS 

SMR .  .00  GMS 

S6R .  .00  OMS 

LZQOID  ...  3.49  OMS 

VAPC»  ....  .00  OMS 

PROJXCTILX  LAYKR  HO.  2  ...  SOLZD  ....  .00  OMS 

UMSH  ....  .00  OMS 

SMR .  .00  OMS 

S6R .  .00  GMS 

LZgUZD  ...  10.08  GMS 

VAPOR  ....  .00  OMS 

PROJKCTZLK  LAYKR  MO.  3  ...  SOLID  ....  6.94  GMS 

UMSH  ....  5.93  GMS 

SMR .  .00  GMS 

S6R .  1.01  OMS 

LZQUZD  ...  17.73  OMS 

VAPOR  ....  .00  GMS 

TAROKT  MATKRIAL  .  SOLID  ....  15.15  GMS 

PRAO  ....  15.15  GMS 

S6R .  .00  OMS 

LIQUID  ...  4.36  GMS 

VAPOR  ....  .00  GMS 


SOMK  UMSHOCRKD  PROJ  MATL  RKMAINS  ... 

TOT  MPROJfUMSH  ...  5.932  GMS  (-  .155NPROJ) 

AVO  VP,UMSH  .  11.053  KM/S  (>1.105V) 

ALL  OTHXR  SOLID  PROJKCTZLK  MATKRIAL  (IP  ANY)  IS  LIXXLY  TO  BK  PRAGMKNTBD. 
AMY  SOLID  TAROKT  MATKRIAL  RKMAINING  IS  ALSO  LIKKLY  TO  BK  PRAGMKNTKD. 

THK  DKBRIS  CLOUD  CONSISTS  OP  SHOCXKD  AND  RKLKASKD  PROJBCTILK  MATKRIAL  AMD 
ALL  KJKCTKD  TAROKT  MATKRIAL. 

DKBRIS  CLOUD  VBLOCITY  SUMMARY  ... 

DKB  CLD  CKNTBR-OP-MASS  VKL  (VCON)  ...  6.114  KM/SBC  («  .611V) 

DKB  CLD  LKADING  KDOK  VKL  (VL)  .  10.652  KM/SKC  (-1.065V) 

DKB  CLD  KXPANSION  VKL  (VKXP)  .  4.539  KM/SKC  («  .454V) 


DKBRIS  CLOUD  HALP-ANOLK 


36.589  DBG 


DISTRIBUTION  UST 
(WL-TR-94-7039) 


E^ia  AFB  Olficef: 

WL/MNSA  10 

WL/MNOI  (STINFO  Facility)  1 
WL/CA-N  1 

HQ  Defenae  Nudear  Agency/SPSP 
Attn:  MrTony  Fredeiidcaon 
6601  Td^tnph  Rd 
Alexandria  VA  22310^396 
2 

Space  and  Miaaile  Systems  Center/MG  WE 
Attn:  Mr  Mdidi  Kavaiy 
P.O.  Box  92960 
Los  Angeles  CA  90009-2960 
1 

BMIXVAQT 
The  Penta^n 

Attau  Lt  Col  Charles  Martin 
Washington  DC  20301-7000 
1 

USA  Space  A  Strategic  Defense  Command 
CSSD-SL-L 

Attau  Dr  Robert  Becker 
PO  BOX  1500 
HnmsviUeAL  35807-3801 
1 


4 

♦ 


DASIAC 

2560  Huntington  Ave,  Suite  500 
Alexandria  VA  22303-1490 
1 


DTIC/DDAC 
Cameron  Station 
Alexandria  VA  22304-6145 


148 


